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O r Abstract. Neutrinoless double beta decay, which is a very old and yet elusive process, 

is reviewed. Its observation will signal that lepton number is not conserved and the 
neutrinos are Majorana particles. More importantly it is our best hope for determining 
the absolute neutrino mass scale at the level of a few tens of meV. To achieve the last 
goal certain hurdles have to be overcome involving particle, nuclear and experimental 
physics. 

Nuclear physics is important for extracting the useful information from the data. 
One must accurately evaluate the relevant nuclear matrix elements, a formidable 
task. To this end, we review the sophisticated nuclear structure approaches recently 
been developed, which give confidence that the needed nuclear matrix elements can 
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. be reliably calculated employing different methods: a) the various versions of the 

Quasiparticle Random Phase Approximations, b) the interacting boson model, c) the 
energy density functional method and d) the large basis Interacting Shell Model. It 
is encouraging that, for the light neutrino mass term at least, these vastly different 
approaches now give comparable results. 
& ■ From an experimental point of view it is challenging, since the life times are long 

and one has to fight against formidable backgrounds. One needs large isotopically 
enriched sources and detectors with high energy resolution, low thresholds and very 
low background. 

If a signal is found, it will be a tremendous accomplishment. Then, of course, the real 
task is going to be the extraction of the neutrino mass from the observations. This 
is not trivial, since current particle models predict the presence of many mechanisms 
other than the neutrino mass, which may contribute or even dominate this process. 
We will, in particular, consider the following processes: 

(i) The neutrino induced, but neutrino mass independent contribution. 

(ii) Heavy left and/or right handed neutrino mass contributions. 

(iii) Intermediate scalars (doubly charged etc). 

(iv) Supersymmctric (SUSY) contributions. 

We will show that it is possible to disentangle the various mechanisms and 
unambiguously extract the important neutrino mass scale, if all the signatures of the 
reaction are searched in a sufficient number of nuclear isotopes. 
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1. A brief history of double beta decay 

A brief history of the double-beta decay is presented. 

1.1. The early period 

Double beta decay (DBD), namely the two-neutrino double-beta decay (2^/9/3-decay) 

(A, Z) (A, Z + 2) + e~ + e~ + V e + V e , (1) 

was first considered in publication [1] of Maria Goeppert-Mayer in 1935. It was Eugene 
Wigner, who suggested this problem to the author of [1] about one year after the 
Fermi weak interaction theory appeared. In the work of Maria Goeppert-Mayer [T] 
an expression for the 2z//3/3-decay rate was derived and a half-life of 10 17 years was 
estimated by assuming a Q-value of about 10 MeV. 

Two years later (1937) Ettore Majorana formulated a new theory of neutrinos, 
whereby the neutrino u and the antineutrino V are indistinguishable, and suggested 
antineutrino induced /3~-decay for experimental verification of this hypothesis [2j. Giulio 
Racah was the first, who proposed testing Majorana's theory with real neutrinos by chain 
of reactions 

(A, Z) ->(A,Z + l) + e- + v, v + (A', Z') ->• (A 1 , Z' + 1) + e~, (2) 

which is allowed in the case of the Majorana neutrino and forbidden in the case of Dirac 
neutrino j3]. In 1939, Wolfgang Furry considered for the first time neutrinoless double 
beta decay (O^/3/3-decay), 

(A,Z) -> (A,Z + 2) +e~ +e" (3) 

a Racah chain of reactions with virtual neutrinos ((A, Z + 1) = (A', Z')) pQ. Here A, A' 
are the nuclear mass numbers and Z, Z' the charges of the nuclei involved. The available 
energy A is equal to the Q-value of the reaction, i.e. the mass difference of the ground 
states of the two atoms involved. 

In 1952 Henry Primakoff [5] calculated the electron-electron angular correlations and 
electron energy spectra for both the 2z//3/9-decay and the O^/S/S-decay, producing a useful 
tool for distinguishing between the two processes. 

At that time nothing was known about the chirality suppression of the 0vj3{3- 
decay. It was believed that, due to a considerable phase-space advantage, the Ou/3/3- 
decay mode dominates the double beta decay rate. Starting 1950 this phenomenon 
was exploited in early geochemical, radiochemical and counter experiment ijjl It was 
found that the measured lower limit on the /3/3-decay half-life far exceeds the values 
expected for this process, TfY 2 ~ 10 12 — 10 15 years. In 1955 the Raymond Davis 
experiment [7], which searched for the antineutrinos from reactor via nuclear reaction 
V e + 37 CI — > 37 Ar + e~, produced a zero result. The above experiments were interpreted 

| For more detailed historical review of the experimental activities in the field of double-beta decay 
we recommend reader a recent review on this subject by A. Barabash [B]. 
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as proof that the neutrino was not a Majorana particle, but a Dirac particle. This 
prompted the introduction of the lepton number to distinguish the neutrino from its 
antiparticle. The assumption of lepton number conservation allows the 2z//3/3-decay but 
forbids the Oz//3/3-decay, in which lepton number is changed by two units. 

In 1949 Fireman reported the first observation of the /3/3-decay of 124 Sn in a 
laboratory experiment [8J, but he disclaimed it later |9j. The first geochemical 
observation of the /3/3-decay, with an estimated half-life T 1 / 2 ( 130 7 1 e) = 1.4 x 10 21 years, 
was announced by Ingram and Reynolds in 1950 (TO]. Extensive studies have been made 
by Gentner and Kirsten [Til H2] and others [131 H3] on such rare-gass isotopes as 82 Kr, 
128 Xe, and 130 Xe, which are /3/3-decay products of 82 Se 128 Te, and 130 Te, respectively, 
obtaining half lives around 10 21 y for 130 Te. 

1.2. The period of scepticism 

Shortly after its theoretical formulation by Lee and Yang, parity violation in the weak 
interaction was established by two epochal experiments. In 1957 Wu et al. discovered 
the asymmetry in the angular distribution of the /3-particles emitted relative to the spin 
orientation of the parent nucleus 60 Co. A year later Goldhaber et al [15] discovered 
that the neutrinos are polarized and left handed by measuring the polarization of a 
photon, moving back to back with the neutrino, produced by the de-excitation of a 
152 Eu* nucleus after K-capture. In 1958 the seemingly confused situation was simplified 
in the form of the vector-axial vector (V -A) theory of weak interactions describing 
maximal parity violation in agreement with available data. In order to account for the 
chiral symmetry breaking of the weak interaction only left handed fermions participate 
and the mediating particles must be vectors of spin 1, which are left handed in the sense 
that they couple only to left handed fermions. 

The maximal parity violation is easily realized in the lepton sector by using the 
two-component theory of a massless neutrino, proposed in 1957 by L. Landau, T.D. 
Lee, C.N. Yang and A. Salam (This idea was first developed by H. Weyl in 1929, but 
it was rejected by Pauli in 1933 on the grounds that it violates parity). In this theory, 
neutrinos are left handed and antineutrinos are right-handed, leading automatically to 
the V - A couplings. 

With the discovery of parity violation, it became apparent that the Majorana/Dirac 
character of the electron neutrino was still in question. The particles that participate 
in the O^/3/3-decay reaction at nucleon level are right-handed antineutrino V e and left 
handed neutrino v e : 

n ->■ p + e" +vf, + n ->■ p + e~ . (4) 

Thus even if the neutrino is a (massless) Majorana particle, since the first neutrino has 
the wrong helicity for absorption by a neutron, the absence of the Oi//3/3-decay implies 
neither a Dirac electron neutrino nor a conserved lepton number. 

The requirement that both lepton number conservation and the 75 invariance of 
the weak current had to be violated, in order the O^/3/3-decay to occur, discouraged 
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experimental searches. 

1.3. The period of Grand Unified Theories 

The maximal violation of parity (and of charge-conjugation) symmetry is accommodated 
in the Standard Model (SM), which describes jointly weak and electromagnetic 
interactions. This model was developed largely upon the empirical observations 
of nuclear beta decay during the latter half of the past century. Despite the 
phenomenological success of the SM, the fundamental origin of parity violation has 
not been understood. In spite of the fact that the SM represents the simplest and the 
most economical theory, it has not been considered as the ultimate theory of nature. 
It was assumed that, most likely, it describes a low energy approximation to a more 
fundamental theory. 

With the development of modern gauge theories during the last quarter of the 
previous century, perceptions began to change. In the SM it became apparent that the 
assumption of lepton number conservation led to the neutrino being strictly massless, 
thus preserving the 75-invariance of the weak current. With the development of Grand 
Unified Theories (GUT's) of the electroweak and strong interactions, the prejudice has 
grown that lepton number conservation was the result of a global symmetry not of a 
gauge symmetry and had to be broken at some level. In other words modern GUT's and 
supersymmetric (SUSY) extensions of the SM suppose that such conservation laws of 
the SM may be violated to some small degree. The lepton number may only appear to be 
conserved at low energies because of the large grand unified mass scale A.qut governing 
its breaking. Within the proposed see-saw mechanism one expects the neutrino to 
acquire a small Majorana mass of a size ~ (light mass) 2 / Aqut, where "light mass" is 
typically that of a quark or charged lepton. The considerations of a sensitivity of the 
0^/3/9-decay experiments to a neutrino mass m u ~ 1 eV became the genesis of a new 
interest to double beta decay. Thus the interest in O^/3/3-decay was resurrected through 
the pioneering work of Kotani and his group [16], which brought it again to the attention 
of the nuclear physics community. 

Neutrino masses require either the existence of right-handed neutrinos or require 
violation of the lepton number (LN) so that Majorana masses are possible. So, one is 
forced to go beyond the minimal models again, whereby LF and/or LN violation can be 
allowed in the theory. A good candidate for such a theory is the left-right symmetric 
model of Grand Unification inaugurated by Salam, Pati, Mohapatra and Senjanovic 
[T7J HU [19] and especially models based on SO(10), which have first been proposed by 
Fritzsch and Minkowski [20] , with their supersymmetric versions [2U |22j [23] . The left- 
right symmetric models, representing generalization of the SU(2)l eg) U(l) SM, predict 
not only that the neutrino is a Majorana particle, that means it is up to a phase 
identical with its antiparticle, but automatically predict the neutrino has a mass and a 
weak right-handed interaction. 

In the left-right symmetric models the LN conservation is broken by the presence of 
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the Majorana neutrino mass. The LN violation is also inherent in those SUSY theories 
whereby R-parity, defined as R p = (— i) 3 - B + L + 25 ' ) with S, B, and L being the spin, 
baryon and lepton number, respectively is not a conserved quantity anymore. 

The O^/3/3-decay, which involves the emission of two electrons and no neutrinos, 
has been found as a powerful tool to study the LN conservation. Schechter and Valle 
proved that, if the Oz//3/3-decay takes place, regardless of the mechanism causing it, the 
neutrinos are Majorana particles with non-zero mass [2U [25] . It was recognized that 
the GUT's and R-parity violating SUSY models offer a plethora of the Oz//3/3-decay 
mechanisms triggered by exchange of neutrinos, neutralinos, gluinos, leptoquarks, etc. 
[2611271 [28]. 

The experimental effort concentrated on high Qpp isotopes, in particular on 48 Ca, 
76 Ge, 82 Se, 96 Zr, 100 Mo, 116 Cd, 130 Te, 136 Xe and 150 Nd j^JEHlET]. In 1987 the first actual 
laboratory observation of the two neutrino double beta decay (2z//3/3-decay) was done 
for 82 Se by M. Moe and collaborators [32], who used a time projection chamber. Within 
the next few years, experiments employing counters were able to detect 2z//3/3-decay of 
many nuclei. In addition, the experiments searching for the signal of the Oz//3/3-decay 
pushed by many orders of magnitude the experimental lower limits for the 0i//9/9-decay 
half-life of different nuclei. 

1.4- The period of massive neutrinos - the current period 

Various early measurements of neutrinos produced in the sun, in the atmosphere, and by 
accelerators suggested that neutrinos might oscillate from one "flavor" (electron, muon, 
and tau) to another, expected as a consequence of non-zero neutrino mass. Non-zero 
neutrino mass can be accommodated by fairly straightforward extensions of the SM of 
particle physics. Starting 1998 we have a convincing evidence about the existence of 
neutrino masses due to SuperKamiokande [33], SNO, [31] KamLAND [35] and other 
experiments. 

Thus neutrino oscillations have supplied additional information in constructing 
Grand Unified Theories of physics. It also has provided additional input for cosmologists 
and opened new perspectives for observation of the O^/9/3-decay. 

So far the 2i/ / 3 / 3-decay has been recorded for eleven nuclei ( 48 Ca, 76 Ge, 82 Se, 
96 Zr, 100 Mo, 116 Cd, 128 Te, 130 Te, 150 Nd, 136 Xe, 238 U) [291 EH HEP- In addition, the 
2i//3/3-decay of 100 Mo and 150 Nd to + excited state of the daughter nucleus has been 
observed and the two-neutrino double electron capture process in 130 Ba has been 
recorded. Experiments studying 2v (3/3- decay are presently approaching a qualitatively 
new level, when high-precision measurements are performed not only for half-lives but 
also for all other observables of the process. As a result, a trend is emerging towards 
thorough investigation of all aspects of 2i//3/3-decay, and this will furnish very important 
information about the values of nuclear matrix elements, the parameters of various 
theoretical models, and so on. In this connection, one may expect advances in the 
calculation of nuclear matrix elements and in the understanding of the nuclear-physics 
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aspects of double beta decay. 

Neutrinoless double beta decay has not yet been confirmed. The strongest limits 
on the half-life of the O^/3/9-decay were set in Heidelberg-Moscow [36], NEM03 [37] 138]. 
CUORICINO [39] and KamLAND-Zen 00] experiments: 

Ty 2 ( 76 Ge) > 1.9 x 10 25 y, T?fc{ 100 Mo) > 1.0 x 10 24 y, 

T°y 2 ( 130 Te) > 3.0 x 10 24 y, T^^Xe) > 5.7 x 10 24 y. (5) 

There is also a claim for an observation of the Oi//3/3-decay of Ge with half-life 
Ty 2 = 2.23t[}; 44 x 10 25 years (3D [32]. One of the the goals of the upcoming GERDA 
experiment |43j is to put this claim to a test by improving the sensitivity limit of the 
detection by more than an order of magnitude. The next generation experiments, which 
will be using several other candidate nuclei, will eventually be able to achieve this goal 
as well [31] . 

1.5. The period of Majorana neutrinos? 

There is a hope that the period of Majorana neutrinos is not far. This period should 
start by a direct and undoubtable observation of the 0^/3/9-decay. It would establish that 
neutrinos are Majorana particles, and a measurement of the decay rate, when combined 
with neutrino oscillation data and a reliable calculation of nuclear matrix elements, 
would yield insight into all three neutrino mass eigenstates. 

2. An overview 

The question of neutrino masses and mixing is one of the most important issues of 
modern particle physics. It has already been discussed in a number of excellent reviews 
[331 US I3EJ H7J HE] and its relevance to the Oz//3/3-decay will also be briefly discussed in 
this report. 

Today, seventy five years later, Oi//3/3-decay ([3]), continues to be one of the most 
interesting processes. The experimental status and prospects regarding this process will 
be reviewed in section El The corresponding non exotic 2vj3j3 (CQ) has been observed in 
many systems, see section [HJ 

If the neutrinos are Majorana particles other related processes in which the charge 
of the nucleus is decreased by two units may also occur, if they happen to be allowed 
by energy and angular momentum conservation laws, e.g. 

(A, Z) — y (A, Z — 2) + e + + e + (0z/ positron emission). (6) 

Here the available energy is A = Q — Am e c 2 , i.e., everything else being equal, it is 
somewhat kinematically disfavored compared to the usual two electron emission for 
which the available energy A is equal to the Q value. 
Electron positron conversion: 

(A, Z) + — > (A, Z — 2) + e + (Ou electron positron conversion). (7) 
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Nuclear Femto Laboratory 
for (3 (3 experiments 




Figure 1: Schematic diagrams of f3f3 decays in nuclear femto(10 _15 m) laboratories, where 
single /9-decay is not allowed and neutrinoless DBD is much enhanced. Left hand side: 
DBD decay scheme. Right hand side: Neutrinoless DBD with the Majorana v exchange 
between two nucleons [49J. 

A = Q — 2m e c 2 — 6b, where is a binding energy of the absorbed atomic electron. 
The resonant neutrinoless double electron capture (Oz/ECEC), 

(A, Z)+e^+e^ — > (A, Z—2)**, (resonant Ov double electron capture ), (8) 

was first considered very long time ago by Winter [50]. It is always allowed, whenever 
(J7J) is. This reaction was expounded in more detail later [5TJ [52] as a two step process: 
In the first step the two neutral atoms, (A,Z) and the excited (A,Z-2), get admixed 
via the lepton number violating interaction. In the second step the (A,Z-2) atom and, 
possibly, the nucleus de-excite. The available energy is A = Q — B 2 h, B 2 h being the 
energy of two electron holes in the atomic shells of daughter nucleus. 
Decays to excited states are in some cases possible. Then in addition to x rays one has 
various decay modes with emission of a single 7, a pair of 7's, internal electron-positron 
pair formation and emission of electron by internal conversion [53], i.e. 

(A, Z) + e 6 - + e 6 - -> (A, Z - 2) + X, X = 7, 2 7 , e+e", e~ nt (9) 

The life time expected was very long, since the above mixing amplitude was tiny 
compared to the energy difference of the two atoms involved. It has recently, however, 
been gaining in importance [511 ES] after ion Penning traps [56] made it possible to 
accurately determine the Q values, which gave rise to the the presence of resonances. 
This, in turn, could lead to an increase of the width by many orders of magnitude, see 
section [T21 for details. 

Another lepton number violating process, not hindered by energy conservation in 
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any nuclear system, involves the neutrinoless bound muon capture {571 158] , 

(A, Z) + Hi — » (A, Z — 2) + e + (Ou muon-positron conversion). (10) 

The best experimental limit on the muon to positron conversion branching ratio has been 
established at PSI [59] for the 48 Ti nuclear target. The muonic analogue of neutrinoless 
double beta decay [601 El] , 

(A, Z) + Hb ~ * (A Z — 2) + fi + (0^ muon-muon conversion), (11) 

has never been searched for experimentally. For bound muon in atom energy 
conservation of the process is very restrictive. It is satisfied only for three isobars, 
namely 44 Ti, 72 Se and 82 Sr. 

The above processes are expected to occur whenever one has lepton number 
violating interactions. Lepton number, being a global quantity, is not sacred, but it is 
expected to be broken at some level. In short, these processes pop up almost everywhere, 
in every theory. On the other hand since, if there exist lepton number violating 
interactions, the neutrinos have to be Majorana particles, all the above processes can, in 
principle, decide whether or not the neutrino is a Majorana particle, i.e. it coincides with 
its own antiparticle. This is true even if these processes are induced not by intermediate 
neutrinos but by other mechanisms as we will see below. 

Neutrinoless double beta decay (Eq. (E}) seems to be the most likely to yield the 
information [HI E21 E21 EH [651 ES E7J [68] we are after. For this reason we will focus 
our discussion on this reaction, but we will pay some attention to resonant neutrinoless 
double electron capture, which has recently been revived [5T1 [521 EH E5J Ell HOI ITT] , 
since its observation seems to be a realistic possibility [54l 155] . We will only peripherally 
discuss the other less interesting processes [5T] . 

From a nuclear physics [651 ESI [721 E2 EH EH1 EH] point of view, calculating the 
relevant nuclear matrix elements is indeed a challenge. First almost all nuclei, which 
can undergo double beta decay, are far from closed shells and some of them are even 
deformed. One thus faces a formidable task. Second the nuclear matrix elements are 
small compared to a canonical value, like the one associated with the matrix element 
to the (energy non allowed) double Gamow- Teller resonance or a small fraction of some 
appropriate sum rule. Thus, effects which are normally negligible, become important 
here. Third in many models the dominant mechanism for O^/3/3-decay does not involve 
intermediate light neutrinos, but very heavy particles. Thus one must be able to cope 
with the short distance behavior of the relevant operators and wave functions (see section 
HO] for details). 

From the experimental point of view it is also very challenging to measure the 
slowest perhaps process accessible to observation. Especially, if one realizes that even, 
if one obtains only lower bounds on the life time for this Oz//5/3-decay, the extracted 
limits on the theoretical model parameters may be comparable, if not better, and 
complementary to those extracted from the most ambitious accelerator experiments. 

The recent discovery of neutrino oscillations J77J [7H1 [791 EQ] have given the first 
evidence of physics beyond the Standard Model (SM) and in particular they indicate 
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that the neutrinos are massive particles. They were able to show that the neutrinos are 
admixed, determined two of the mixing angles and set a stringent limit on the third 
(for a global analysis see, e.g., [81]). Furthermore they determined one square mass 
difference and the absolute value of the other. Neutrino oscillations, however, cannot 
determine: 

• Whether the neutrinos are Majorana or Dirac particles. 

It is obviously important to proceed further and decide on this important issue. 
Neutrinoless double beta decay can achieve this, even if, as we have mentioned, 
there might be processes that dominate over the conventional intermediate neutrino 
mechanism of O^/3/3-decay. It has been known that whatever the lepton number 
violating process is, which gives rise to Oz//3/3-decay, it can be used to generate a 
Majorana mass for the neutrino [23]. This mechanism, however, may not be the 
dominant mechanism for generating the neutrino mass |82j . 

• The scale of the neutrino masses. 

These experiments can measure only mass squared differences. 
This task can be accomplished by astrophysical observations or via other 
experiments involving low energy weak decays, like triton decay or electron capture, 
or the Oz//3/3-decay. It seems that for a neutrino mass in meV, (10 -3 eV), region, the 
best process to achieve this is the Oi//3/3-decay. The extraction of neutrino masses 
from such observations will be discussed in detail and compared with each other 
later (see section HJ). 

• The neutrino hierarchy 

They cannot at present decide which scenario is realized in nature, namely the 
degenerate, the normal hierarchy or the inverted hierarchy. They may be able to 
distinguish between the two hierarchies in the future by observing the wrong sign 
muons in neutrino factories [831 E3] ■ 

For details on such issues see a recent review [85J. 

The study of the Oz//3/5-decay is further stimulated by the development of GUT's and 
supersymmetric models (SUSY) representing extensions of the SU(2) L <g) U(l) standard 
model. The GUT's and SUSY offer a variety of mechanisms which allow the O^/S/S-decay 
to occur [86] . 

The best known mechanism leading to Oz//3/3-decay is via the exchange of a Majorana 
neutrino between the two decaying neutrons [331 E21 E31 E31 EHJ [66J, [87] . Nuclear physics 
allows us to study the light (m v <C m e ) and heavy (m„ 3> m p ) neutrino components 
separately. In the presence of only left handed currents and for the light intermediate 
neutrino components, the obtained amplitude is proportional to a suitable average 
neutrino mass, which vanishes in the limit in which the neutrinos become Dirac particles. 
On the other hand in the case of heavy Majorana neutrino components the amplitude is 
proportional to the average of the inverse of the neutrino mass, i.e. it is again suppressed. 
In the presence of right handed currents one can have a contribution similar to the one 
above for heavy neutrinos but involving a different (larger) average inverse mass with 
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some additional suppression due to the the fact the right handed gauge boson, if it 
exists, is heavier than the usual left handed one. 

In the presence of right handed currents it is also possible to have interference 
between the leptonic left and right currents, ji — ]r interference. In this case the 
amplitude in momentum space becomes proportional to the 4-momentum of the neutrino 
and, as a result, only the light neutrino components become important. One now has two 
possibilities. First the two hadronic currents have a chirality structure of the same kind, 
i.e. Jl — Jr- Then one can extract from the data a dimensionless parameter A, which 
is proportional to the square of the ratio of the masses of the L and R gauge bosons, 
K = (mL/mji) 2 . Second the two hadronic currents are left handed, which can happen 
via the mixing e of the two bosons. The relevant lepton number violating parameter rj 
is now proportional to this mixing e. Both of these parameters, however, also involve 
the neutrino mixing. They are, in a way, proportional to the mixing between the light 
and heavy neutrinos. 

In gauge theories one has, of course, many more possibilities. Exotic intermediate 
scalars may mediate Oz//3/3-decay [44J. These are not favored in current gauge theories 
and are not going to be further discussed. In superstring inspired models one may have 
singlet fermions in addition to the usual right handed neutrinos. Not much progress has 
been made on the phenomenological side of these models and they are not going to be 
discussed further. 

In recent years supersymmetric models are taken seriously and semirealistic 
calculations are taking place. In standard calculations one invokes universality at 
the GUT's scale, employing a set of five independent parameters, and uses the 
renormalization group equation to obtain all parameters (couplings and particle masses) 
at low energies. Hence, since such parameters are in principle calculable in terms of the 
five input parameters, one can use experimental data to constrain the input parameters. 
One, then, can use the Oi//3/3-decay experiments to constrain the R-parity violating 
couplings, which cannot be specified by the theory [26l [881 EHl EOl EH [92l [93l [9^1 [95] . 
Recent review articles [3TJ [651 ESI EE] also give a detailed account of some of the latest 
developments in this field. 

From the above discussion it is clear that one has to consider the case of heavy 
intermediate particles. One thus has to tackle problems related to the very short ranged 
operators in the presence of the nuclear repulsive core. If the interacting nucleons 
are point-like one gets negligible contributions. We know, however that the nucleons 
are not point like, but that they have structure described by a quark bag with a size 
that can be determined experimentally. It can also be accounted for by a form factor, 
which can be calculated in the quark model or parametrized in a dipole shape with a 
parameter determined by experiment. This approach, first considered by Vergados |96j, 
has now been adopted by almost everybody. The resulting effective operator has a range 
somewhat less than the inverse of the proton mass (see sect. 4 below). 

Another approach in handling this problem consists of considering particles other 
than the nucleons present in the nuclear soup. For + — > + the most important such 
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particles are the pions. One thus may consider the double beta decay of pions in flight 
between nucleons, like 

7r~ — > 7r + + e~ + e~ , n — > p + n + + e~~ + e~ . (12) 

Recognition of such contribution first appeared as a remark by the genius of Pontecorvo 
[97] in the famous paper in which he suggested that the ratio of the lifetimes of the 
128 Te and 130 Te isotopes, which merely differ by two neutrons, is essentially independent 
of nuclear physics. He did not perform any estimates of such a contribution. Such 
estimates and calculations were first performed by Vergados [98] in the case of heavy 
intermediate neutrinos, .i.e. vector and axial vector currents. It was found that it yields 
results of the same order as the nucleon mode with the above recipe for treating the 
short range behavior. It was revived by the Tuebingen group [261 [271 1281 193] in the 
context of R-parity violating interactions, i.e. scalar, pseudoscalar and tensor currents 
arising out of neutralino and gluino exchange, and it was found to dominate. 

In a yet another approach one may estimate the presence of six quark clusters in 
the nucleus. Then, since the change of charge takes place in the same hadron there is 
no suppression due to the short nature of the operator, even if it is a 5-function. One 
only needs a reliable method for estimating the probability of finding these clusters in 
a nucleus [99] . 

All the above approaches seem reasonable and lead to quite similar results. The 
matrix elements obtained are not severely suppressed. This gives us a great degree of 
confidence that the resulting matrix elements are sufficiently reliable, allowing double 
beta decay to probe very important physics. 

The other recent development is the better description of nucleon current by 
including momentum dependent terms, such as the modification of the axial current 
due to PCAC and the inclusion of the weak magnetism terms. These contributions 
have been considered previously j87J 1100] , but only in connection with the extraction 
of the rj parameter mentioned above. Indeed these terms were very important in this 
case since they compete with the p-wave lepton wave function, which, with the usual 
currents, provides the lowest non vanishing contribution. Since in the mass term only 
s-wave lepton wave functions are relevant such terms have hitherto been neglected. 

It was recently found [101] . however, that for light neutrinos the inclusion of 
these momentum dependent terms reduces the nuclear matrix element by about 25%, 
independently of the nuclear model employed. On the other hand for heavy neutrinos 
the effect can be larger and it depends on the nuclear wave functions. The reason 
for expecting them to be relevant is that the average momentum (q) of the exchanged 
neutrino is expected to be large [102] . In the case of a light intermediate neutrino 
the mean nucleon-nucleon separation is about 2 fm which implies that the average 
momentum (q) is about 100 MeV/c. In the case of a heavy neutrino exchange the 
mean inter nucleon distance is considerably smaller and the average momentum (q) is 
supposed to be considerably larger. 

Since Ovfift decay is a two step process, in principle, one needs to construct and 
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sum over all the intermediate nuclear states, a formidable job indeed in the case of 
the shell model calculations. Since, however, the average neutrino momentum is much 
larger compared to the nuclear excitations, one can invoke closure using some average 
excitation energy (this does not apply in the case of 2v(3/3 decays). Thus one need 
construct only the initial and final + nuclear states. This is not useful in Quasiparticle 
Random Phase Approximation (QRPA), since one must construct the intermediate 
states anyway. In any case, it was explicitly shown, taking advantage of the momentum 
space formalism developed by Vergados |103j . that this approximation is very good 
[104[ 1105] . The same conclusion was reached independently by others [106J via a more 
complicated technique relying on coordinate space. 

Granted that one takes into account all the above ingredients in order to obtain 
quantitative answers for the lepton number violating parameters from the results of 
0^/9/3-decay experiments, it is necessary to evaluate the relevant nuclear matrix elements 
with high reliability. The most extensively used methods are the large basis Interacting 
Shell Model (ISM) calculations, (for a recent review see [65]) and QRPA( for a recent 
review see [66, 65J). The ISM is forced to use few single particle orbitals, while this 
restriction does not apply in the case of QRPA. The latter suffers, of course, from the 
approximations inherent in the RPA method. So a direct comparison between them is 
not possible. 

The shell model calculations have a long history pH [75J [Ml ttDSJ EDSl HUB QH] in 
double beta decay calculations. In recent years it has lead to large matrix calculations 
in traditional as well as Monte Carlo types of formulations [721 [73j [74j II 12|. II 13|. 1114] . 
For a more complete set of references as well as a discussion of the appropriate effective 
interactions see Ref. [65]). 

There have been a number of QRPA calculations covering almost all nuclear targets 
[TT51 ITT61 ITT71 ITT81 [1191 H201 IT2T1 IT221 IT231 [TM [T25l IT26] . These involve a number of 
collaborations, but the most extensive and complete calculations in one way or another 
include the Tuebingen group. We also have seen some refinements of QRPA, like proton 
neutron pairing and inclusion of renormalization effects due to Pauli principle corrections 
[1271 1128] . Other less conventional approaches, like operator expansion techniques have 
also been employed [66] . 

Recently, calculations based on the Projected Hartree-Fock-Bogoliubov (PHFB) 
method [129] , the Interacting Boson Model (IBM) [130] and the Energy Density 
Functional (EDF) method |131] entered the field of such calculations. The above 
schemes, in conjunction with the other improvements mentioned above, offer some 
optimism in our efforts for obtaining nuclear matrix elements accurate enough to allow 
us to extract reliable values of the lepton violating parameters from the data. 

As we have mentioned neutrinoless double beta decays (DBD) are concerned with 
fundamental properties of neutrinos. These properties arise out of interactions involving 
high energy scales, which are of great interests from view points of particle physics 
and cosmology. On the other hand, DBD processes are nuclear rare-decays in the 
low energy scale, which are studied experimentally by low-energy and low-background 
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nuclear spectroscopy, as given in review articles [301 EE] • 

Double beta decays are low-energy second-order weak processes with Qpp ~2 - 3 
MeV. Decay rates of 2i//3/3-decay within the SM are of the order of 10~ 20 /y, and the 
rates of 0z//9/9-decay beyond the SM are even many orders of magnitudes smaller than 
2i//3/3-decay rates, depending on the Qpp value and the effective Majorana neutrino mass 
(m„) (for definition see Eq. (I33jl ) in case of the light neutrino mass process. Then the 
O^/3/3-decay half lives are of the orders of Ty 2 « 10 27 y and 10 29 y in cases of the IH 
(inverted hierarchy) mass of (m„) ~ 30 meV and the NH (normal hierarchy) mass of 
(m u ) ~ 3 meV, respectively. 

For experimental studies of such rare decays, large detectors with ton-scale DBD 
isotopes are needed to get 0^(3 (3- decay signals in case of the IH v mass. Here the signals 
are very rare and are as low as Epp ^2 - 3 MeV. Background (BG) signal rates, however, 
are huge in the energy region of Eb < 3 MeV. Thus it is crucial to build ultra low BG 
detectors to find the rare and small Oz//3/3-decay signals among huge BGs in the low 
energy region. We are going to review this (see section [8]) in the case of most of the 
nuclear targets of experimental interest [301 EB H32] ( 76 Ge, 82 Se, 96 Zr, 100 Mo, 116 Cd, 
128 Te, 130 Te, 136 Xe, 150 Nd). 



3. The neutrino mass matrix in various models 



Within the SM of elementary particles, with the particle content of the gauge bosons 
A^, and W^'°, the Higgs scalar isodoublet $ = (0°, cj>~) (and its conjugate $*) and 
the fermion fields arranged in: 

• Isodoublets: (u a L, d a i) and (v a L, ^ai) for quarks and leptons respectively and 

• Isosinglets: u a R, d a R and e a R 

where a is a family index taking three values, the neutrinos are massless. They can not 
obtain mass after the symmetry breaking, like the quarks and the charged leptons, since 
the right handed neutrino is absent. 



3.1. Neutrino masses at tree level 

The minimal extension of the SM that would yield mass for the neutrino is to introduce 
an isosinglet right handed neutrino. Then one can have a Dirac mass term arising via 
coupling of the leptons and Higgs as follows: 

h a A^ L e aL ) f 0- ^ v m -»■ 

h Q A^a L e aL ) ( V ^ J vp R or = K,frj= (13) 

Thus one can have: 

^=(^,^( (m v f )(:i) (") 
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In the above expression, as well as in analogous expressions below, we have explicitly 
indicated not only the 6x6 matrix, but the states on which this matrix acts and m D is 
the above 3x3 matrix. One then obtains 6 Majorana eigenvectors which pair wise can be 
associated to opposite eigenvalues. Their sum and their difference may equally well be 
selected as physical states which correspond to three Dirac neutrinos and their charged 
conjugate (antineutrinos). This is fine within the above minimal extension. In grand 
unified theories, however, one is faced with the problem that these neutrinos are going to 
be very heavy with a mass similar to that of up quarks, which is clearly unacceptable. So 
such a model is inadequate^. Besides, such neutrinos viewed as Majorana with opposite 
CP eigenvalues or as Dirac particles cannot contribute to Quj3j3 decay. 

The next extension is to introduce a Majorana type mass involving the isosinglet 
neutrinos and an additional isosinglet Higgs field, which can acquire a large vacuum 
expectation value, an idea essentially put forward by Weinberg |134] long time ago. 
Thus the neutrino mass matrix becomes: 



Thus, provided that the Majorana mass matrix has only very large eigenvalues, one 
obtains an effective Majorana 3x3 matrix: 



which can provide small neutrino masses provided that the eigenvalues of the matrix Mr 
are sufficiently large. Mr can be arbitrarily large, since the the new scale, associated 
with the vacuum expectation of the isosinglet, does not affect the low energy scale 
arising from the vacuum expectation value of the standard Higgs particles. This is the 
celebrated see-saw mechanism. More precisely the type I see-saw mechanism, since, as 
we will see below, there exist other see-saw types (for a summary see, e.g., Abada et 



Thus with the above mechanism the neutrino flavors get admixed, the resulting 
eigenstates are Majorana particles and lepton number violating interactions, like 0vf3f3 
decay, become possible. 

Other extensions of the SM are possible, which do not require the addition of the 
right handed neutrinos but additional exotic scalars or fermions [18| Wf\. e.g. 

• An isotriplet A of Higgs scalars whose charge decomposition is 5 , 8~, 5°. 
Then this leads to the coupling: 



§ There may exist light Dirac neutrinos in theories formulated in extra dimensions, see e.g. the recent 
review by Smirnov [133] . If these neutrinos do not couple to the usual leptons they are of little interest 
to us. If they do and it so happens that the standard neutrinos are Majorana, they also become 
Majorana, except in the case of very fine tuning. 




(15) 




(16) 



q/[T35]l 
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Figure 2: The tree level contribution to the neutrino mass mediated by an isotriplet 
scalar. 



This after the isotriplet acquires a vacuum expectation value becomes 
yielding the neutrino Majorana mass matrix 

m a,p = ( h T)a,P^Va L Vp R (17) 

If, for some reason, the introduction of such an isotriplet, acquiring a small vacuum 
expectation value, is not preferred, the Majorana mass matrix can be obtained 
assuming that the isotriplet A possesses a cubic coupling /za with two standard 
Higgs doublets [4"4" j 11361 [1371 1138j (see Fig. [2]). Then one finds an effective Majorana 
neutrino mass from Eq. (fT7|) via the substitution 

v A v 2 /i A 



y/2 2 ml' 

where v/y/2 is the vacuum expectation value of the standard Higgs doublet (see 
Eq. ( TIB"]) ), niA is the mass of 5°. This is sometimes refer to as see-saw mechanism 

ii nam]. 

As we will see later this mechanism may lead to a new contribution to neutrinoless 
double beta decay via the direct decay of 5 into two electrons. 

An isotriplet of fermions with hypercharge zero (S + , S°, 

In this case the leptons couple to the isotriplet via the higgs doublet ( see [139| 1140] 
and references therein): 




-TPjJy/2 J 
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Figure 3: The tree level contribution to the neutrino mass mediated by a fermion 
isotriplet 



Then a subsequent coupling of the isotriplet to the leptons yields an effective 
Majorana coupling of the form: 

= (18) 

where is the mass of the neutral component of the isotriplet (see Fig. [3] and 
ref. [139J). It is sometimes referred as mechanism see-saw III. This mechanism, 
however, by itself cannot constitute a viable neutrino mass generator since it leads 
to two eigenstates with zero mass. This can be circumvented [1391 114U] . but then 
the model becomes more complicated. 



3.2. Neutrino masses at the loop level 

There many ways to obtain neutrino masses at the 1-loop level |141] , which have nicely 
been summarized by Smirnov [133] . We will only discuss one such case, which arises in 
the presence of R-parity violating supersymmetry, which leads to a viable neutrino mass 
spectrum |142j . in the sense that it can yield three massive neutrinos, if one includes 
not only the tree level contribution arising from the bilinear terms p. 43 J , but both type 
the 1-loop contributions |144[ 1145] shown in Fig. |H This is interesting, since, in such 
models, as we will see below, one can have particles other than neutrinos contributing to 
0i//9/3-decay. The above Majorana matrices are symmetric, in general complex, matrices. 



3.3. SUSY, GUT's and Family symmetries 

In many models, like the standard see-saw, the smallness of neutrino mass requires 
the existence of a heavy mass scale. The coexistence of two mass scales can naturally 
be accommodated in supersymmetry. In minimal supersymmetric extensions of the SM 
|146] one construct the see-saw matrix of Eq. ( fl5|) . see e.g. the review [IB]- Furthermore 
this can be extended to larger symmetries, e.g. two commuting symmetries, a Grand 
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Figure 4: The squark-quark (lower panel) and slepton-lepton (upper panel) 1-loop 
diagrams contributing to the neutrino mass in the presence of R-parity violating 
supersymmetric theories. 



Unified Symmetry G GUT and a family symmetry Gf. The family symmetry could be 
continuous, like SU (3), SU (2) or U (1) or discreet Z^, S3 or S3 x S3 etc. We are not going 
to further elaborate on such situations, which have been summarized in recent reviews 
[4U| I147[ 1148] to which we direct the interested reader. We should also mention that 
U(l) flavor symmetries arise naturally in superstring inspired models |149[ 1150] 1151] , in 
particular for the heterotic string and the 4-d fermionic constructions. D-brane models 
have also paved the way for completely new structures |152[ 11531 11541 [155J and in 
particular a very interesting formulation D-brane inspired mass textures |156j . 

All the models considered in this section lead to a light effective neutrino mass acting 
jointly or separately. Almost all of them involve parameters, which can be adjusted to 
fit phenomenology. It is not clear which one, if any, is going to ultimately be the 
theoretically preferred one. 

3.4- Neutrino mixing 

We have seen above that in general the neutrino mass matrix is a complex symmetric 
matrix. It can, however, be diagonalized by separate left and write unitary 
transformations, which can take the form [44] : 




where we have added the superscript to stress that they are the states entering the 
weak interactions. S^' are 3x3 matrices with S^ 11 ^ and approximately unitary, 
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while S"( 12 ) and are very small. (u' L , N' L ) and {y' R , N' R ) are the eigenvectors from 
the left and the right respectively. Thus the neutrino mass in the new basis takes the 
form: 

M u = jz ( m j 9 'jL u 'jR + MjN' jL N' jR ) + H.C. (20) 
j'=i 

This matrix can be brought into standard form by writing: 



and defining: 



v 



m . = \ mi \e~ ia \ Mj = \Mi\e 



Then 



M V = Y, ( m, + [M^NjNj) (21) 



Note, however, that: 

u c = v' jR + e^v' jL = e ia > (u' jL + e- ia >v' jR ) = e ia v 



j 



N c = N' jR + e^N' jr = (u' jR + e~^N' jR ) = N 3 (22) 
i.e. they are Majorana neutrinos with the given Majorana phases. Furthermore 



iR 



The second of Eqs. (119]) can now be written as 

where e ta and e** are diagonal matrices containing the above Majorana phases. The 
neutrinos interact with the charged leptons via the charged current (see below). So 
the effective coupling of the neutrinos to the charged leptons involves the mixing of the 
electrons S e . Thus the standard mixing matrix appearing in the absence of right-handed 
neutrinos is: 

U PMNS = U = = (S^) + (24) 
The other entries are defined analogously: 

Um = (s^Y S {i] \ (ij) = (12), (21), (22) (25) 
In particular the usual electronic neutrino is written as: 

^e^EK'VC^) ( 26 ) 



v, 



:» = EK%<^) (27) 
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In other words the left handed neutrino may have a small heavy component, while the 
right handed neutrino may have a small light component. Note also that the neutrinos 
appearing in weak interactions can be Majorana particles in the special case that all 
Majorana phases are the same. 

The Pontecorvo-Maki-Nakagawa-Sakata neutrino mixing matrix Upmns is 
parametrized by 



UpMNS — R23R13R12, 

where the matrices Rij are rotations in ij space, i.e., 



(28) 



R 
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( 1 \ 

C12 S12 
-S12 C12 y 

Cl2 S12 \ 

-S12 C12 
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s 13 e~ 
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where 



Cjj = cos 



Sjj = sin 1 



(29) 



(30) 



#12, #13 and #23 and three mixing angles and S is the CP- violating phase. Then, we get 



PMNS 



C13S12 



e~ iS s 13 



C12C13 

-C23S12 - e'"Ci 2 Sl3S 2 3 C12C23 - < '"•\!2-\l:;- s '2:l 'M:v s 2:i 

,5„ „ „ „i5 



AS, 



AS, 



,(31) 



\ S12S23 — e C12C23S13 — e""C23Si2Si3 — C12S23 C13C23 y 

If neutrinos are Majorana particles Upmns i n Eq. (!3T|) is multiplied by a diagonal phase 
matrix P, which contains two additional CP-violating Majorana phases a\ and a.2- 



P = diag(e iai ,e ia2 ,e lS ). 



(32) 



4. The elusive absolute scale of the neutrino mass 

With the discovery of neutrino oscillations quite a lot of information regarding the 
neutrino sector has become available (e.g., for recent reviews see [HHJ 1157] ). More 
specifically we know: 

• The mixing angles 812 and 623 and we have both a lower and an upper bound on 
the small angle #13 

• we know the mass squared differences: 

a Ln = A L = m\- ml, and A^ TM = |A^ 3 | = \m 2 3 - m\\ 

entering the solar and atmospheric neutrino oscillation experiments. Note that we 
do not know the absolute scale of the neutrino mass and the sign of A23 
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For determination of an absolute scale of the neutrino mass the relevant neutrino 
oscillation parameters are the MINOS value Am 2 TM = (2.43 ± 0.13) x 1(T 3 eV 2 
[158] . the global fit value Am 2 SVN = (7.65ig$j) x 10" 5 eV 2 (STJ, the solar-KamLAND 
value tan 2 # 12 = 0.452^;^ |159] and the recent Daya Bay observation sin 2 2# 13 = 
0.092 ± 0.016 (stat) ± 0.005 (syst) with a significance of 5.2 standard deviations 
[160J. We note that non-zero value of mixing angle #13 was already observed also 
by the T2K (0.04 < sin 2 2fl 13 < 0.34) [1ST], the DOUBLE CHOOZ (sin 2 2fl 13 = 
0.085 ± 0.029 (stat) ± 0.042 (syst) (68% CL)) [162] and the RENO (sin 2 2#i 3 = 
0.103 ±0.013 (stat) ±0.011 (syst)) [163] collaborations. 
Based on the above we have the following scenarios: 

• Normal Spectrum (NS), ni\ < m 2 < m 3 : 

A 2 2 2 A 2 22 

m suN = m 2 -m 1 , Am ATM = m 3 -m 1 



m = mi, m 2 = \J Am| UN ± m 2 , , m 3 = yj Am 2 ATM + m. 
Inverted Spectrum (IS), m 3 < mi < m 2 : 

A2 22a2 22 

Am SUN = m 2 -m l , Am ATM = m 2 - m 3 
m = m 3 , m 2 = Am 2 ATM + m 2 ,, 



m x = yj Am2 TM - Am2 UN + ml 

The absolute scale m of neutrino mass can in principle be determined by the following 
observations: 

• Neutrinoless double beta decay. 

As we shall see later (section [5]) the effective light neutrino mass (m„) extracted in 
such experiments is given as follows [HJ 1164] : 

k> = e (uiiy m k 

k 

2 2 2iai 2 2 2iao 2 / O o\ 

= Ci 2 q 3 e x m\ + c n s 12 e 2 m 2 + s 13 m^. (33) 

• The neutrino mass extracted from ordinary beta decay, e.g. from triton decay 
[T65 lfT66lfT67] . 



{m v ) 



decay 



\ k 



en) i 

ek 



-mr 



'cl 2 cl 3 ml + c\ 3 s\ 2 m\ + s 2 13 m 2 3 . (34) 
assuming, of course, that the three neutrino states cannot be resolved. 
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From astrophysical and cosmological observations (see, e.g., the recent 
summary |168j ) . 



m v =^2m k < m astro (35) 

k 

The current limit m astro depends on the type of observation |168] . Thus CMB 
primordial gives 1.3 eV, CMB+distance 0.58 eV, galaxy distribution and and lensing 
of galaxies 0.6 eV. On the other hand the largest photometric red shift survey 
yields 0.28 eV |169j . For purposes of illustration we will take a world average of 
Castro = 0.71 eV. 

The above mass combinations can be written as follows: 

(i) Normal Hierarchy (NH), mi <C m 2 <C m 3 . 
In this case we have 

A 2 2 2 a 2 2 2 

m suN = m 2 ~ m i , atm =m 3 -m 1 



Thus: 

• Triton decay. 



'c\ 2 c\^n\ + sf 2 c? 3 (Am| UN + m 2 ) + s\ z (Am\ Tm + mf) (36) 
Cosmological bound: 



{m v ) 



decay 



m v = m 1 + y] Am 2 SVN + m\ + sj (Am^ + mj) (37) 

Ov(3(3 decay: 

(m v ) = c\ 2 c\ z m x e 2 ^ + s 2 12 c 2 13 e 2ia * J Am 2 UN + mj 



+ S ? 3V /(Ami TM + mf), (38) 

where ot\ and a 2 (a 3 = 0) are Majorana CP violating phases of the elements 
U el and U e2 with U ej = \U ej \e ia ^ (i = 1,2). 

By assuming NH, i.e., that mi is negligibly small (mi <C ^Am^, m 2 — 
'Amg UN , and m 3 ~ ^Am^ TM ), we obtain 



< 4 • 10~ 3 eV. (39) 
(ii) Inverted Hierarchy (IH), m 3 <C rri\ < m 2 : In this case we have 

a2 2 2 a 2 22 

Am SUN = m 2 -m l , Am ATM = m 2 - m 3 

• Triton decay. 

Wdccay = (s 13 m 2 + s 2 2 c 2 3 (Ami TM + m 2 ) 

+c 2 2 c 2 3 (Ami TM - Am 2 UN + m 2 )) V2 (40) 
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Figure 5: (Color online) The neutrino mass limits in eV as a function of mass of the 
lowest eigenstate m also in eV, extracted from cosmology (left panel) triton decay (right 
panel). From the current upper limit of 2.2 eV of the Mainz and Troitsk experiments we 
deduce a lowest neutrino mass of 2.2 eV both for the NS and IS. From the astrophysical 
limit value of 0.71 eV the corresponding neutrino mass extracted is about 0.23 eV for 
the NS and IS. It is assumed: Ami TM = (2.43 ± 0.13) x 10~ 3 eV 2 [158] . Amf UN = 
(7-651^) x 10~ 5 eV 2 [El], tan 2 # 12 = 0.4521^1 [H3] and sin 2 2^ 3 = 0.092 ± 0.016 
[T60] . 



• Cosmological bound: 

m u = m 3 + ^j (Am 2 ATM + m 2 3 ) + yj (Am 2 TM - Am 2 UN + m§) (41) 

• 0vf3(3 decay: 

(m u ) = s 2 13 m 3 + s 2 12 cl 3 e 2ia ^(Aml TM + ml) 

+ cl 2 c 2 13 e 2 ^ ^(Aml TM -Aml UN + m 2 ). (42) 
Since in IH scenario m 3 is negligibly small (mi ~ m 2 — J Am 2 tm and 
m 3 < yjAm 2 ATM ), we find 

\{m v )\ ~ ^Am2 tm c 2 3 (l-sin 2 2^ 12 sin 2 « 12 )^, (43) 

where ayi = a 2 — a±. The phase difference cti 2 is the only unknown parameter 
in the expression for |(m„)|. From (T4"3"j) we obtain the following inequality |170] 

1.5 • 10~ 2 eV < Km*) | < 5.0 • 10~ 2 eV. (44) 

(iii) Quasi-degenerate (QD) spectrum, mo = mi ~ m 2 — m 3 . Then 

• Triton decay and Cosmology: 

(m v ) Aecay = m , m u = 3m (45) 

• 0vf3(3 decay: 

The effective Majorana mass is relatively large in this case and for both types 
of the neutrino mass spectrum is given by the expression 

m |l - 2c 2 3 c 2 2 | < \{m u ) \ < m . (46) 
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The above results are exhibited in Fig. E]fbr the tritium /3-decay and cosmological 
limits as a function of the lowest neutrino mass and in Fig. [6] for the case of the 
decay both for the NS and the IS scenarios. The allowed range values of \(m v )\ as a 
function of the lowest mass eigenstate m is exhibited. For the values of the parameter 
sin 2 2# 13 new Double Chooz data are used |162j . The IH allowed region for \(m u )\ is 
presented by the region between two parallel lines in the upper part of Fig. [6j The 
NH allowed region for | (m u ) | ~ few meV is compatible with m.o smaller than 10 meV. 
The quasi-degenerate spectrum can be determined, if mo is known from future /3-decay 
experiments KATRIN |165[ 1166] and MARE |167] or from cosmological observations. 
The lowest value for the sum of the neutrino masses, which can be reached in future 
cosmological measurements |171l 11721 1173] . is about (0.05-0.1) eV. The corresponding 
values of m are in the region, where the IS and the NS predictions for \(m u )\ differ 
significantly from each other. 

From the most precise experiments on the search for O^/3/3-decay [36| |39| HQ] by 
using of nuclear matrix elements of Ref. |174] the following stringent bounds were 
inferred (see Table CP 

\(m u )\ < (0.20 - 0.32) eV ( 76 Ge), 

< (0.33 - 0.46) eV ( 130 Te), 

< (0.17-0.30) eV ( 136 Xe). (47) 

These bounds we obtained using the Ozz/S/S-decay NMEs of |174] calculated with 
Brueckner two-nucleon short-range correlations. There exist, however, a claim of the 
observation of the Oz//3/3-decay of 76 Ge made by some participants of the Heidelberg- 
Moscow collaboration [12]. Their estimated value of the effective Majorana mass 
(assuming a specific value for the NME) is | (m u ) | ~ 0.4 eV. This result will be checked 
by an independent experiment relatively soon. In the new germanium experiment 
GERDA (131 1175] , the Heidelberg-Moscow sensitivity will be reached in about one year 
of measuring time. 

In future experiments, CUORE[I76], EXOpTTO EE], MAJORANA [179], Su- 
perNEMO [HQ], SNO+ [M], KamLAND-Zen and others [6j EE CEE2] , a sensitivity 

\(m v )\ ~ a few 10~ 2 eV (48) 

is planned to be reached, what is the region of the IH of neutrino masses. In the case of 
the normal mass hierarchy | (m u ) | is too small in order to be probed in the O^/3/3-decay 
experiments of the next generation. 

Recently, however, for the explanation of the Reactor Antneutrino Anomaly |185j . 
a light sterile neutrino has been introduced with mass squared difference: 

Am 2 2A = \m\ - ml\ w Am 2 4 = \m\ - m\\ > 1.5(eV) 2 . (49) 

which couples to the electron neutrino with a mixing angle: 

sin 2 29 u = 0.14 ± 0.08 (95% C.L.). (50) 
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m [eV] 

Figure 6: (Color online) We show the allowed range of values for | (ttt,^) | as a function 
of the lowest mass eigenstate m using the three standard neutrinos for the cases of 
normal (NS, mo = mi) and inverted (IS, mo = m.3) spectrum of neutrino masses. Also 
shown are the current experimental limits and the expected future results [183j (QRPA 
NMEs with CD-Bonn short-range correlations and g e J^ = 1.25 are assumed [1741 1184] ). 
Note that in the inverted hierarchy there is a lower bound, which means that in such 
a scenario the 0^/3/3-decay should definitely be observed, if the experiments reach the 
required level. The same set of neutrino oscillation parameters as in Fig. |5]is considered. 



On the other hand in a recent global analysis more than one sterile neutrino are needed 
[186] , with somewhat smaller mass squared differences, but similar couplings. Due to 
such a mixing, even if their couplings are of the usual Dirac type, the resulting mass 
eigenstates are of the Majorana type due to their coupling to the usual neutrino. 

The U(4 x 4) neutrino mixing matrix in the presence of one sterile neutrino with a 
small mixing becomes |187j : 

U = -R34-R24-Rl4-R23-Rl3-Rl2-P- (51) 

It depends on 6 mixing angles (6*14, 6*24, #34, 6*12, 6*13, 6*23, 3 Dirac (#24, Su, 5 13) and 3 
Majorana (a%, «2, 0:3) CP- violating phases entering the diagonal P matrix : 

P = diag U ia \ e ia \ e l{a3+5l3 \ e i5li ) . (52) 

Similarly, one can parametrized the (7(5 x 5) mixing matrix for two sterile neutrinos as 
(10 mixing angles and 5+4 CP violating phases) 

U = -R45-R35-R34-R25-R24-R23-R15-R14-R13-R12-P, (53) 
where P = diag (e ia \ e ia \ e l{a3+Sl3 \ e^+^i^ e iSi 5 \ 
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Figure 7: (Color online) The same as in Fig. [61 if one considers one (version 3+1) or 
two (version 3+2) sterile neutrinos, which are heavier than the standard neutrinos. Best 
fit points for the 3+1 (Am^ = 1.78 eV 2 , U e4 = 0.151) and 3+2 (Am^ = 0.46 eV 2 , 
= 0.108 and Am^ = 0.89 eV 2 , U e 5 = 0.124) scenarios from reactor antineutrino 
data are taken into account |186j . In addition, best fit values Am 2 TM = 2.43 x 10~ 3 eV 2 
[158] . Am 2 UN = 7.65 x 10~ 5 eV 2 [81], tan 2 # 12 = 0.452 [159] and sin 2 2fl 13 = 0.092 + 0.016 
|160] are assumed. 

If the presence of one sterile neutrino, the effective neutrino mass in the Oz^/3/3-decay 
is given by [187J 

I / \ I I 2 2 2 2iai ,2 2 2 2iao 

3+i\ = \c 12 c 13 c u e 1 m 1 + c 13 c u s 12 e 2 m 2 

+ c 2 4 s 2 3 e 2iQ3 m 3 + s 2 4 m 4 |. (54) 

We note that from three additional angles (rn v ) \ depends only one of them, namely #14. 
If there are two sterile neutrinos we end up with 

I / \ 1 12 2 2 2 2iqi 2 2 2 2 2i«9 

K^)3+2| = |c 12 c 13 c 14 c 15 e 1 m l + c 13 c 14 c 15 s 12 e 2 m 2 

+ c 2 u c 2 15 s 2 13 e 2ia *m 3 + c 2 5 s 2 4 e 2iQ4 m 4 + sj 5 m 5 \. (55) 

Due to the extra terms in Eqs. (|54|) and (|55|) and their couplings, depending on the extra 
majorana phases, the sterile could dominate, increase or deplete (m u ). By assuming best 
fit values for \U e4 \, m 41 and \U e s\, m 2 x from reactor antineutrino data |186j the allowed 
range values of (m u ) as a function as a function of the lowest mass eigenstate mo in the 
presence of one and two sterile neutrinos is shown in Fig. [71 
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5. The intermediate Majorana neutrino mechanism in Ou/3/3 decay 



To proceed further in our study of the neutrino mediated Oz//3/3-decay process it is 
necessary to study the structure the effective weak beta decay Hamiltonian. In general it 
has both left handed and right handed components. Within the SU(2) L x SU (2) R x U (1) 
at low energies it takes the current- current form: 



\/2 



J 



eJ 



kJ> 



R 



h.c. 



Here, e is mixing of Wl and Wr gauge bosons 

Wl = cos eWi - sin eW 2 , Wr = sin eWi + cos eW 2 



(56) 



(57) 



where W\ and W 2 are the mass eigenstates of the gauge bosons with masses Jl%j 
and J% 2 , respectively. The mixing is assumed to be small, sine ~ e, cose ~ 1, and 
mw 1 ~ ^w L i fnw 2 ~ m w H - « is the mass squared ratio 



rn 



K 



Hi 



rn 



(5f 



H' 2 



is the standard hadronic current in V-A theory: 



T/4 



£ 



SV7 + ipAf 



a 



2m. 



u, 



(59) 



where u p (i) and u n (i) are the spinors describing the proton and neutron i. m p is the 
nucleon mass and g M is the momentum transfer, gy = gy(g 2 ), #m = 9Ai(q 2 ), 9a = 9a{q 2 ) 
and (7p = gp(q 2 ) are respectively the vector, weak-magnetism, axial- vector and induced 
pseudoscalar form-factors. 

We will see below that it is necessary to also consider the right handed current Jr 
of the form V+A, 



J: 



R 



(7 



fly^ + l 9M^~qu + 9 ATI* + gpTlB 



u. 



(0, 



(60) 



even though normally one expects its contribution to ordinary beta decay to be 
suppressed by a factor k or e. Some relations among form factors in J^' and Jr 
are considered because the strong and electromagnetic interactions conserve parity. 

^l(^r) and rj^L^ep) are ^- e ^ operators representing the left (right) handed electrons 
and electron neutrinos in a weak interaction basis in which the charged leptons are 
diagonal. We have seen above the the weak neutrino eigenstates can be expressed in 
terms of the propagating mass eigenstates [44] (see Eqs. ( l26l) and ( |27l) ). The mass 
eigenstates Vk,N}. satisfy the Majorana condition: f^Cfe = C v\t AfcHfc = C N k , where 
C denotes the charge conjugation and £, H are phase factors, which guarantee that the 
eigenmasses are positive = e 



see Eq. fl22 
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Figure 8: The neutrino mass contribution at the nuclear level in the presence of left 
handed currents for light intermediate neutrino (a) and heavy neutrino (b). 



Before proceeding further we should mention that in the context of the above 
interaction neutrinoless double beta decay is a two step process. The neutrino is 
produced via the lepton current in one vertex and propagates in the other vertex. If the 
two current helicities are the same one picks out of the neutrino propagator the term: 

_P_ [ m^, m) «? (61) 
q z — rrij { —l/rrij, wf- S> q 

where q is the momentum transferred by the neutrino. In other words the amplitude 
for light neutrino becomes proportional to its mass, but for a heavy neutrino inversely 
proportional to its mass. 

If the leptonic currents have opposite chirality the one picks out of the neutrino 
propagator the term: 

a ~ 2 -> 4' m ) < I 2 ( 62 ) 
q z — rrij q z 

i.e. in the interesting case of light neutrino the amplitude does not explicitly depend on 
the neutrino mass. The kinematics becomes different than that for the mass term. 



5.1. The Majorana neutrino mass mechanism 

This mechanism is the most popular and most commonly discussed in the literature (see 
Fig. From this figure we can see read off the couplings and the phases. We have also 
seen that for currents of the same chirality one picks out the mass of the propagating 
neutrino. Thus the lepton violating parameter is defined as (m u )/m e with (m u ) given 
byEq. AMD- 

We will consider only 0+ — > 0^ transitions. Then both outgoing electrons are in 
the S1/2 state. Thus for the ground state transition, restricting ourselves to the usual 
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Figure 9: The neutrino mass contribution at the nuclear level in the presence of right 
handed currents proceeding via the right handed boson (a) (k 2 suppression factor), the 
combination of the right handed boson and the mixing between the two bosons (b) and 
(c) (kxc suppression factor) and via the light gauge boson via its mixing with its heavy 
partner (d) (e 2 suppression factor). 



light left handed neutrino mass mechanism, we obtain the following expression for the 
O^/3/9-decay inverse half-life: 

2 



Ois 1 
1/2J 



G, 



01 



m: 



(V 



(63) 



Extraction of {m v ) from the above life time will have a wide range of implications in 
physics, it can, e.g., constrain the baryon asymmetry in the universe [188] \Yb\ etc. 

Another less popular possibility is the mass contribution arising from DBD in the 
presence of right handed currents (see Fig. [9]) or heavy neutrinos in general fl 89 j . The 
above expression can be generalized to include many mass mechanisms [HI [62], [63] EU 
as follows: 



[T° u 



/2J 



Goi[|^l| 2 + \Xr\ 2 — C X X L X R + .. 



(64) 
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The coefficient C[ is negligible, because these terms do not interfere to leading order due 
to the different helicity of the 2 electrons. Here ... indicate other non traditional modes 
(SUSY etc.). The nuclear matrix elements entering the above expression are given in 
units of M GT and are denoted [63] by x- 

X L = ^M U » + V L N M° N », X R = r l *M$', (65) 

where the nuclear matrix elements M® u and Mjtf will be discussed later (see section 
[9]). The subscripts L(R) indicate left (right) handed currents respectively. The lepton- 
number non- conserving parameters, i.e. the effective neutrino mass (m u ) given by Eq. 
fl33|) and rj^ ,rj^ are given as follows |44j : 

4 = E {UT"? —A ^, (66) 

k Mk 

V R N = (K 2 + e> + 2e K )j: {Utlf (67) 

k M k 

with m p (m e ) being the proton (electron) mass. Gqi is the integrated kinematical factor 
[53j EU [63], [HT] . The nuclear matrix elements associated with the exchange of light and 
heavy neutrino must be computed in a suitable nuclear model. The ellipses ... mean 
that Eq. fl64"j) can be generalized to the mass term resulting from any other intermediate 
fermion. 

At this point we should stress that the main suppression in the case of light neutrinos 
comes from the smallness of neutrino masses. In the case of heavy neutrino not only 
from the large values of neutrino masses but the small couplings, for the left handed 
neutrinos and k and e for the right-handed ones. 

5.2. The neutrino mass independent mechanism (leptonic left-right interference, A and 
i] terms). 

As we have already mentioned in the presence of right handed currents one can have 
interference between the leptonic currents of opposite chirality (see Fig. [TOl) . The 
elementary amplitude is now proportional to the 4-momentum transfer. We thus have 
a space component and a time component in the relevant amplitude. This leads to 
different kinematical functions and to two new lepton number violating parameters [H] 
(A) and (77) defined by 

(V) = eVRL, W = KVRL, VRL = j: (U^U^)^. (68) 

1 

The parameters (A) and (77) are small not only due to the smallness of the parameters 
k and e but, in addition, due to the smallness of U^ 2l \ As we have already mentioned 
the (A) can also have a different origin (see section l6.2.4p . 

All the above contributions, even though the relevant amplitudes are not explicitly 
dependent on the neutrino mass, vanish in the limit in which the neutrino is a Dirac 
particle. 
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Figure 10: The Feynman diagrams at the nucleon level when the leptonic currents are of 
opposite chirality leading to the dimensionless lepton number violating parameters (A) ( 
panels (a) and (b) of the figure) and (rj) (panels (c) and (d) of the figure) of Oi//3/3-decay. 
Note that in part (a) the process proceeds via the right handed vector boson, while in 
part (b) through the mixing of the left and right handed bosons. 



The above expression for the lifetime is now modified to yield \%5[ 1190j : 

[T1J2V = G^i\M° G v T \ 2 {\X L \ 2 + lA^I 2 - C[X L X R + ... 
+C 2 \ (A) \X L cos^ + C 3 \( V ) \X L cosi, 2 + C 4 \ (A) | 2 
+C 5 1 ( V ) | 2 + C 6 1 (A) 1 1 ( V ) [cos^ - ^2) + Re(C 2 (\)X R 
+C 3 (r ] )X R )}, 

(69) 

where X L and X R are defined in Eq. fl65l) . ip\ and ip 2 are the relative phases between X L 
and A and X^ and rj respectively. The coefficient C[, representing the mixing between 
the left and the right handed currents is kinematically suppressed [191j . The ellipses {...} 
indicate contributions arising from other particles, e.g., intermediate SUSY particles or 
unusual particles which are predicted by superstring models or exotic Higgs scalars etc 
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(see below section E]). 

Many nuclear matrix elements appear in this case, but they are fairly well known 
and they are not going to be reviewed here in detail (see e.g. jHJ [62] and [631 EH EHJ [66] ) . 
For the reader's convenience we are only going to briefly indicate them in our notation 
[87J) the additional nuclear matrix elements, not encountered in the mass mechanism. 
These are: XFu>, Xgtuj, Xr, Xi±, X2± Xf, Xgt, Xt, Xp where 

XFu = ( — ) "TT ? (70) 
g A M GT 

Xgtu, = -rj — , (71) 



M GT 

R 

M GT 



Mr . . 

XR = TT (72) 



and 



Xi± = ±3 X ' F + x'gt - QXt, (73) 

X2± = ±XFuj + XGTu; ~ gXl± (74) 

where [87] (x'f = M' F /M GT etc for the space part) and (xfui = M Ful /M GT etc for the time 
component). In the limit in which the average energy denominator can be neglected 
we obtain 



Xf = Xf = Xfu,, (75) 
Xgt = x'gt = Xgtcj = 1- (76) 

The quantities Gqi have been tabulated [53l EH EHl [871 H92j . see also [65] for a 
recent review. The coefficients C[, Ci, i = 2—6 are combinations of kinematical functions 
and the nuclear matrix elements have been previously discussed [192j . The coefficients 
Ci,i = 2, ...,6 with and without p-n pairing can be found in the literature [87j. For a 
more conventional formulation, restricted, however, in the light neutrino sector, see [65J. 

It is worth mentioning that in the case of the rj, in addition to the usual Fermi 
Gamow- Teller and tensor terms, we have additional contributions coming from the 
nucleon recoil term (xr) and the kinematically favored spin antisymmetric term (xp)- 
Due to these two effects the limit extracted for 77 is much smaller than that for A [87J. 
Effective operators of a similar structure also appear in the context of R-parity violating 
interactions when a neutrino appears in the intermediate states (see below). 

There seem to be significant changes in the nuclear matrix elements, when the 
p — n pairing is incorporated. This point needs special care and further exploration is 
necessary. It has only been examined in some exactly soluble models, e.g. SO (8), or 
better approximation schemes [1931 1194] . but only in connection with the 2i//3/3-decay, 
or shell model calculations but for systems, which do not double beta decay [195J. 

Returning back to the question of the availability of nuclear matrix elements 
relevant for neutrinoless double beta decay, we refer once again to existing excellent 
recent reviews [651 ESI I192j . These reviews also provide a more detailed description of 
the nuclear models employed. 
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Figure 11: The Feynman diagrams at quark level leading to majoron emission in the 
0i//9/3-decay instead of the more well known mass term. Here x° stands for the majoron, 
not to be confused with the neutralino, which we will encounter later in connection with 
supersymmetry. 



5.3. Another neutrino mass independent mechanism (majoron emission) 

It is well known that in some theories lepton number is associated with a global, not a 
local, symmetry. When such theories are broken spontaneously, one encounters physical 
Nambu-Goldstone bosons, called majorons. These bosons only couple to the neutrinos. 
So in any model which gives rise to mass term for the light neutrino (mass insertion 
in the neutrino propagator), one may naturally have a competing majoron-neutrino- 
antineutrino coupling. Such a mechanism is shown at the quark level in Fig. [TTJ 
The majoron, which couples to the left handed neutrinos, comes from the neutral 
member of the isotriplet. Such a multiplet, however, cannot easily be be accommodated 
theoretically. So this type of majoron is not present in the usual models. On the other 
hand there is a majoron x°, the imaginary part of an isosinglet scalar, which couples 
to the right handed neutrino with a coupling g\-. This gives rise to the mechanism 
shown in Fig. [12] at the nucleon level. The right handed neutrino, however, has a small 
component of light neutrinos (see Eq. (127]) ). 

£vu X ° = ZJ 9ij [ViLlsVjL] X°, (77) 

i<j 

with 

9^ = T.U^ufpg%^ (78) 

i<j 

with g a p the coupling of the majoron to the corresponding neutrino flavors. The 
expression for the half-life takes the form 

[Tifc]' 1 = G oi\(g)M° u \ 2 , (79) 

with (g) = J2i<j ~U ei U e j gij- Notice that, even if g^ takes natural values, the coupling 
gij is very small due to the smallness of the mixing matrix U^ 2l \ Thus the effective 
coupling(g) is very small. So, even though we do not suffer in this case from the 
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Figure 12: The same process as in Fig [IT], but written at the nucleon level in the case 
of the isosinglet majoron, which couples to the right handed neutrinos. 



suppression due to the smallness of the mass of the neutrino, the majoron emission 
mechanism is perhaps unobservable. There exist, however, exotic models, which, in 
principle, may allow majoron emission with a large coupling, like the bulk majoron [196J 
and others [1971 1198] , which we are not going to pursue this theoretically any further. In 
any case it is straightforward to extract the limits on the effective coupling (g) , since the 
nuclear matrix elements are the same as in the light neutrino mass mechanism. Note, 
however, that the spectrum of the summed energy of the two electrons is continuous 
and the kinematical function is different. 



6. Mechanisms in 0z//9/3-decay not involving intermediate neutrinos 

6.1. The direct decay of doubly charged particles to leptons 

Such a candidate is the isotriplet scalar, which can generate Majorana neutrino mass 
as see-saw mechanism II. The doubly charged component can directly decay into two 
leptons (see Fig. [TBI .The coupling to the quarks is achieved via the charged Higgs 
isodoublet in models where it survives the Higgs mechanism (e.g. SUSY) or through 
gauge bosons (Figs [T3a and [T3b respectively). For consistency between the two we 
slightly change the notation here implying that the cubic coupling /i^ originates from a 
dimensionless quartic coupling A of the Higgs particles of the figure with an isosinqlet 
scalar acquiring a vacuum expectation value Vr. 

The lepton number violating parameter associated with Fig. [13] is given by: 

2y/2 Xv R m 2 u m p g ee v R m 2 w m p g ee 

VAH ~ rn*rn* r ' VAW - 4 2 , r («U) 

Taking g ee = 1, a natural value, and Xvr ~ 1 TeV, = m& = 100 GeV and m u = 10 
MeV, one finds jH] tjah ~ 3 x 10~ 8 . Similarly taking vr = mw R = ^A fl = rn-w I = 
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(a) 



(b) 



Figure 13: We show the direct decay of the isotriplet into two leptons (a). The coupling 
to the quarks is achieved via the charged Higgs isodoublet in models where it survives 
the Higgs mechanism (e.g. SUSY). In (b) we show the direct decay into two leptons of an 
isosinglet doubly charged Higgs, present in left-right symmetric models. The coupling 
to the quarks is now achieved via the right handed gauge bosons. 

10m\y, we find t/avk ~ 10~ 6 . 

Note that the term 7]&h adds to Xl of Eq. ( 1641) . while rj^w adds to the Xr. 

6.2. The R-parity violating contribution to 0^(3 (3 decay. 
In SUSY theories R-parity is defined as 



with B =baryon, L =lepton numbers and s the spin. It is +1 for ordinary particles and 
-1 for their superpartners. R-parity violation has recently been seriously considered in 
SUSY models. It allows additional terms in the superpotential given by 



where a summation over the flavor indices i, j, k and the isospin indices a,b is understood 
( Xijk is antisymmetric in the indices i and j). The last term has no bearing in our 
discussion, but we will assume that it vanishes due to some discreet symmetry to avoid 
too fast proton decay. The first term is a lepton number violating bilinear and, since 
it cannot be rotated away, it can lead to neutrinoless double beta decay. The A's are 
dimensionless couplings not predicted by the theory. The couplings are assumed to be 
given in the basis in which the charged fermions are diagonal. In the above notation 
L,Q are isodoublet and E c , D c isosinglet chiral superfields, i.e they represent both the 
fermion and the scalar components. 

The above R-parity violating superpotential can lead to Majorana neutrino masses 
without the need of introducing the right-handed neutrino and invoking the see-saw 
mechanism [951 H42j . One then can have contributions to neutrinoless double beta 
decay in the usual way via intermediate massive neutrinos as discussed above. 





W = e i L a i H 2a + XijkLiL^E^eab + \' ijk LtUjD c k e ab + \" jk UiUjD c k , 



(82) 
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Figure 14: The R-parity violating contribution to 0v(3f3 decay mediated by neutralinos 
arising from the bilinear terms in the superpotential. For comparison we give the 
neutrino mediated process of Fig. [8] expressed at the quark level. 



6.2.1. The contribution arising from the bilinears in the superpotential. 
The first term in the superpotential, Eq. (1821) . can cause mixing between the neutrino 
and the neutralinos as soon as the s-neutrino develops a vacuum expectation value. As a 
result it can directly lead to neutrinoless double beta decay j9U [95j H99] via the effective 
W-charged lepton-neutralino interaction: 

c = ~72 Kn W ^ L ^ L - (83) 

where K n is a dimensionless quantity, associated with each of the four neutralinos, which 
arises due to neutrino neutralino mixing |94j. This term then gives rise to a diagram 
analogous to that of Fig. [8] with the intermediate particle now being the neutralino, 
which is heavy with mass M x o and leads to a short ranged operator. For the reader's 
convenience this is shown in Fig. [141 at the quark level. One thus obtains an analogous 
lepton number violating parameter: 

V L N rfr = £ Kl^. (84) 

6.2.2. The contribution arising from the cubic terms in the superpotential. It has also 
been recognized quite sometime ago that the second and third terms (cubic terms) in 
the superpotential could lead to neutrinoless double beta decay [26j [881 ES] . Typical 
diagrams at the quark level are shown in Fig. [15j Note that as intermediate states, 
in addition to the sleptons and squarks, one must consider the neutralinos, 4 states 
which are linear combinations of the gauginos and higgsinos, and the colored gluinos 
(supersymmetric partners of the gluons). Whenever the process is mediated by gluons 
a Fierz transformation is needed to lead to a colorless combination. The same thing is 
necessary whenever the fermion line connects a quark to a lepton. As a result one gets 
at the quark level not only scalar (S) and pseudoscalar (P) couplings, but tensor (T) 
couplings as well. This must be contrasted to the V and A structure of the traditional 
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Figure 15: The R-parity violating contribution to 0iyj3j3 decay mediated by sfermions 
and neutralinos (gluinos) arising from the cubic terms in the superpotential. The dots 
indicate the lepton violating R-parity interactions. 



mechanisms. One, therefore, must face the problem of how to transform these operators 
from the quark to the nucleon level. 



6. 2. 3. The lepton number violating parameters from the cubic terms without intermediate 

neutrinos. As we have mentioned the effective lepton number violating parameter 

arising from the bilinear terms in the superpotential is analogous to that arising from 

the heavy intermediate neutrinos and, thus, it will not be further discussed. We will 

concentrate on the cubic terms in the superpotential [261 E3 I200[ I201[ 1202] . Then the 

effective lepton number violating parameter arising from these terms, assuming that the 

pion exchange mode dominates, as the authors of Refs. [26, 66J find, can be written as 

3 

Vsusy = (A' m ) 2 -(XPS Vps + Vt), (85) 

with J]ps{t]t) associated with the scalar and pseudoscalar (tensor) quark couplings given 
by 

VPS = Vx,e + Vx,<j + VxJ + Vg + 7V§, (86) 
Vt = Vx,<i - Vxj + % - V'g- (87) 

These authors find xps — (5/3), but, as we shall see below, it depends, in general, on 
ratios of nuclear matrix elements. For the diagram of Fig. fT5h one finds 



(G F m 2 w Y 



K~) 2 A>e~e~. 
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For the diagram of Fig. [ToT o one finds 

= 2^kr liKif{ 'i hi + (K ' )2< 5 )aa1, (89) 

* = l a 'W^f^ 2 + (Kfi)2 '5' (90) 

For the diagram of Fig. [T5b one finds 

7r 1 m p 



m. 



where 



^ = 12^ (G^F^-' (92) 

^ = (-fJ = e^i , >v=(— ) 2 - (93) 

m 4 777. 

" 4 X i=l 

where ? and e- are the couplings of the neutralino to the relevant fermion and 

sfermion. These are calculable (see, e.g., Ref. [H]). Thus ignoring the small Yukawa 

couplings coming via the Higgsinos and taking into account only the gauge couplings, 
we find 

Z21 + tanOwZu 

C *' 5 ~ s~u^ ' l9bj 



Z 2i + (tan6 w /3)Z U 
Z\i 



€ Xi,U = —Q > 96 

sm6 w 



3cos9 w ' (97) 

where Zu, Z 2 i are the coefficients in the expansion of the B, W3 in terms of the neutralino 
mass eigenstates. Note that in this convention some of the masses m^. may be negative. 

We should mention here that, if the gluino exchange is dominant, the lepton number 
violating parameter 77 SUSY simplifies and becomes r) x > with 77 A / given in section EO 

6.2.4- The case of light intermediate neutrinos It is also possible to have light neutrino 
mediated Ou (3(3 decay originating from R-parity violating interactions. In this case one 
has the usual (3 decay vertex of the V — A type in one vertex and the sfermion mediated 
vertex, of the S — P type, in the other end [95j H99[ 1200] (see Fig. [16]). The lepton 
number violation is achieved via the mixing of isodoublet and isosinglet sfermions. The 
simplest diagram, which involves intermediate sleptons, can arise from the following 
interactions: 

Clqdc -> A m V; k L u L d R t k (98) 

k 
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Figure 16: The light neutrino mediated Oz//3/3-decay in R-parity violating SUSY models. 
In addition to the usual gauge vertex one has a scalar vertex mediated by sleptons or 
down-type squarks. The lepton number violation proceeds via the mixing between the 
isodoublet and isosinglet sfermions. This is not indicated by an x on the scalar line, as 
it is customary, since we use the sfermion mass eigenstates. 



for the d-quark vertex and 

Cllec -> A m E V£U ej D jL e c J k (99) 

for the neutrino vertex. In the above expression V L (V R ) are the mixing matrices which 
express the doublet (singlet) selectrons in terms of the mass eigenstates. U is the usual 
neutrino mixing matrix. The effective s-lepton propagator is 

^ = £%^«E^1^ *Ml = Ml-Ml (100) 

k=X 1V1 k k=l 1V1 e 

The last equation follows from the orthogonality condition on the mixing matrices and 
the fact that the splitting 6M% is small compared to the average s-lepton mass M*. If, 
further, the mixing between generations can be ignored we get 

- ™ 2 °> AM1 (101) 



2 Ml 

Combining the above results with the usual V — A coupling one gets: 

M = ^)*A& L d R e L ^^e%ii L >y x d L , (102) 

with 

Ae - ^A m A m ^ U ej 2 M , (103) 

In the case of squark exchange of Fig. [16] the above expressions become 

C LQDC -> A' ni E V d fu L e c R dl (104) 

k 

for the u-quark vertex and 

£lqd< -> A' ni E V d ^U ej u jL d R d k (105) 

k 

41 



for the neutrino vertex. Combining them we get 



with 



A< = (_^)^A^ L e c H e L7 A ^-^, u Ljx d L (106) 

in a rather obvious notation. 

In the case e and d contributions, in the context of perturbation theory, one can 
simplify the above expressions using explicitly the coupling between the singlet and the 
doublet sfermions of the lower charge. In this case 
sin 28~ 



—AM- = + A tan f3)m x x = e,d. (108) 



2 

Before proceeding farther we have to perform a Fierz transformation: 

I 

u L e c R e L ^ x k a -i a d R = - -[u L d R e L e c R k x + u L d R e L ia Xu k v e c R (109) 
+ u L ia Xu k u d R e L e c R ] - \u L ia a pd R e L ia af} e c R k x . 

o 

We must now go to the nucleon level and perform a Fourier transform to the coordinate 
space. For + — > + transitions the space component yields: 

M = (^) 2 (/a) 2 (-A [M' T + M GT + r F M;]e 7o q.7(l + 7 5 K 



^ f M GT + (1 - 2^)-^Mf ) g e(l + 7s )e c (110) 



with 



A = A d -/96 , rF = —(-2^ + l). (Ill) 



The parameter A as well as the quantities M and M have the same meaning as in the 
mass independent contribution in the conventional approach (see section I5.2j) . Note, 
however, that in the present mechanism there is no term analogous to the rj of section 
15.21 We should stress that this novel mechanism can lead to transitions J + , J ^ 0. So, 
contrary to conventional wisdom, from the observation of such transitions one cannot 
definitely infer the existence of right handed currents. 



7. Handling the short range transition operators 

We have seen that there exist many mechanisms contributing to neutrinoless double 
beta decay, involving the exchange of only heavy particles. These result to short range 
transition operators. 
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7.1. The mode involving only nucleons 



If the nucleons are treated as point like particles, then the effective transition operator 
essentially behaves like a 5 function in the inter nucleon distance. Thus their 
contribution vanishes, due to the presence of a nuclear hard core. This can be cured, 
if the nucleons can be treated as extended objects. This can be done by introducing 
into the nucleon current a nucleon form factor |96j, e.g. like a dipole shape with a 
characteristic mass m A « 850 MeV. This can also be accomplished, if one utilizes a 
quark model for the nucleon [203J. 

• V-A theories 

In this case the approach has become pretty standard, i.e. the spin isospin structure 
is similar to that for the light neutrino mass term except for the radial part 
which now becomes: 



-F N (x A ) — , x A = r ke m A , F N (x) = — [x + 3x + 3 )- 
r kl 48 v 



e 



and we are not going to discuss them further. We only mention that it can also 
proceed via the 2-pion mode with the nuclear operator associated with a 2 7r = 0.1 
(see next section) 

• S, PS theories 

In the scalar case we find that the operator is spin independent and has the same 
radial part as in the previous case. 

In the case of the pseudo scalar part we find that the operator becomes 

m e m p 3 \2m p /3J ^ 

/ 2 dF N {x A ) ( dF N {x A ) dF 2 N {x A )\\ 

o-k-o-e h 1 [x A , a k , a e ) 1 — = — 

\ x A ax A \ x A ax A ax A J J 

(112) 

where T(x A , (Je) is the tensor component. 

This process can also proceed via the 2-pion mode (see next section). In this case 
one can use the lepton number violating parameter t]^h with nuclear matrix element 
as given in the next section with c^tt = 0.20. 

7.2. The pion mode in R-parity induced Oz//3/3 decay. 

Even though the pion model, (j!2p . may be important in other cases when the 
intermediate particles are heavy, giving rise to short range operators, in this section we 
will elaborate a bit further on its application in the extraction of the R-parity violating 
parameters associated with the processes discussed above. The nuclear matrix elements 
can now be calculated using the effective transition operators 

ME k = (^) 2 a kn ™E [M& + M*-]. (113) 

m p m e 
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Figure 17: The pion mediated Oz//3/3-decay as a contact interaction (a). It arises as an 
Itt, (b) (an analogous 7r-vertex in the other nucleon is understood), and 2tt exchange 
contributions (c). The lepton number violation occurs either in one nucleon (b) or in 
the pions (c). 



Where the two above matrix elements are the usual GT and T matrix elements with 
the additional radial dependence given by 

Fgt — e ~ x > f t = (3 + 3x + x 2 ) — , (114) 

x 

Fg T =(x-2)e- x , Ffi = (1 + x) e~ x , (115) 

In the case of the pion exchange mechanism, in particular for Itt exchange, it is important 
to include the nucleon form factors [2021 1204] . We are not going to elaborate further on 
this point. The complete expressions for the transition operators are given in section |9j 
We will, instead, concentrate on the coefficients an d ai n . We will begin by 
considering the elementary particle treatment [27]. 

• The coupling coefficients a^- One finds 

X ( TYl \ 

®2n = -^glK — ( 116 ) 
6 /i \ m pJ 

Obtained under the factorization approximation [27] : 

5 

< 7r + | JpJ P \n~ > — - < 7r + | J P \0 >< 0\jp\n~ >, 
o 

<0|Jp|7r^> =mlK (117) 

with the parameter h n given by 

^ = 0.668^1 — — — (118) 

m d + m u 

Thus using the current quark masses these authors [27] find = 0.20. 

• The coupling coefficients a l7r 
One finds 

4 



1 h ( — 

36/1 T 77 \ m P, 



"itt = -Fp^T^grK [ — ) (119) 
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The needed parameters were obtained using the factorization approximation in the 
case of 1 — 7r mode 

5 

< p\JpJp\n ii > = — < p\ J P \n >< 0| Jp\n^ >, 

<p\J P \n> =F P w4.41 (120) 

The matrix element < 0|Jp|7r~ > was given above (see Eq. f )117p ). Thus these 
authors [27] find: 

a ln = -4.4 x 10~ 2 (121) 

In order to provide a check of the approximations involved in the above treatment 
and to explore the new possibilities appearing in the microscopic treatment, e.g. the role 
played by the non local terms or the three possibilities entering in Fig. [181 which are not 
distinguished in the standard treatment, we will consider the quark structure of the pion 
and the nucleon |203j . Thus we will evaluate the relevant amplitude by making a non 
relativistic expansion of the hadronic current employing a constituent quark mass equal 
to 1/3 of the nucleon mass. Furthermore for the pion and the nucleon internal relative 
quark wave functions we will employ harmonic oscillator wave functions, adjusting the 
size parameters to fit related experiments |203j . Then we will compare this amplitude 
to that obtained by more standard techniques [93] . 

a) The In mode. 

Let us begin with the second process of Eq. ( fT2l (see diagrams (b) of Fig. [TT1) . 
which is further analyzed at the quark level in Fig. [T8J In this case it is clear that 
the amplitude must be of the PS type only. The tensor contribution cannot lead to a 
pseudoscalar coupling at the nucleon level. Such a coupling is needed to be combined 
with the usual pion nucleon coupling in the other vertex to get the relevant operator for 
a + — > + decay. 

Let us begin with diagram [18] (a), which is simply the decay of the pion into two 
leptons with a simultaneous change of a neutron to a proton by the relevant nucleon 
current. Then one finds that, if the non local terms, which lead to new type operators 
not studied up to now, are ignored, the "direct" diagram makes no contribution. 

The "exchange" contribution, Fig. d8] (b), in which the produced up quark of the 
meson is not produced from the "vacuum" but comes from the initial nucleon, is a bit 
more complicated. The result is: 

c l7r = 1.37 f^(x), a l7T = 0.071 fg>{x) « 0.071, (122) 

which is in size almost a factor of 2 larger than that obtained in elementary particle 
treatment [27] (see Eq. (I121j) ) . Note, however, that our results depend on the pion size 
parameter. Finally diagram [18] (c), in which the qq is produced by the weak interaction, 
leads to the expression: 

c lw = 2 ^y^ f^x) = 3.4 xl0~ 2 ff (x) (constituent masses)(123) 



3g r Jb% 



m 
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Figure 18: The pion mediated 0^/9/9-decay in the so-called In mode. In diagram (a) 
the quarks of the pion are spectators , i.e. the heavy intermediate heavy fermion f is 
exchanged between the other two quarks. In (b) one of the interacting quarks is in the 
pion. In (c) the qq pair is produced by the weak interaction, while in (a) and (b) the qq 
is produced by the strong interaction out of the vacuum in the context of a multigluon 
exchange (a la 3po mode) indicated by x . / stands for an intermediate neutral fermion 
(heavy Majorana neutrino, gluino or neutralino). 



with 

/£(*) = x^, x= b -f (124) 

The corresponding coefficient that must multiply the nuclear matrix element for x = 2 
is 

f 2 

aiw = c l7r ^g^ = 5.0 x 1(T 3 (125) 

J A 

b). The 2tc mode. 

The first process of Eq. [12] is described by diagram (c) of Fig. [T7] and is further 
illustrated in Fig. [19] at the quark level [203 ] . In Fig. [19] / stands for an intermediate 
fermion, heavy Majorana neutrino, neutralino or gluino. 

In the case of V — A theories we find 



2 , 2 1 16 

The actual value critically depends on the harmonic oscillator size parameter. For 
= 0.5 fm we get a 27r ~ 0.1. 
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Figure 19: The Oz//3/3-decay of pions in flight (2tt mode of Fig. [T7l) illustrated at the 
quark level. / stands for a effective exchange of a heavy Majorana fermion ( heavy 
neutrino or, as in R-parity violating supersymmetry, a neutralino, gluino etc). The arcs 
around the pion line merely indicate that the pion is a bound state of two quarks. 

In supersymmetry one encounters scalar, pseudoscalar and tensor interactions. For 
the pseudoscalar interaction the situation is a bit more complicated [203]. One finds 



Taking m u = ma = (l/3)m p , 6jv = lf m and b n = 0.5fm we find that is much 
smaller than the value obtained in the elementary particle treatment [93] , i.e. we find 
«27r = —0.05. A larger value can be obtained, if one uses the current quark masses. 
For m u ~ 2md=10 MeV one finds — —1-3, but then the validity of non relativistic 
expansion may be questionable. 

8. Experimental aspects of double beta decays 

8.1. Progress of DBD experiments 

8.1.1. Experimental aspects of neutrinoless double beta decays Neutrinoless double beta 
decays are concerned with fundamental properties of neutrinos and weak interactions, 
which bear some signatures of the high energy scale, and are of great interest from 
the view point of particle physics and cosmology, as it has been described in previous 
sections. On the other hand, DBD processes are nuclear rare-decays in the low energy 
scale, which are studied experimentally by low-energy and low-background nuclear 
spectroscopy, as given in review articles [3H [62j |63j [30] . 

Since double beta decays are low-energy second-order weak processes, the decay 
rates of 2z//3/3 within the SM are of the order of 10 _20 /y, and the rates of Oz//3/3-decay 
beyond the SM are even much smaller than 2z//3/3-decay rates, depending on the assumed 




(127) 
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light neutrino mass in the case of the Majorana neutrino mediated process. Then the 
expected 0v/3(3 half-lives are of the order of Ty 2 ~ 10 27 y and 10 29 y in the case of the 
IH (inverted hierarchy) mass of 30 meV and the NH (normal hierarchy) mass of 3 meV, 
respectively. 

For experimental studies of such rare decays, large detectors with ton-scale DBD 
isotopes are needed to get Oz//3/3-decay signals in case of the IH neutrino mass. Here, 
the signals are very rare and occur as low as Epp ^2-3 MeV. Background (BG) signal 
rates, however, are huge in the energy region of Eb < 3 MeV. Thus it is crucial to build 
ultra low BG detectors to find the rare and small (V/3/3-decay signals among huge BGs 
in the low energy region. 

In spite of this, double beta decays have several unique features that make it realistic 
to search for the low-energy ultra-rare Qv(3(3 signals among huge BGs. 

i. Since (3(3 half-lives of the order of 7i/ 2 10 19 - 10 21 y are 10 orders of magnitude 
longer than the age of the earth, the (3(3 isotopes are available as almost stable isotopes 
on the earth, and it is possible to get ton-scale (10 28 ) (3(3 isotopes in order to observe 
some rare Oz//3/3-decay events with Ty 2 w 10 27 y. 

ii. There are even-even nuclei, where the double beta decays are allowed due to the 
pairing interaction, but single /3-decays are energetically forbidden. Using such DBD 
nuclei, one can be free from huge single (3 BGs, which would be larger than (3(3 rates by 
factors around 10 30 . 

iii. The 0^/3/3-decay process with a virtual Majorana neutrino exchange between 
two nucleons in a nucleus is greatly enhanced because the nucleons are close to each 
other in the nucleus. Then it is feasible to access the small neutrino masses of the orders 
of Sm slJN (solar v) - 5m ATM (atmospheric v). 

iv. High energy- resolution and/or correlation studies of (3(3 rays select low-energy 
rare (3(3 signals from huge BG events. There are several DBD nuclei to be studied to 
confirm ultra-rare Of (3 (3 '-decay events and possible DBD processes. 

Then (3(3 experiments studies neutrino properties in nuclei, which are around 10~ 15 
m in diameter. Thus the (3(3 nuclei are regarded as excellent femto (10 -15 ) laboratories 
with a large filtering power to reject single (3 and other RI BG signals and a large 
enlargement factor to enhance the Oz//3/3-decay signal of the neutrino physics interest. 
In the nuclear femto laboratory, two nucleons collide with each other. The luminosity 
is around L = 10 48 cm _2 s _1 in a single DBD isotope. Then the summed luminosity for 
one ton (10 28 ) isotopes is around 10 76 cm~ 2 s _1 . The cross section for exchange of the 
Majorana neutrinos with IH mass of 40 meV is of the order of 10~ 83 cm 2 . Then the 
event rate is around 2-3 per year. The huge luminosity of the DBD femto collider 
enables one to search for the utra-rare 0u(3(3 event and the very small neutrino mass 
09]. 

8.1.2. Progresses of DBD experiments. Progresses of DBD experiments are well 
described in review papers |3TJ [62j [631 USB 1205] and refs. therein. Here brief remarks on 
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the progresses are given. 

Early experimental studies of /3/3-decays were made by geochemical methods by 
measuring the number of the /3/3-decay products. They are inclusive measurements of 
both the 2v(3(3 and Oi//3/3-decay rates to the ground and excited states. In most realistic 
cases, the geochemistry methods give the 2z//3/3-decay rate to the ground state because 
O^/3/3-decay rates are much reduced due to the small neutrino mass and decays to excited 
states are disfavored due to the small phase space volume. 

Advanced mass spectroscopy was used to measure the number of (3(3 product 
isotopes accumulated for the long geological period of around 10 6-7 years in old ores, 
and to evaluate the long half-lives of the orders of 10 20 y. Extensive studies have been 
made on such rare-gass isotopes as 82 Kr, 128 Xe, and 130 Xe, which are /3/3-decay products 
of 82 Se 128 Te, and 130 Te, respectively, as described in the review articles. Among them 
three groups have obtained half lives of around 10 21 for 130 Te [131 EJ EEl IT4"] . 

Direct counter experiments of (3(3 decays are exclusive measurements to identify 
0i//9/9-decay signals. They have been made by coincidence counter measurements of 
two (3 rays or by two (3 measurements with tracking chambers, as described in review 
articles. 

High-sensitivity counter experiments were made by using detectors as (3(3 sources. 
Der Mateosian and Goldharber have got a limits on the 48 Ca Oz//3/3-decay rate by 
using large CaF 2 scintillators |206j . Stringent limits on the 76 Ge O^/3/3-decay rate 
were obtained with high energy-resolution Ge semiconductor detectors by Fiorini et al 
|207} 1208] . Coincidence measurements of two (3 tracks with spark and streamer chambers 
were made to get limits on the 48 Ca and 82 Se Ovfifi rays [2091 E3HI ECU EE]. 

In 1980's, the Ge experiments for 76 Ge O/3/3-decays are improved much by the 
Milano group [2T3] and Avignone et al [211] • Ejiri et al used the ELEGANT III with 
a Ge detector surrounded by Nal scintillators to get limits on the 76 Ge Oi//3/3-decays to 
the ground and excited states |215l 12161 1217] . 

The first measurement of the 21/(3(3 rays by the direct counting method was made 
on 82 Se by Elliott, Hahn and Moe [32]. The observed rate agrees with the (3(3 rate 
measured previously by the geochemical method [21811219] . The first measurement of the 
2v(3f3-dec&y rate by the direct counting method alone on 100 Mo, where no geochemical 
measurement was made beforehand, was carried out by the ELEGANT group (Ejiri et 
al.) [2201 122T] . 

So far, high-sensitivity counter experiments with the mass sensitivities of the orders 
of eV and sub eV have been made on several /3/3-decay nuclei, and half lives of 2u(3(3- 
decay rates on many nuclei have been measured by direct counting methods, as reported 
in the review papers [3U ESS 1205] . Recent experiments are discussed in 18.31 

8.2. Methods and detectors for DBD experiments 

8.2.1. Methods for DBD experiments. Double beta decays proceed normally through 
the 2z//3/3-decay process within the SM. Transition rates of the Oz//3/3-decay processes 
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beyond SM are much rarer than those of the 2z/ / 9 / 3-decay process in most cases. It 
is thus necessary to separate experimentally the 0^/9/9-decay processes from the 2u/3/3- 
decay process. 

Geochemical methods counts the number of the decay product isotopes in ores of 
DBD nuclei for geological time of TW2 ~ 10 6-7 years, which are mostly due to the 
2i/ / 9 / 9-decay process and Ozz/S/S-decay processes are not separated from the 2i//3/3-decay 
process. 

Direct counting methods have been extensively used for measuring various DBD 
processes. The Oz//3/3-decay processes are identified from the sum energy spectrum 
of E/3/3 = E(fii) + E(fo) for two rays, as shown in Fig. 2 in Ref. [30]. They 
show a sharp peak characteristic of the 2-body kinematics at Epp = Qpp, while 2v(3(} 
shows a continuum spectrum characteristic of the 4-body kinematics. Neutrinoless DBD 
followed by the Majoron (see section [5731) . which is the Goldstone boson associated with 
spontaneous breaking of B — L symmetry, shows the spectrum characteristic of the 
3-body kinematics. 

The 0i//9/3-decay processes due to the left handed weak current and the right- 
handed one (RHC) are experimentally identified by measuring the energy and angular 
correlations of the two (3 rays, as shown in Fig. 4 in ref. |30] . 

The left handed weak current Oupft process includes several modes such as the light 
v exchange, the heavy v exchange, the SUSY particle exchange, and others, as discussed 
in the previous sections. Relative contributions of these modes to the 0^/3/9-decay rate 
may be investigated by observing several DBD isotopes as well as those for the ground 
and excited states, provided that the matrix elements are properly evaluated. Then 
experimental studies of several DBD isotopes are necessary. 



8.2.2. Sensitivity of DBD experiment DBD event rates are so low that DBD 
experiments are necessarily be carried out by using high-sensitivity detectors at low- 
background underground laboratories. 

In case of the O^/3/3-decay process with the light Majorana neutrino exchange, the 
transition rate T 0u per year(y) per 1 ton (t) of DBD isotopes is expressed in terms of the 
nuclear sensitivity S n and the effective mass of the light Majorana neutrinos of \(m u )\ 
as 

r°» = \(m v )\ 2 S n . (128) 

The nuclear sensitivity is written as 

S n = (78 meV)- 2 |M^| 2 G^O.OIA)- 1 , (129) 

where M 0v is the nuclear matrix element, A is the mass number, G 0u is the phase space 
factor in units of 10~ 14 /y. 

The mass sensitivity 1(^)1 is defined as the minimum mass to be measured by the 
O^/3/3-decay experiment. It is expressed in terms of the detector sensitivity D as follows: 

\{m m )\ = S- l ' 2 D- l l\ D = {eNT){5)-\ (130) 
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where e is the Ov/3/3 peak efficiency, N is the number of the DBD isotopes in unit of 
ton, T is the run time in unit of year and S is the minimum counts required for the peak 
identification with 90 % CL (confidence level). It is given as 5 pal.Q + 1.7 (BNT) 1 ^ 2 
with B being the BG rate /t/y, and 5 pa2.3 for BNT > 2 and < 2, respectively. 

The nuclear sensitivity S n is proportional to the phase space factor G 0u and |M^| 2 . 
Thus DBD nuclei with large G Qu and large M® u are selected for the high nuclear 
sensitivity. DBD detectors with a large efficiency e, a large amount (N) of isotopes, 
and a small BG rate (B) are used for the high detector sensitivity. 

8.2.3. DBD detectors Neutrinoless f3~f3' decays of (A, Z) -» (A,Z + 2) are studied 
by measuring two f3~ rays, while neutrinoless f3 + f3 + decays of (A, Z) — > (A, Z — 2) 
proceed through (3 + (3 + , (3 + EC, and EC EC 7, where EC (electron capture) is detected 
by measuring the X ray, and f3 + by measuring the /3 + and the two 511 keV annihilation 
7 rays. 

High-sensitivity DBD experiments require DBD isotopes with high nuclear 
sensitivity S n , i.e. the large Qpp and the large phase space factor G 0u , as given by 
Eq. (EEZHD - 

There are several such DBD isotopes of decays. However, no nuclei 

with large S n are available, although BG rates for the (3 + (3 + and (3 + EC are quite small 
by measuring both /3 + and the annihilation 7 rays and/or the X ray. Accordingly most 
of high-sensitivity DBD experiments are concentrated on the Oi//3~/3~ decays with the 
large Qpp and G 0u . Thus hereafter, we discuss mainly (3~(3~ decays. The large Qpp ~3 
MeV helps reduce BG rates since most of BGs from natural RIs are below 3 MeV. 

The 0vj3~j3~ experiments with the IH (30 meV) mass sensitivity are carried out by 
using low BG (B pal/t y) detectors and ton-scale (N ~0.5 - 1) DBD isotopes with the 
large nuclear matrix element of \M® U \ pa3 (see Tabled]) and the large phase space factor 
of G 0u ~ 5 in unit of 10 _14 /y for a long (T w 2 - 4 year) run time. On the other hand 
one needs DBD isotopes of around iV w 50 — 100 ton and ultra-low BG detectors with 
B pa 0.01/t/y to reach the NH (2-4 meV) mass sensitivity. 

Possible DBD isotopes to be used for high-sensitivity 0uf3~f3~ experiments are given 
in Tabled Among them, 82 Se, 100 Mo, 116 Cd, 130 Te, 136 Xe have the large Q w w 3 MeV 
and the large G 0u ~ 5 in unit of 10~ 14 /y, as shown in Table [H 130 Te has the largest 
abundance ratio of 34.5%, while others have the abundance ratios of around 10%, and 
are enriched to 85-90% by means of centrifugal separation. 

76 Ge is a special case with smaller Qpp pa 2 MeV and the smaller G 0l/ ~ 0.71, 
but high sensitivity experiments are possible by using low-BG 76 Ge detectors with high 
energy-resolution [203 [208]. 150 Nd has a very large phase space factor of G 0u ~ 23.2 
in units of 10 _14 /y, but the natural abundance ratio is only 5.6%, and the enrichment 
is hard. The matrix element may be reduced because of the large difference in nuclear 
shape between the initial and final nuclei. 

DBD experiments are carried out by using either calorimetric detectors or 
spectroscopic detectors. Calorimetric detectors are made partly of DBD isotopes, and 
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thus the detection efficiency is as large as e ~ 0.6-0.9. Cryogenic detectors with Te [222] , 
semiconductor detectors with 76 Ge and 116 Cd, and scintillation detectors with 116 Cd and 
136 Xe are used as calorimetric detectors. 

Cryogenic bolometers of ZnS, CaMo04, ZnMo04 and others are shown to be used 
as high energy-resolution DBD detectors. BG rates of these detectors are reduced by 
measuring both thermal signals as well as scintillation signals [2231 12241 225J. 

Spectroscopic tracking detectors with DBD isotopes outside the detectors are used 
to measure individual /3 rays from the isotopes. Then the low BG measurements are 
possible even though the efficiency is low and the energy resolution is modest. 82 Se 
and 100 Mo are studied by spectroscopic detectors with PL (plastic scintillator) arrays. 
Spectroscopic detectors measure (3 — (3 energy and angular correlations, which are used 
to confirm and identify the 0^/9/9-decay process. 

So far, high-sensitivity DBD experiments have been made mainly on the ground 
state + transition because of the large phase space factor. Transitions to the excited 
+ states are experimentally measured in coincidence with 7 rays. Then BG rates are 
much reduced, even though the phase space factors are smaller by almost one order of 
magnitude. Measurements of both the ground and excited transitions are interesting 
to confirm the 0u(3(3— decay signal and to study the 0v(3/3 mechanism. The excited + 
states in 82 Se, 100 Mo, 136 Xe, and 150 Nd are located in low excitation region. Some excited 
states are different in shape from the ground state. Thus O^/3/3-decay measurements for 
both the ground and excited states are very interesting. 

8.3. Present status and future projects of DBD experiments 

8.3.1. Neutrinoless double beta decays Experimental studies of O^/3/3-decay have been 
carried out on several (3(3 nuclei [30] |3T| |4*9] . Some of them by counter measurments 
are listed in Table [U 128 Te was studied by a geochemical method, and the hallife and 
mass limits are 7.7 ± 10 24 y and 1.1 - 1.5 eV |226j . Here the effective mass given in 
the 6th column shows a range of the values evaluated from the experimental half-life by 
using various matrix elements. In fact the mass depends much on the matrix element as 
discussed in section [TUJ Calorimetric detectors have been used for isotopes such as 48 Ca, 
76 Ge, 116 Cd and 130 Te. Among them, 76 Ge experiments (Heidelberg-Moscow, IGEX) 
with Ge semiconductors [361 [223 [2281 USE [230] and the 130 Te experiment (CUORICINO) 
with 41 kg Te02 cryogenic bolometers |231[ 1232] have good energy-resolution and give 
stringent limits on the absolute value of effective Majorana neutrino mass of the order 
of 0.3 - 0.5 eV. Currently, the most stringent limit on \ (m u ) \ comes from the lower limit 
on the T 1 ° / 1/ 2 ( 136 Xe) measured in KamLAND-Zen experiment [40] (see Tabled]). 

Spectroscopic detectors (ELEGANT V, NEMO III) have been used for 82 Se, 100 Mo, 
116 Cd and other isotopes with large Qpp values [38] I233[ I234[ 1235] . They are /3-ray 
tracking detectors with (3(3 sources separated from detectors. NEMO III provides 
stringent limits on the Qv(3(3 half-lives for 82 Se, 100 Mo, and other isotopes [381 12341 1235] . 

Recently, a claim for the Oi/S/3-decay peak, corresponding to the effective Majorana 
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Table 1: Limits on neutrino less double beta decays Ty^ exp (claim for evidence is denoted 
with upper index c). Qpp : Q value for the + — > + ground state transition. G 0u : 
kinematical (phase space volume) factor (g^ = 1.25 and R = 1.2 fm A 1 / 3 ). (m„): The 
upper limit on the effective Majorana neutrino mass, deduced from T®J^ exp by assuming 
the ISM [236] (g e J f = 1.25, UCOM src), the EDF [131] (g e J f = 1.25, UCOM src), the 
(R)QRPA (1.00 < g e j[ < 1.25, the modern self-consistent treatment of src), and the 
IBM-2 [130] (1.00 < g e J f < 1.25, Miller-Spencer src), nuclear matrix elements (see 
section [TUj) . src means short-range correlations. 
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neutrino mass of 0.32 eV, was made by a part of the Heidelberg-Moscow collaboration 
[4"T| |4*2] . The result depends on the off-line analysis method. It should be checked by 
future GERDA/MAJORANA experiments with lower BG rates. 

Neutrino-mass sensitivities of the CUORICINO and NEMO III detectors are limited 
to be around 300 - 500 meV because of the limited (3(3 isotopes of 11 - 7 kg and the 
large BG rates. Thus future experiments with higher mass sensitivity are necessary to 
prove/disprove the Heidelberg- Moscow claim and to search for the Majorana neutrino 
in the lower v mass regions. 

The neutrino oscillation studies have given a strong impact to high-sensitivity 
studies of (3(3 experiments since the effective mass suggested is of the order of ^5m 2 ~ 
2 meV - 50 meV, which next-generation (3(3 detectors can access if the z/'s are Majorana 
particles. Future experiments with higher mass sensitivities are in progress (see below). 

8.3.2. Two neutrino double beta decays The 2z//3/3-decay is a process fully consistent 
with the SM of electroweak interaction. The inverse half-life of 2i//3/3-decay is free of 
unknown parameters on the particle physics. It can be expressed as a product of an 
accurately known phase-space factor G 2u (E , Z), which includes fourth power of gA, and 
the double Gamow Teller transition matrix element M 2u (A, Z), which is a quantity of 
the second order in the perturbation theory: 



Here M 2u includes nuclear effects due to the nuclear residual interactions, the 
nuclear medium and the nuclear quenching. It is obtained from the experimental decay 
rate. Experimental studies of 2z//3/3-decay half-lives have been made on some nuclei by 
geochemical methods, and several nuclei by direct counting methods |241j . The 2i/(3(3- 
decay has been observed so far in 12 nuclides ( 48 Ca, 76 Ge 82 Se, 96 Zr, 100 Mo, 116 Cd, 
128 Te, 130 Te, 136 Xe, 150 Nd, 130 Ba and 238 U) and in two excited states [21]. Recent 
NEMO III experiments provide high-statistic spectroscopic studies of the 2i//3/3-decay 
rates [38 ll240j . Spectroscopic measurements of two (3 rays are useful to reduce BG rates 
and energy correlations of two (3 rays are used to identify the 2v(3(3 mechanism. 

The measurement of 2i//3/3-decay rates gives us information about the product of the 
squared effective axial-vector coupling constant and 2z//3/3-decay matrix elements. They 
are presented in Table [2j The 2z//3/3-decay matrix elements are sensitive to nuclear-spin 
isospin correlations. The observed values for M 2v are used to investigate the nuclear 
structure and the nuclear interactions associated with the 0^/9/3-decays. 

8.3.3. High sensitivity experiments In the case of the inverted mass hierarchy, (3(3 
detectors with the IH mass sensitivity of (m m ) ^ 20 - 50 meV can be used to study 
the Ov (3 /3-decay, while in the case of the normal hierarchy one needs higher sensitivity 
detectors with (m m ) ^2-4 meV. Several groups are now working for next-generation 
(3(3 experiments with the IH mass sensitivities of 20~50 meV, as discussed in the reviews 
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Table 2: The 2v/3(3 matrix elements \M 2u \ deduced from the measured half-life Ty 2 by 
counter experiments [241] [4"0] . g& = 1.269 is assumed. 
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Experimental proposals for future f3f3 experiments have been made on several f3/3 
isotopes, and some of them are listed in Table [H They are mostly experiments 
because of large kinematical (phase space) factors. DBD experiments with different 
isotopes and different methods are indispensable to confirm and identify the 0^/9/3-decay 
event and the O^/3/3-decay mechanism. Some of them are briefly described below. 

76 Ge experiments with low-BG high resolution 76 Ge detectors are of special interest 
for proving or disproving the Heidelberg-Moscow claim of the large O^/9/3-decay peak 
[4"2] . and for further high-sensitivity ton-scale experiments. 

GERDA: It aims at high energy-resolution studies of 76 Ge O^/3/3-decays by using 
high-purity 76 Ge detectors to check the Heidelberg-Moscow claim and the possible 
Majorana neutrinos in the QD region at LNGS (Gran Sasso). GERDA uses 18 kg 
76 Ge detectors in Phase I, and add 20 kg detectors in Phase II. The Ge detectors are 
immersed into high purity liquid Ar in order to avoid BG contributions from cryostats 
[175] . GERDA is now running. 

MAJORANA: The MAJORANA demonstrator uses 40 kg Ge detectors at the 
Sanford underground lab. to test /investigate the half life of the Heidelberg-Moscow 
claim and the QD mass regions and to prove the feasibility for a future ton-scale IH 
(10 27 y, 20-40 meV) experiment. The Ge detectors are PPC (P-type Point Contact) 
detectors with excellent PSA(Pulse Shape Analysis). They are cooled by using ultra- 
pure electro-formed Cu cryostat [179J. The BG goal is 4/t y, which scales to 1/t y for 
the 1-ton experiment. The enriched detectors will be on-line in 2013 - 2014. 
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These detectors can also be used to study DM and neutrino scattering in the low 
energy region. The GERDA and the Majorana collaboration will be merged for one 
ton-scale future experiment by selecting the best techniques developed and tested by 
GERDA and MAJORANA. Recent developments are given in the report [242] . 

There are other experimental plans for QD masses. Among them, CANDLES is for 
48 Ca /3/3-decays with an array of CaF2 crystals |243j . which is based on the ELEGANT 
VI experiment with CaF 2 . The Qpp is large, but the natural abundance of 48 Ca is only 
0.2%. Thus the efficient isotope enrichment is crucial. 

Several groups are working for future experiments with the IH mass sensitivities of 
around (m u ) ^ 20 - 50 meV. DBD isotopes required are those with the large = 
2.5 - 3 MeV and G 0u = 3-5 10~ 14 /y to get large nuclear sensitivities of the order of 
(SV) -1 / 2 =15 - 20 meV. The experiments use large-scale low-BG detectors with ton-scale 
enriched isotopes and B wl/ty. 

MOON (Molybdenum Observatory Of Neutrinos): This is an extension of 
ELEGANT V [233]. It is a hybrid (3(3 and solar v experiment with 100 Mo to study 
the Majorana v masses with the QD - IH mass sensitivities of 100 - 20 meV and the low 
energy solar us |244[ I245j I246j I247j 1248] . Double beta decays to both the ground and 
the 1.132 MeV excited + states are studied to confirm the Ov(3(3 events and to study 
the Oi//3/3-decay mechanism. MOON can be used for supernova neutrinos as well |249] . 

Detectors under considerations are : A) the super-module of PL plate and fiber 
scintillators |245[ I246[ 1247] and B) the cryogenic bolometer of ZnMoO^ The PL 
scintillator module (A) is used for spectroscopic study of two /3-ray energy and angular 
correlations to identify the 0z/ (3 /3-decay process. Here one module consists of a plate (PL) 
scintillator for the (3 energy and two sets of X-Y fiber scintillator planes for the vertex 
identification, between which a thin 100 Mo film is interleaved. The energy resolution is 
o~ ^2.2% at E = Qpp to reduce the 2z//3/3-decay tail in the O^/3/3-decay window. The half 
life (mass) sensitivity is 3 x 10 26 y (45 meV) with 480 kg 100 Mo for 5 years. Proto-type 
detectors were built to show the energy resolution as required |245[ [246J . 

The ZnMo04 bolometer (B), which is under discussion with the Milano-Rome group 
is for calorimetric study of the sum of two /3-ray energy. The high energy-resolution 
(AE ~ 5 keV), the particle identification by the scintillation and/or the pulse-shape 
analysis and the high efficiency(e ~0.8) make the high-sensitivity study possible. Thus 
it is good to start with the detector B, and proceed to A to confirm the 0z//3 (3- decay 
process by two /3-ray measurement. The half-life (mass) sensitivities are 4 x 10 25 y (120 
meV) for 3 y run with 12 kg 100 Mo and 7 x 10 26 y (30 meV) for 5 y run with 220 kg 
100 Mo. 

SuperNEMO: The goal of SuperNEMO is to reach a sensitivity of 10 26 y, which 
corresponds to the IH mass of 40-110 meV [250] • The detector consists of huge tracking 
chambers and scintillation detectors with 100 kg of (3(3 isotopes of 150 Nd or 82 Se to 
search for the v mass below 0.1 eV. The efficiency is 30% and the resolution is 4% in 
FWHM. The BG impurities are deduced to be 208 T1 < 2 and 214 Bi < 10 in unit of /i 
Bq/kg It plans to use 20 modules, each module with 5 kg (3(3 isotopes. The first module 
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is a demonstrator with 7 kg of 82 Se. The mass sensitivity of the demonstrator with 15 
kg y is 210-570 meV, while that of the full detector array with 500 kg y is 53-145 meV, 
depending on the nuclear matrix element. It is based on NEMOIII, and thus it is crucial 
to improve the energy resolution and the efficiency 

AMoRE Advanced Molybdenum based Rare process Experiment Large volume 
CaMoC>4 crystals with enriched material have been developed to study the 0^/3/3-decays 
of 100 Mo and to search for cold dark matter |251j . Pilot experiments of 1 kg with 
scintillation technique and Cs(I) active veto are in preparation. In order to improve the 
energy resolution, cryogenic CaMo0 4 detectors are being developed. To avoid BGs from 
the 2i//3/3-decays of 48 Ca, depletion of Ca in 48 Ca < 0.001% is made by using ALSIS 
(Advanced Laser Stable Isotope Separation). Additional light sensor and time constant 
of phonon signal are effective to select signals. The goal of AMoRE is to study 0^/3/3 
decays of 100 Mo in the region of IH mass of 50 meV (3 1 26 y ) by using 100 kg CaMo04 
cryogenic detectors. 

COBRA: It uses a large amount of high energy-resolution CdZnTe semiconductors 
at room temperature |252j . The modular design makes coincidence measurements 
possible to reduce BG rates. The crystal includes several (3(3 isotopes to be studied. 
The collaboration now tests 64 CZT 1 cm 3 detectors at LNGS. The goal is to study the 
Majorana neutrino in the IH mass region by measuring the 0^(3(3 from 116 Cd with Qfts 
= 2.809 MeV. The detector is composed by 64 K crystals with 0.42 ton Cd isotopes 
enriched in 116 Cd. Reduction of BGs from RI impurities inside and around detectors 
are important. Pixelisation(Semiconductor tracker) can be a major step forward. 

CUORE (Cryogenic Underground Observatory for Rare Event): This is an 
expansion of CUORICINO. It is a high energy resolution bolometer array to measure the 
130 Te OvPP decays with Q = 2.529 MeV at LNGS [239]. It uses natural Te0 2 crystals 
with the natural 130 Te abundance of 34%. It consists of 988 Te02 crystals with the net 
130 Te mass of 203 kg. The detector array is under construction since 2005. 

The experiment emphases 20 times more massive than CUORETINO, better energy 
resolution of 5 keV, high granularity, and thus low BG rates. The CUORICINO BG 
rate, which is around 0.16/keV/kg/y (4 10 3 /5keV/y/ton of 130 Te), is expected to be 
reduced to 0.02 - 0.01/keV/kg/y in CUORE. Then, in cases of the BG rates of B = 0.01 
- 0.001/keV/kg/y (2.5 10 2 - 2.5 1075keV/y/ton of 130 Te), the half-life and the v mass 
sensitivities are expected to be around 2.1 - 6.5 10 26 y and 50 - 25 meV, which depend 
on the nuclear matrix element. The first CUORE tower is CUORE-0. 

EXO (Enriched Xenon Observatory): The (3(3 experiment of 136 Xe with Q = 
2.467 MeV [253J is made by using the laser tagging technique to select the residual 
nuclei of 136 Ba to suppresses all kinds of RI BGs. The energy-resolution of around 
a ~ 2 % is achieved by measuring both the ionization and scintillation signals. The 
lton enriched Xe detector with the energy-resolution of o = 1.6% gives the i/-mass 
sensitivity of 50-70 meV for a 5 y run. The 10 ton Xe detector with the improved energy 
resolution of a = 1 % will give the sensitivitjjJJ] of 11-15 meV. The 200 kg 136 Xe liquid 

|| In a recent paper the EXO collaboration[254! reported that no signal has appeared in a search for 
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detector is used at WIPP to study the 2^/3/3-decay and the quasi-degenerate z/-mass, as 
the first step without the Ba tagging. EXO observed the 2u(3(3 half life of 2.1 x 10 21 
[255] . The key point of this experiment is the tagging efficiency of the 136 Ba nuclei. 

KamLAND-Zen (Kamioka Large Anti Neutrino Detector Zenon): It studies the 
136 Xe DBD by means of the KamLAND detector with the 1 kton liquid scintillator 
at Kamioka [256] . A mini balloon with 3.2 m in diameter is set at the center for the 
136 Xe-loaded liquid scintillator. It includes 136 Xe isotopes around 400 kg. The energy- 
resolution and the vertex-resolution are 6.8 %/y/E and 12.5 cm/s/E. The collaboration 
neasured the 2v(3(3 half-life [40J, and hopes to reach the sensitivity of around 50 meV. 

NEXT (Neutrino Experiment with a Xe TPC): A Xe TPC with 100-150 kg enriched 
136 Xe is used at LSC [257J. It is a low BG and good E-resolution TPC with separate 
readout planes for tracking and energy. The NEXT- 100 sensitivity for 5 y run is about 
5.9 10 25 y (better than 100 meV). 

Borexino with 2 ton 136 Xe: The sensitivity is around 100 meV [258J. 

DCBA (Drift Chamber Beta-ray Observatory): It uses a tracking chamber in a 
magnetic field to study 150 Nd (3(3 decays [259.J. The (3 energy is obtained by the /3-ray 
trajectories. DCBA-T3 is now under construction. The good energy-resolution and 
efficient enrichment of 150 Nd isotopes are necessary. 

SNO+ (Sudbury Neutrino Observatory +): It uses the 1 k ton scintillation detector 
with 0.1 % Nd isotopes to study QD-IH v masses by using natural (5.6 %) and enriched 
(50%) 150 Nd isotopes [260J. The mass sensitivities are 100 meV with the natural (5.6 
% 150 Nd 56 kg) and 40 meV with enriched isotopes (50 % enriched 150 Nd 500 kg). The 
collaboration is trying to find a realistic way of the Nd isotope separation, which is of 
great interest to study the IH neutrino mass. It aims at the scintillator filling at the 
beginning of 2013. 

DBD experiments for NH mass: Higher sensitivity DBD experiments for NH masses 
of {m v ) =2-4 meV require a large amount of high nuclear-sensitivity (S n ~15 meV) 
DBD isotopes of the order of iV f» 10 — 50 tons, and extremely low BG detectors with 
B <0.1-0.01/t y. DBD isotopes to be studied are 82 Se, 100 Mo, and 136 Xe. High energy- 
resolution cryogenic detectors such as ZnSe, ZnMo04 with pulse shape analyzes and/or 
scintillation signals |223j 1224] may be used to search for the NH masses by 0z//3,9-decay 
experiments of 82 Se and 100 Mo. 

8.4- Experimental studies of DBD matrix elements 

8.4-1- Experimental probes for DBD matrix elements 

Nuclear matrix elements (Mj") for O^/3/3-decay are crucial for extracting the 
effective Majorana v mass and other parameters, relevant to particle physics models 

neutrinoless double-beta decay of 136 Xe with an exposure of 32.5 kg-yr and a background of w 1.5 x 10 -3 
kg _1 keV^ 1 y~ 1 in the ±er region of interest. This implies a lower limit on the half life, > 1.6 x 10 25 

y (90% CL), corresponding to an effective Majorana mass of less than 140-380 meV, depending on the 
nuclear matrix element. 
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Figure 20: Nuclear spin isospin responses for weak interactions. They are studied by v 
probes via weak interactions, EM (7) probes via EM interactions and by nuclear probes 
via strong interactions jl9], [301 H32| 



beyond the SM, from 0i>(3(3 experiments, while nuclear matrix elements (M 2u ) for 2vf3f3- 
decay can be derived experimentally from the observed 2z//3/3-decay half-lives. Extensive 
calculations of M^ T and M® u have been made in terms of QRPA, RQRPA, shell model, 
and so on, as given elsewhere in the theoretical sections. 

Most (3(3 strengths are located in (3(3 (double Gamow- Teller) giant resonances, i.e., 
in the high-excitation region [261J . Thus the (3(3 matrix elements get very small in 
comparison with single particle estimates and are sensitive to nuclear structures, nuclear 
spin-isospin correlations, nuclear deformations, nuclear interactions, nuclear medium 
effects on the weak coupling constant qa and others. The theoretical evaluations for 
them are hard. Experimental studies of nuclear structures and nuclear interactions, 
which are relevant to 2v(3(3-dec&y and O^/3/3-decay matrix elements, are very interesting 
to get reliable evaluations for them gSJ EQJ ElEJ 1262] . 

Nuclear matrix elements of M 2ugt and M® u are expressed in terms of the successive 
single (3 processes through intermediate \J n ) states. Among the intermediate states, 
low-lying single particle- hole states play dominant roles |261[ 1263] . The single (3 matrix 
elements are given by spin-isospin responses for Qtslj = r ± [i L r L Y[ J x a s }. They are 
studied experimentally using hadron, photon, and lepton probes as shown in Fig. [20j 

Lepton probes 

Nuclear weak responses for low-lying intermediate states are obtained from single 
(3 decay rates and EC rates. However, they are limited to /3 ± decays from the ground 
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state in the intermediate nucleus. The decays are mostly allowed Gamow- Teller (GT) 
decays with r^cr, and are first forbidden decays with r ± irY 1 and T ± [ir L Yi x o~ s ] in some 
nuclei. 

Muon capture reactions of are used to get the (3 + strengths in the 

intermediate nucleus [2621 1264] . Excitation energies and angular momenta involved 
in this reaction are E— 0-50 MeV and O 1 * 1 , l^ 1 and 2^. The capture cross section is 
quite large, and most muons are stopped and captured in the case of medium heavy 
nuclei. Then the r + weak strength distribution in the intermediate nucleus is derived by 
measuring the decaying neutrons and 7-rays from excited states and radioactive isotopes 
produced by the muon capture. 

One direct way to get the weak responses is to use v beams. Since v nuclear 
cross sections are as small as a = io~ 40 ~ 42 cm 2 , one needs high-flux v beams and large 
detectors [4"9| l30| 1132] . Low energy v beams with E < 100 MeV can be obtained from 
pion decays. Intense pions are produced by nuclear interaction with GeV protons. Weak 
decays of stopped 7r + give low energy neutrinos as 

7T+ ->■ fi + + fi + -> u e + + e + , (132) 

where and u e , are well separated by the decay time. Intense 1 GeV protons from 
SNS at ORNL provides intense neutrinos around 10 15 per second |265] and the J-PARC 
booster synchrotron with 3 GeV protons produces neutrinos around 3xl0 14 per second 

[266]. 

Photon probes 

Weak responses for (3 + decays are studies by using photo-nuclear (7,X) reactions 
through isobaric analogue states (IAS) as shown for the first time by (p/y) reactions 
[267[ 1268] . The /3 and 7 matrix elements are related as 

< f\g v mP\i >^ — (2T) 1 ' 2 < f\e m^IAS >, (133) 

e 

where \IAS >= (2T)- l ' 2 T-\i >, T is the isospin of the parent nucleus and mP 1 and m 7 
are analogous /3 and 7 transition operators. Thus one can obtain the /3 matrix element 
for \i >— > \f > by observing the analogous 7 absorption |/ >— > \IAS > through the 
IAS of \i >, where \f > and \i > are the final state and the intermediate state in the 
/3/3-decay, as shown in Fig. [2TJ These photo-nuclear reactions through IAS are used to 
get the /3 + matrix elements to excited states in the intermediate nucleus. 

In medium heavy nuclei, IAS is located on the El giant resonance (GR) at the high 
excitation region. Accordingly IAS shows up as a sharp isobaric analogue resonance 
(IAR), and the photo- nuclear reaction includes IAR, GR, and the interference term as 
given by |269] 

dai ] Q N) = k[\A{I)j,\ 2 + Ej\A(G)j\ 2 + 2Re(A(I)j, A{G)j, j% (134) 

where A(I)jr,A(G)ji and are the IAR and GR amplitudes and the relative phase at 
IAR. Then one can get the phase of the matrix element from the interference. 
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Figure 21: Level and transition schemes of (3/3, single (3 + , and 7 via IAS. T is the 
isospin lowering operator. 

Laser electron photons, which are obtained from laser photons scattered off GeV 
electrons, are used for the photonuclear reaction. The polarization of the photon can be 
used to study El and Ml matrix elements separately [262] . HI7S (High Intensity 7-ray 
Source), New SUBARU and other electron synchrotrons with E e =l - 3 GeV provide 
laser electron photons used for the photo-nuclear reaction. 

Nucleon and nuclear probes 

Nuclear charge exchange reactions with nuclear (hadron) probes are used to study 
nuclear spin-isospin responses. Charge exchange spin-nonflip reactions are used for 
vector weak responses, while charge exchange spin-flip reactions for axial-vector weak 
responses [30l 126 lj . Extensive studies of charge exchange reactions have been made to 
get GT(1 + ) responses with r^a. The reactions studied are (p,n) (n,p), (d, 2 He), ( 3 He,t) 
(t, 3 He), and ( 7 Li, 7 Be) at IUCF, KVI, MSU, RCNP, Triumf and others. Some charged 
particle experiments are in |270[ I271[ I272[ I273[ 1274] . Medium-energy projectiles with 
Ei/A= 0.1 - 0.3 GeV are used for studying rcr responses because of the relatively large 
spin-isospin interaction (V Tcr ) and the small distortion interaction (Vo) at the medium 
energy. 

The charge exchange reaction with the medium-energy projectile is mainly due 
to the central isospin and spin-isospin interactions. Then the cross section with the 
transferred momentum (q) and energy (oj) is given as 



where K(E i ,cu), , J a , and B(a) are the kinematical factor, the nuclear distortion 




(135) 
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\J a \ 2 B(a), a = F, GT. (136) 



factor, the volume integral of the spin-isospin interaction, and the nuclear spin-isospin 
response, respectively, a denotes the isospin and spin channel; a=F for isospin Fermi 
and a= GT for spin-isospin GT. 

The cross section at deg for q ps and u ~ is corrected for the kinematical and 
normalization factors, and is expressed as 

da a (0°) 1 
dVt K(E i} 0)ND 

The proportionality of the cross section at the forward angle of 9 ~0 deg to the ra 
response B(a) has been studied for charge exchange reactions with medium-energy light 
projectiles. In fact, the proportionality is good for spin-flip reactions with B{GT) > 0.1, 
but some deviation from the proportionality is found in the reactions with smaller 
B(GT) due to contributions from tensor/non-central interactions. 

Charge exchange ( 3 He,t) reactions (ChER) relevant to the /3/3-decays were studied 
at RCNP by using the 420 MeV 3 He beam and the high energy-resolution beam line and 
spectrometer system [2701 EZH [2751 EZSl [2771 EZSl (2791 [2801 12SH- The beam stability 
and the beam line system have been improved to give the fantastic energy-resolution 
of AE/E f» 5 x 10~ 5 and AE ps 25 keV, and ChERs have been studied for several (3/3 
nuclei. 

The GT strengths are mostly located at the GR resonances in the high excitation 
region, and accordingly, the strengths for the low-lying states are small. In cases of 
96 Zr, 100 Mo, and 116 Cd, where valence neutrons and valence protons are in the different 
major shells, there is only one strong GT state (ground state) of (g7/2) n (g9/2) p . Thus 
the ( 3 He, t) reaction shows a strong GT strength only at the ground state. The high 
resolution GT^ strength measurement can give significant insight into the details of the 
nuclear structure and can help to determine the O^/3/3-decay NME less nuclear model 
dependent. 

Nucleon transfer reactions were studied to get the occupation and vacancy 
probabilities of nucleons involved in double beta decays [282] . They are used to evaluate 
j3/3 matrix elements. 

8.4-2. DBD matrix elements via low lying intermediate states In this subsestion, 
we briefly discuss experimental GT strengths for low-lying states and 2v(3(3 matrix 
elements [1 32 1 [30] . The matrix element M 2u is given approximately by a sum of products 
of two single GT matrix elements, Mqj, (m) from the direction of (3 + and Mq? (m) from 
/3~. It is expressed as M 2v ps Mq u Tj and 

M 2u = V M^(m)M^(m) 

GT ^Qw/2 + m e + E x (l+)-E ' 1 1 

Here, (_E X (1+) — E ) is the energy difference between the m th intermediate 1 + state and 
the initial ground state. Qpp is the Q-value of the /3/3-decay. 

The single GT matrix elements can be derived from charge exchange ( 3 He, t), 
(d, 2 He) and other reactions and /^/EC-decay rates. Measurements of GT ± strengths 
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have been performed for A = 100, 116 [2701 EH], 48 [2711 ESS], 76 [277], A = 96 
[2721 1279] . 116 [275], A=136 [280J and 150 [281]. For A = 128 and 130 the data are 
expected to be available soon. 

The coherent sum of the weighted products of the measured M c t T (m) and Mqj, (m) 
are in accord with the observed M^ [270| 1279] . Actually it has been shown that the 
2vf3f3 proceeds through low-lying (Fermi Surface) states, and not much through the GT 
giant resonance [1321 [263] . in agreements with observations by charge exchange reactions. 
The measured 2vf3f3 half-lives and the matrix elements [24 lj are given in Table [21 

SSD (Sigle state dominance) hypothesis [283] suggests that 2z//3/3-decay transitions, 
where the first 1 + state of the intermediate nucleus is the ground state (e.g., 100 Mo, 
116 Cd, 128 Te) are governed by the transition through that J n = 1 + state. The SSD 
hypothesis looks fine in case of nuclei like 100 Mo, where there is only one J 71 " = 1 + state. 

Recently the nuclear matrix elements M G V T are shown to be expressed in terms of 
single (3 matrix elements via Fermi-surface quasi particle states (FSQP) [263[ I284[ I285J. 



M% T = Y,M{k){A{k))-\ M(k) = M { G ^(k)M { G ^(k), (138) 

k 

where the sum is over the FSQP (low-lying) states (k) in the intermediate nucleus, 
A(fc) is the energy denominator and the matrix elements of M^(fc) and M G ^(k) are 
the experimental single GT matrix elements deduced from charge exchange reactions, 
f3 decay rates and EC rates. Thus it includes no adjustable parameters. The possible 
deviation of qa from 1.26 is embedded in the observed matrix elements of Mq^ and 
similarly in M$(k),M$(k). 

In the quasi-particle representation, the experimental matrix elements are given by 

Mct (A;) = k e l h m(j k , J k )P t {k), M { G ^(k) = kf s m( Jk , J k )P f (k), (139) 

where m(j k J n k) is the single particle matrix element for the GT transition between 
single particle j k and J k states with j k and J k being the neutron and proton spins, and 
Pi{k) = U p (J k )V n (j k ) and Pf{k) = V p (J k )U n (j k ) are the pairing reduction factors for 
transitions from the initial ground state to k-th intermediate state and from this state to 
final state, respectively. The effective coupling constants, k e J^ and k A ^ f , in unit of g A , 
represent the nuclear core effects such as the spin-isospin correlations, the short-range 
correlations, and others, in addition to the quenching effect, while the nuclear surface 
(shell) effects are given by the pairing factors of Pi(k) and Pf{k). In fact values for k A 
and k A do not depend much on individual states, and thus one can evaluate the single 
(3 matrix elements, if not available experimentally, by using the k e A ^ and k A i for other 
states in neighboring nuclei and the caluculated values for m(j k , J k )Pij(k). 

As shown in Fig. [22} the observed 2z/ / 9 / 9-decay matrix elements are well reproduced 
by the sum of the matrix elements through the low-lying FSQP 1 + states in intermediate 
nucleus. 

In general there are several low-lying 1 + states, and thus 2z//3/3-decay proceeds not 
only through the lowest state but also other FSQP states. The product of the matrix 
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Figure 22: Nuclear matrix elements for 2v{3/3 decays . Open squares: observed values 
(M 2u ) , closed squares: FSQP values (M = EM (k) /A (k)). The matrix elements for the 
excited 0+ state in 100 Ru are plotted at ,4=102 [28411285] . 



elements is M$(k) M^(k) = k e J fi k e J f) (m(j k , J k )) 2 Pi{k)P f {k), which is positive. 
Thus contributions from the FSQP states are not cancelled with each other, but are 
constructive. 

The 0^/3/3 transition operator is a two-body operator. It is shown that the matrix 
element is approximately given by the sum of the matrix elements for successive 
processes via intermediate states [301 [132J . Then the matrix elements for the successive 
(3 transitions via low-lying intermediate states are derived from the single (3 matrix 
elements as in case of the 2v/3/3 process. 

Since the 0i//9/9-decay process is a virtual neutrino exchange between 2 nucleons 
in the nucleus, intermediate states with higher spins are involved. They are studied 
by measuring angular distributions of charge exchange reactions to intermediate 1*, 2 ± , 
states and photo-excitations to IAS of the intermediate l^ 1 states [262J. Charge exchange 
reactions and photo-excitations of IAS are under progress by Muenster, MSU, NC, 
RCNP, and other groups to study nuclear structures relevant to 2v(3/3 and 0v/3{3 
processes. 
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8.5. Two-neutrino double beta decay and bosonic neutrino 

Neutrinos may possibly violate the spin-statistics theorem, and hence obey Bose 
statistics or mixed statistics despite having spin half. A violation of the spin-statistics 
relation for neutrinos would lead to a number of observable effects in cosmology and 
astrophysics. In particular, bosonic neutrinos might compose all or a part of the cold 
cosmological dark matter (through bosonic condensate of neutrinos) and simultaneously 
provide some hot dark matter. A change of neutrino statistics would have an impact 
on the evolution of supernovae and on the spectra of supernova neutrinos. The idea of 
bosonic neutrinos has been proposed independently in Ref. |286j . where cosmological 
and astrophysical consequences of this hypothesis have been studied. 

If neutrinos obey at least partly the Bose-Einstein statistics the Pauli exclusion 
principle (PEP) is violated for neutrinos. The parameter sin 2 x can be introduce to 
characterize the bosonic (symmetric) fraction of the neutrino wave function |287J. A 
smooth change of sin 2 \ from to 1 transforms fermionic neutrinos into bosonic ones. 
The assumption of violation of the PEP leads to a number of fundamental problems 
which include loss of a positive defmiteness of energy, violation of the CPT invariance, 
and possibly, of the Lorentz invariance as well as of the unitarity of S-matrix. No 
satisfactory and consistent mechanism of the Pauli exclusion principle violation has 
been proposed so far. 

The lepton number conserving 2i//3/3-decay) can be used to study the statistical 
properties of neutrinos |287j . In the presence of bosonic neutrinos two contributions 
to the amplitude of the decay from diagram with permuted neutrino momenta have 
relative plus sign instead of minus in the Fermi-Dirac case. 

The PEP violation strongly changes the rates of the 2^/3/3-decays and modifies the 
energy and angular distributions of the emitted electrons. The effect of bosonic neutrinos 
is different for transitions to + ground states and 2 + excited states. In Fig. [23]the energy 
spectra of two electrons for different values of the bosonic-fraction sin 2 \ is presented 
for the 2t/ / 9 / 9-decay of 100 Mo to ground state of final nucleus. With increase of sin 2 x the 
spectra shift to smaller energies. We note that substantial shift occurs only when sin 2 \ 
is close to 1.0 Pure bosonic neutrinos are excluded by the present data. |287j . In the 
case of partly bosonic (or mixed-statistics) neutrinos the analysis of the existing data 
allows to put the conservative upper bound sin 2 x < 0.6 [287] . 

9. Effective transition operators 

The subject of interest are the lepton number violating (LNV) parameters associated 
with the exchange of light and heavy Majorana neutrinos and with R-parity breaking 
SUSY mechanisms. 

By assuming the dominance of a single mechanism determined by the LNV 
parameter r] K the inverse value of the Oz//3/3-decay half-life for a given isotope (A, Z) 
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T [MeV] 

Figure 23: (Color online) The differential decay rates normalized to the total decay rate 
vs. the sum of the kinetic energy of outgoing electrons T for 2i//3/3-decay of 100 Mo to 
the ground state of final nucleus. The results are presented for different values of the 
squared admixture sin 2 \ of the bosonic component. The spectra have been calculated 
in the single state dominance approximation [288] 1289] . 



can be written as 

\ Vk \ 2 \m ,v ;\ 2 g^(e ,z). (i io) 



\ _ \ n 1 2 |fl / f' 0,/ |2 /Mi/', 

T?Y. 



•1/2 

Here, G 0u (E ,Z) and M'°* are, respectively, the known phase-space factor (E is the 
energy release) and the nuclear matrix element, which depends on the nuclear structure 
of the particular isotopes (A, Z), (A, Z + 1) and (A, Z + 2) under study. 

The phase space factor G 0u (E , Z) includes fourth power of unquenched axial- vector 
coupling constant g& and the inverse square of the nuclear radius R~ 2 , compensated by 
the factor R in M' ^ . The assumed value of the nuclear radius is R = r^A 1 ^ with 
tq = 1.1 fm or tq = 1.2 fm in different publications. The implicit radius and g& 
dependencies in G 0u (E Q , Z) and nuclear matrix element and the problem of the correct 
use of them were discussed in [290J. 

The nuclear matrix element M' 0u is defined as 



eff 



2 



M ,0v = M o,_ ( m) 

Here, g e j[ is the quenched axial-vector coupling constant. This definition of M' ^ 
[29 lj 1292] allows to display the effects of uncertainties in g e j[* and to use the same 
phase factor G 0u (Eq, Z) when calculating the 0i//3/3-decay rate. 
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Before we proceed with the discussion of the nuclear matrix elements we will 
summarize the various types of transition operators entering the neutrinoless double 
beta decay process. We recall that in the various particle models the lightest particle 
exchanged between the two nucleons participating in this process is either much lighter 
than the electron or much heavier than the proton. It is thus possible to separate the 
particle physics parameters from those of nuclear physics. Furthermore the nature of this 
exchanged particle dictates the form of the transition operators. The LNV parameters 
of interest together with corresponding nuclear matrix elements are presented briefly 
below. 



9.1. Transition operators resulting from light neutrino exchange 
Leutri 



In the case of light-neutrino mass mechanism of the O^/3/3-decay we have 

i ni„) 
m e 



(142) 



The nuclear matrix element associated with the light Majorana-neutrino exchange 
M® u consists of the Fermi (F), Gamow- Teller (GT) and tensor (T) parts as 



Ml 



e .ff\2 



(o+| E 

kl 



+ M, 



GT 



M% 



H F {r kl ) 

" (g e A ff ) 2 



+ H GT (r kl )a M - H T (r k i)S k i]\0p. 



Here 

S k i = 3(<7fc • r k i)(ai ■ r k i) — a k i, o ki = a k ■ 
The radial parts of the exchange potentials are 



2 r c 

HF,GT,T(r k l) = —R I 
7T Jo 



jofi,2{q r ki)h F ,GT,T{ ( l 2 ) ( l 



dq. 



(143) 
(144) 

(145) 



q + E 

where R is the nuclear radius and E is the average energy of the virtual intermediate 
states used in the closure approximation. The closure approximation is adopted in the 
calculation of the NMEs relevant for Oi//3/3-decay with the exception of the QRPA. The 
functions hp,GT,T(q 2 ) are given by [293J 



h F (q 2 ) =f 2 v {q 2 ) 



hcriq 2 ) 



(/i P -/in) 2 q 2 
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(146) 



For the normalized to unity vector and axial-vector form factors the usual dipole 
approximation is adopted: f v (q 2 ) = 1/(1 + q 2 /M^) 2 , f A (q 2 ) = 1/(1 + q 2 /M\) 2 . M v = 
850 MeV, and Ma = 1086 MeV. qa = 1-254 is assumed and the difference in magnetic 
moments of proton and neutron is — /i n ) = 4.71. 

The above definition of the includes contribution of the higher order terms of 
the nucleon current and for the induced pseudoscalar the Goldberger-Treiman PCAC 
relation, gp(q 2 ) = 2m p gA(q 2 )/{q 2 + m 2 ) was employed [101] . 

9.2. Transition operators resulting the heavy neutrino exchange mechanism 

We assume that the neutrino mass spectrum include heavy Majorana states N k with 
masses M k much larger than the typical energy scale of the Oz//3/3-decay. These heavy 
states can mediate this process as the previous light neutrino exchange mechanism. 
The difference is that the neutrino propagators in this case can be contracted to points 
and, therefore, the corresponding effective transition operators are local unlike in the 
light neutrino exchange mechanism with long range internucleon interactions. The 
corresponding LNV parameters are rjj^ and 77^. 

Separating the Fermi (F), Gamow- Teller (GT) and the tensor (T) contributions we 
write down for the NME 

_ M f 

"JV 



U 0v _ 1V1 F(N) M M 
•M „ = 2 r MgT(N) - M T {N) 



9 a 

= (OtlErkWl^P 1 + H { G N T \r kl )a kl - H?\r kl )S kl )\0+), 
ki 9a 

(147) 

where S k i and a k i are given in Eq. (j!44p . The radial parts of the exchange potentials 
are 

HfA TT (r k i) = / j 0fi ,2{qrki)h FyGT {q 2 )q 2 dq. (148) 

7r m p m e Jo 

9.3. Transition operators resulting from the R-parity breaking SUSY mechanism 

Assuming the dominance of gluino exchange, we obtain for the LNV parameter the 
following simplified expression 



% a s X' 2 n m p 
6 G 2 p m 4 7 m 5 



1+P^ 

\ m UL. 



(149) 



Here, Gp is the Fermi constant, a s = <?f/(47r) is SU(3) C gauge coupling constant. m Si , 
mj^ and m- g are masses of the u-squark, d-squark and gluino, respectively. 

We should mention again that for this type of interactions the pion-exchange 
mechanism (Itt and 2n) discussed in subsection 17.21 dominates over the conventional 
two-nucleon mechanism. Thus, denoting the Oz//3/3-decay NME associated with gluino 
and neutralino exchange as M°y , we have |28l [93] 

M°J = c lT (M% T - M^) + c 2T (m 2 J t - M 2 ?) (150) 
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with 

^ 2 y^K 9,F P ^ 1 flm% g 2 s _ 

9 m z p m e {m u + m d ) g\ ' 18 mpn e (m u + m d ) 2 g 2 A 

Here, gs and Fp stand for the standard pion-nucleon coupling constant (g s = 13.4) and 
the nucleon pseudoscalar constant (we take the bag model value Fp 4.41 from Ref. 
[204J), respectively. f n = 0.668771^ and m n is the mass of pion. m u and denote 
current quark masses. The partial nuclear matrix elements of the Tji v SUSY mechanism 
for the Qv/3/3 process are: 
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E rjtr? 
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(152) 



with 



H&, T {r u )= --Rj Q jofl fr H ) ^±L^ f*tf)dq, 
H 2 G t, t (tm) =--R joMrki) ( ,l'\, 2 f A {^)dq. (153) 



1 + q 2 jm\ 
4 D i" 00 , N <r7 m t 

The two-nucleon exchange potentials are expressed in momentum space as the 
momentum dependence of normalized to unity nucleon formfactors (/a(? 2 )) is taken 
into account. 



9.4- Transition operators resulting from the squark-neutrino mechanism 

In the case of squark-neutrino mechanism [202] due to the chiral structure of the R-parity 
breaking ($ p ) SUSY interactions, the amplitude of 0z//3/3-decay does not vanish in the 
limit of zero neutrino mass unlike the ordinary Majorana neutrino exchange mechanism 
proportional to the light neutrino mass. Instead, the squark-neutrino mechanism is 
roughly proportional to the momentum of the virtual neutrino which is of the order 
of the Fermi momentum of the nucleons inside of nucleus pp ~ lOOMeV. This is a 
manifestation of the fact that the LNV necessary for Oi//3/3-decay is supplied by the 
]fip SUSY interactions instead of the Majorana neutrino mass term and therefore this 
mechanism is not suppressed by the small neutrino mass. The corresponding SUSY 
LNV parameter is defined as 



^ilfe^ifcl „ ; _ nQ d 

V m l(fc) m \(k)l 

(154) 



q V 2y/2G F {k) 1 ™ 2 ™ 2 



Here we use the notation da.) — d,s,b. This LNV parameter vanishes in the absence of 
q~L — qn - mixing when 9 d = 0. 

At the hadron level we assume dominance of the pion-exchange mode. Then, the 
nuclear matrix element associated with squark-neutrino mechanism can be written as a 
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sum of GT ad tensor contributions [21)2] 



Mf = M Gm -M m . 
The exchange potentials are given by 



(155) 
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10. Oz//3/3-decay nuclear matrix elements 

Interpreting existing results as a measurement of 1(^)1 and planning new experiments 
depends crucially on the knowledge of the corresponding nuclear matrix elements 
(NMEs) that govern the decay rate. The NMEs for Oz//3/3-decay must be evaluated 
using tools of nuclear structure theory. There are no observables that could be directly 
linked to the magnitude of Oz//3/3-decay nuclear matrix elements and, thus, could be used 
to determine them in an essentially model independent way. A reliable calculation of 
NMEs will be of help in predicting which are the most favorable nuclides to be employed 
for O^/5/3-decay searches. 

The evaluation of the nuclear matrix elements can be separated in two steps. 

• The evaluation of the transition matrix elements between the two interacting 
particles (two body ME). 

Each particle is assumed to occupy a set of single particle states, determined by the 
assumed model. The spin as well the orbital structure of the operator as it has been 
discussed in sections [5] and El The operators discussed in section [5] are long ranged, 
except when the intermediate neutrino is heavy leading to short range transition 
operators. In section El except for the case of Eq. (11 101) . all operators are short 
ranged. The way of dealing with the short range operators has been discussed in 
section [71 Taking all these into account the effective transition operators have been 
constructed in section [91 

• The second step involves the construction of the many body wave functions. 

One needs the wave functions of the initial and final nuclear systems. If closure 
is not employed, as in the case of QRPA, one also needs the wave functions of 
all the virtual (intermediate states) allowed by the assumed nuclear model. Some 
many body features arising from the nuclear medium are: i) the renormalization 
effects on the coupling and the modification of the nucleon currents, which have 
already been discussed in section [9] and ii) the short range correlations, which will 
be discussed below. The main techniques of the construction of the many body 
wave functions will be reviewed in this section. We remind the reader that in some 
cases information about the nuclear ME can be extracted from experiments (see 
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The calculation of the O^/3/3-decay NMEs is a difficult problem because ground and 
many excited states (if closure approximation is not adopted) of open-shell nuclei with 
complicated nuclear structure have to be considered. In the last few years the reliability 
of the calculations has greatly improved. Five different many-body approximate meth- 
ods have been applied for the calculation of the O^/3/3-decay NME: 

(i) The Interacting Shell Model (ISM) [2361 12931 [2951 [2961 [297] . 

The ISM allows to consider only a limited number of orbits close to the Fermi level, 
but all possible correlations within the space are included. Proton-proton, neutron- 
neutron and proton-neutron (isovector and isoscalar) pairing correlations in the 
valence space are treated exactly. Proton and neutron numbers are conserved and 
angular momentum conservation is preserved. Multiple correlations are properly 
described in the laboratory frame. The effective interactions are constructed 
starting from monopole corrected G matrices, which are further adjusted to describe 
sets of experimental energy levels. The Strasbourg-Madrid codes can deal with 
problems involving basis of 10 11 Slater determinants, using relatively modest 
computational resources. A good spectroscopy for parent and daughter nuclei 
is achieved. Due to the significant progress in shell-model configuration mixing 
approaches, there are now calculations performed with these methods for several 
nuclei. 

(ii) Quasiparticle Random Phase Approximation (QRPA) [115[ 1116] . 

The QRPA has the advantage of large valence space but is not able to comprise 
all the possible configurations. Usually, single particle states in the Wood-Saxon 
potential are considered. One is able to include to each orbit in the QRPA model 
space also the spin-orbit partner, which guarantees that the Ikeda sum rule is 
fulfilled. This is crucial to describe correctly the Gamow- Teller strength. The 
proton-proton and neutron-neutron pairings are considered. They are treated in the 
BCS approximation. Thus, proton and neutron numbers are not exactly conserved. 
The many-body correlations are treated at the RPA level within the quasiboson 
approximation. Two-body G-matrix elements, derived from realistic one-boson 
exchange potentials within the Brueckner theory, are used for the determination of 
nuclear wave functions. 

(iii) Interacting Boson Model (IBM) [130]. 

In the IBM the low lying states of the nucleus are modeled in terms of bosons. 
The bosons are in either L=0 (s boson) or L=2 (d boson) states. Thus, one is 
restricted to + and 2 + neutron pairs transferring into two protons. The bosons 
interact through one- and two-body forces giving rise to bosonic wave functions. 

(iv) The Projected Hartree-Fock-Bogoliubov Method (PHFB) [129] , 

In the PHFB wave functions of good particle number and angular momentum 
are obtained by projection on the axially symmetric intrinsic HFB states. In 
applications to the calculation of the O^/3/5-decay NMEs the nuclear Hamiltonian 
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was restricted only to quadrupole interaction. The PHFB is restricted in its scope. 
With a real Bogoliubov transformation without parity mixing one can describe only 
neutron pairs with even angular momentum and positive parity. 

(v) The Energy Density Functional Method (EDF) [131] . 

The EDF is considered to be an improvement with respect to the PHFB. The 
density functional methods based on the Gogny functional are taken into account. 
The particle number and angular momentum projection for mother and daughter 
nuclei is performed and configuration mixing within the generating coordinate 
method is included. A large single particle basis (11 major oscillator shells) is 
considered. Results are obtained for all nuclei of experimental interest. 

The differences among the listed methods of NME calculations for the O^/3/3-decay 
are due to the following reasons: 

(i) The mean field is used in different ways. As a result, single particle occupancies of 
individual orbits of various methods differ significantly from each other [298J . 

(ii) The residual interactions are of various origin and renormalized in different ways. 

(iii) Various sizes of the model space are taken into account. 

(iv) Different many-body approximations are used in the diagonalization of the nuclear 
Hamiltonian. 

Each of the applied methods has some advantages and drawbacks, whose effect in the 
values of the NME can be sometimes explored. The advantage of the ISM calculations is 
their full treatment of the nuclear correlations, which tends to diminish the NMEs. On 
the contrary, the QRPA, the EDF, and the IBM underestimate the multipole correlations 
in different ways and tend to overestimate the NMEs. The drawback of the ISM the 
limited number of orbits in the valence space and as a consequence the violation of Ikeda 
sum rule and underestimation of the NMEs. 

In Table |3j recent results of the different methods are summarized. The presented 
numbers have been obtained with the unquenched value of the axial coupling constant 
(d a = 9a — 1-25)0, Miller- Spencer Jastrow short-range correlations |299j (the EDF 
values are multiplied by 0.80 in order to account for the difference between the unitary 
correlation operator method (UCOM) and the Jastrow approach [293] ). the same 
nucleon dipole form-factors, higher order corrections to the nucleon current and the 
nuclear radius R = r Q A 1 ^, with r = 1.2 fm (the QRPA values [29TT [292] for r = 1.1 
fm are rescaled with the factor 1.2/1.1). Thus, the discrepancies among the results of 
different approaches are solely related to the approximations on which a given nuclear 
many-body method is based. 

From Table [3J we can make the following conclusions: 

(i) The ISM values of NMEs, with the exception of the NME for the double magic 
nucleus 48 Ca, practically do not depend on the nucleus. They are significantly 
smaller, by about a factor 2-3, when compared with NMEs of other approaches. 

% A modern value of the axial- vector coupling constant is qa — 1.269. We note that in the referred 
calculations of the O^/3/3-decay NMEs the previously accepted value g A = Qa = 1-25 was assumed. 
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Table 3: The NME of the O^/3/3-decay M® v calculated in the framework of different 
approaches: interacting shell model (ISM) |236[ 1297] . quasiparticle random phase 
approximation (QRPA) [H3 M [3001 EQD [302], projected Hartree-Fock Bogoliubov 
approach (PHFB, PQQ2 parametrization) ) |129j . energy density functional method 
(EDF) [131J and interacting boson model (IBM) [130]. QRPA(TBC) and QRPA(J) 
denote QRPA results of Tuebingen-Bratislava-Caltech (TBC) and Jyvaskyla (J) groups, 
respectively. The Miller-Spencer Jastrow two-nucleon short-range correlations are taken 
into account. The EDF results are multiplied by 0.80 in order to account for the 
difference between UCOM and Jastrow [293] . g e j[ = 9a = 1-25 and R = 1.2 fmA 1 ^ 3 are 
assumed. 
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(ii) The largest values of NME are obtained in the IBM ( 76 Ge and 128 Te), PHFB ( 100 Mo, 
130 Te and 150 Nd), QRPA ( 150 Nd) and EDF ( 48 Ca, 96 Zr, 116 Cd, 124 Sn and 136 Xe) 
approaches. 

(iii) NMEs obtained by the QRPA (TBC) and IBM methods are in a good agreement 
(with the exception of 150 Nd). 

(iv) In the case of 130 Te all discussed methods, with the exception of the ISM, give 
practically the same result. 

(v) The disagreement between IBM-2 and ISM is particularly troublesome, because 
IBM-2 is a truncation of the shell-model space to the S and D pair space and, in 
the limit of spherical nuclei, IBM-2 and ISM should produce the same results. 

(vi) The disagreement between the QRPA(TBC) and QRPA(J) results is not large but 
it needs to be clarified. 

Comparing O^/3/3-decay nuclear matrix elements calculated by different methods 
gives some insight in the advantages or disadvantages of different candidate nuclei. 



73 



10 



29 



10 



28 



> CJ 10 
o 

H 



27 



10 



26 



10 



25 



















▼ (R)QRPA : 


















A LSSM 


- 
















• EDF 


















■ 


1 PHFB 


















♦ IBM-2(I) : 


A 

_ 


A 












• 

v 






- 
- 




V 


▲ 


T 






▲ 






A 






♦ T 














A 


T 




_ 




• 


■ 




T 

• 


• 






• 


• ; 








• 














♦i 










■ 


































■t 










100_ _ 


116_ 


124_ 


128^ 








48 Ca 


76 Ge 


82 Se 


Zr 


Mo Cd Sn 


Te 


13 °Te 


136 Xe 


150 Nd 



Figure 24: (Color online.) The Oz//3/3-decay half-lives of nuclei of experimental interest 
for 1(771^)1 = 50 meV and NMEs of different approaches. The Miller-Spencer Jastrow 
two-nucleon short-range correlations are considered. The axial-vector coupling constant 
qa is assumed to be 1.25. 



However, matrix elements are not quite the only relevant quantities (see section [8] for 
the nuclear sensitivity factor). Experimentally, half-lives are measured or constrained, 
and the effective Majorana neutrino mass (m„) is the ultimate goal. For | | equal to 
50 meV the calculated half-lives for double /3-decaying nuclei of interest are presented 
in Fig. |2H We see that the spread of half-lives for given isotope is up to the factor of 
4-5. 

It is worth to noticing that due to the theoretical efforts made over the last years the 
disagreement among different NMEs is now much less severe than it was about a decade 
before. Nevertheless the present-day situation with the calculation of O^/3/3-decay NMEs 
can not be considered as completely satisfactory. Further progress is required and it 
is believed that the situation will be improved with time. Accurate determination of 
the NMEs, and a realistic estimate of their uncertainty, is of great importance. Nuclear 
matrix elements need to be evaluated with uncertainty of less than 30% to establish the 
neutrino mass spectrum and CP violating phases of the neutrino mixing 

10.1. Uncertainties in calculated NMEs 

The improvement of the calculation of the 0^/9/9-decay NMEs is a very important and 
challenging problem. The uncertainty associated with the calculation of the 0u/3/3- 
decay NMEs can be diminished by suitably chosen nuclear probes. Complementary 
experimental information from related processes like charge-exchange and particle 
transfer reactions, muon capture and charged current (anti) neutrino- nucleus reactions is 
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very relevant. A direct confrontation of nuclear structure models with data from these 
processes improve quality of nuclear structure models (see section [8]). The constrained 
parameter space of nuclear models is a promising way to reduce uncertainty in the 
calculated Oz//3/3-decay NMEs. 

A steady progress in nuclear structure approaches is gradually leading to a better 
understanding and to a reduction of the differences among their results. However, 
even in the most refined approaches, the estimates of remain affected by various 
uncertainties, whose reduction is of great importance. 

i) QRPA calculation of NMEs. 

Due to its simplicity the QRPA is a popular technique to calculate the O^/3/3-decay 
NMEs. One of the most important factors of the QRPA calculation of the (V/3/3-decay 
NMEs is the way the particle-particle strength of the nuclear Hamiltonian g pp is fixed. 
The Tuebingen-Bratislava-Caltech (TBC) group has shown that by adjusting g pp to the 
2i//3/3-decay rates the uncertainty associated with variations in QRPA calculations of the 
Oi//3/3-decay NMEs can be significantly eliminated [2911 EES [2931 US] . In particular, the 
results obtained in this way are essentially independent of the size of the basis, the form 
of different realistic nucleon-nucleon potentials, or on whether QRPA or renormalized 
QRPA (take into account Pauli exclusion principle) is used. This new way of fixing 
parameter space was criticized by the Jyvaskyla group in series of papers maintaining 
the role of single /3-transitions. It was claimed that careful study of single /3 and 2^/3/3 
observables points to serious shortcomings of adopted procedure [304] I305| . These 
objections were refuted in [291] 1292] . In the recent publications [300, 1301] 1302] also 
the Jyvaskyla group adopted the procedure of fixing of g pp proposed by the TBC group 
[303] . 

Usually, two variants of the QRPA are considered. The standard QRPA, which 
is based on the quasiboson approximation, and the renormalized QRPA (RQRPA) 
|127] 1174] 1291] . which takes into account the Pauli exclusion principle. Further 
improvement is achieved within the self-consistent QRPA (SRQRPA) [298] 1306] 1307] 
by conserving the mean particle number in the correlated ground state. The restoration 
of Pauli exclusion principle and of particle number conservation lead to a reduction of 
the Oz//3/3-decay NMEs pH [29H [2921 [298] . 

There is some controversy about the importance of the tensor Mt contribution to 
M® u . According to the ISM [236] the tensor term is small, a fact understood by the 
small model space adopted. The Jyvaskyla group performing the calculation within the 
QRPA claims that M T is negligible [3001 EDI]. M T was neglected in the PHFB [129] 
and EDF [131] calculations of the Oz//3/3-decay NMEs. Contrary, results of the IBM-2 
|130] and the QRPA (TBC) [101] calculations show that Mt can not be neglected and 
its absolute value can be up to 10% of M® v . 

ii) The closure approximation 

The Oz//3/3-decay matrix elements are usually calculated using the closure approximation 
for intermediate nuclear states. Within this approximation energies of intermediate 
states (E n — E{) are replaced by an average value E « 10 MeV, and the sum 
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over intermediate states is taken by closure, J2n \ n >< n \ = 1- This simplifies 
the numerical calculation drastically. The calculations with exact treatment of the 
energies of the intermediate nucleus were achieved within the QRPA-like methods 
[174[ 12911 I292[ 12931 1303] . The effect of the closure approximation was studied in detail 
in [308J. It was found that the differences in nuclear matrix elements are within 10%. 
The the dependence of the NMEs on the average energy of the intermediate states E 
was studied within the nuclear shell model. By varying E from 2.5 to 12.5 MeV the 
variation in the NME was obtained to be less than 5% |297] . 

Hi) The two-nucleon short range correlations and finite nucleon size. 
The physics of finite nucleon size (FNS) and two-nucleon short-range correlations (SRC) 
is different. Both reduce magnitude of the 0i//9/9-decay NME by competing each with 
other. The importance each of them depends on the type of SRC and involved form- 
factor parameters. 

The FNS is taken into account via momentum dependence of the nucleon form- 
factors. For the vector, weak-magnetism (axial-vector) the usual dipole approximation 
with cut-off parameter My = 850 MeV (Ma = 1 086 MeV), which comes come from 
electron scattering experiments (neutrino charged- current scattering experiments), is 
considered. The form-factors suppress high-momentum exchange. We note that in 
the limit of point-like nucleon (M^a ~* oo) the weak magnetism contribution to the 
Oz//3/3-decay would be divergent. 

The SRCs are included via the correlation function f(r), that modifies the relative 
two-nucleon wave functions at short distances: 



Previously, Miller-Spencer Jastrow SRC (a = 1.1 fm~ 2 , b = 0.68 fm~ 2 , c = 1.0) 
have been added into the involved two-body transition matrix elements, changing two 
neutrons into two protons, to achieve healing of the correlated wave functions. A 
suppression of M® u by 20% was found |130|. 11744 1236} 1297] . However, recent work has 
questioned this prescription [174L 13001 1301] . 

The two-nucleon short range correlations were studied within the coupled clusters 
method (CCM) in [174] . The Jastrow function fit for T=l channel reported for Argonne 
V18 and Bonn-CD NN interactions set of parameters (a = 1.59 fm~ 2 , b = 1.45 fm~ 2 , 
c = 0.92) and (a = 1.52 /m" 2 , b = 1.88 fm~ 2 } c = 0.46), respectively [174]. These 
correlation functions were confirmed by exploiting a construction of an effective shell 
model operator for Oz//3/3-decay of 82 Se [309J. The notable differences between the results 
calculated with Miller- Spencer Jastrow and CCM SRC are about of 20%-30% [TO 1297] . 
The previous results with Miller-Spencer treatment of SRC certainly overestimates the 
quenching due to short-range correlations. Of course, the results obtained with the CCM 
SRC are preferable. We note that in the case of Bonn-CD CCM SRC the Oi//3/3-decay 
NMEs are slightly increased [174] . 



®ni(r) -> [1 + f(r)]V nl (r), 
where f(r) can be parametrized as [174] 
f(r) = - ce - ar2 (l-br 2 ). 



(158) 



(159) 
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The two-nucleon short range correlations were treated also through the Unitary- 
Correlation Operator Method (UCOM) [HI HEEl USD], which has the advantages of 
wave-function overlap preservation and a range of successful applications [31 lj . The 
drawback of this approach, when applied to the (V/3/3-decay, is that it violates some 
general properties of the Fermi and Gamow- Teller matrix elements [174] . 

Recently, a question of many-body short-range correlations in the evaluation of the 
0^/3/9-decay NMEs was addressed within a simple model [312J. The existing calculations 
include long-range many-body correlations in model dependent (ISM, QRPA, PHFB, 
EDF, IBM) nuclear wave functions but allow only two particles to be correlated at short 
distances. There are some indications that it is not sufficient. 

iv) The effect of deformation. 

The nuclei undergoing double beta decay, which are of experimental interest, are 
spherical or weakly deformed nuclei with exception of 150 Nd, which is strongly deformed. 
It was found in |313j that deformation introduces a mechanism of suppression of the 
2i//3/3-decay matrix element which gets stronger when deformations of the initial and 
final nuclei differ from each other |313j 1314] . A similar dependence of the suppression 
of both M 2v and M® v matrix elements on the difference in deformations has been found 
in the PHFB jHSJ [316] and the ISM [236]. The NMEs have a well-defined maximum 
when the deformations of parent and daughter nuclei are similar, and they are quite 
suppressed when the difference in the deformations is large. The ISM results suggest 
that a large mismatch of deformation can reduce the matrix elements by factors as large 
as 23. Within the IBM-2 the effects of the deformation are introduced through the 
bosonic neutron-proton quadrupole interaction. For weakly deformed nuclei the effect 
is a reduction by about 20%. 

The QRPA calculation of the O^/3/3-decay NMEs requires a construction of all 
states of the intermediate nucleus, even if closure approximation is considered. The 
results were obtained in spherical limit, which is a significant simplication. Recently, the 
proton-neutron deformed QRPA with a realistic residual NN interaction was developed 
PH 15171 13T81 13T9] . This approach was applied in the case of 76 Ge, 150 Nd and 160 Gd 
and lead to the conclusion that the O^/^-decay of 150 Nd, to be measured soon by the 
SNO+ collaboration, provides one of the best probes of the Majorana neutrino mass 
[TMf3"T9] . 

v) The occupancies of individual orbits. 

The occupancies of valence neutron and proton orbits determined experimentally 
represent important constraints for nuclear models used in the evaluation of the 0^/3/3- 
decay NME. For the 76 Ge and 76 Se they have been extracted by accurate measurements 
of one nucleon adding and removing transfer reactions by J. Schiffer and collaborators 
[320[ 1321] . The main motivation to study these nuclei was the fact that they are the 
initial and final states of 0^/9/5-decay transitions. These measurements offer a possibility 
to compare these experimental results with the theoretical occupations and, if necessary, 
detect which modifications would be required in the mean field or the effective interaction 
in order to obtain improved agreement with the experiment. 
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In a theoretical study [298J measured proton and neutron occupancies were used 
as a guideline for a modification of the effective mean field energies, which resulted 
in a better description of these quantities. The calculation of the (V/3/3-decay NME 
for 7e Ge performed with an adjusted Woods-Saxon mean field combined with the self- 
consistent RQRPA (SRQRPA) method |298j . which conserves the mean particle number 
in correlated ground state, led to a reduction of M® u by 20%-30% when compared to 
the previous QRPA values. 

In the ISM the variation of the nuclear matrix element (NME) for Oz//3/3-decay 
of 76 Ge was studied after the wave functions were constrained to reproduce the 
experimental occupancies of the two nuclei involved in the transition. It was found [322J 
that in the ISM description the value of the NME is enhanced about 15% compared 
to previous calculations. This diminishes the discrepancies between the ISM and the 
QRPA approaches. 

The role of occupancies of the single-particle orbitals in the standard QRPA 
calculation of M® u to ground and 0+ states of final nucleus were studied also in [3231 1324] . 
Unlike the treatment of |298] , whereby occupancies in respect to the correlated SRQRPA 
ground state were considered, the occupancies were evaluated at the level of uncorrelated 
BCS ground state. The basic features of the ground and excited state decays were found 
to be quite different. 

vi) The axial-vector coupling constant qa- 
It is well known that the calculated strengths of GamowTeller /5-decay transitions to 
individual final states are significantly larger than the experimental ones. That effect 
is known as the axial-vector current matrix elements quenching. To account for this, 
it is customary to quench the calculated GT matrix elements up to 70%. Formally, 
this is accomplished by replacing the true value of the coupling constant Qa = 1.269 
(the previously Qa = 1.254 was considered) by a quenched value Qa = 1-0. The 
origin of the quenching is not completely known. This effect is assigned to the A-isobar 
admixture in the nuclear wave function or to the shift of the GT strength to higher 
excitation energies due to the short-range tensor correlations. It is not clarified yet 
whether similar phenomenon exists for other multipoles, besides J = 1 + . 

Quenching is very important for the double beta decay because appears to the 
fourth power in the decay rate. If it occurs also for the O^/3/3-decay, it could significantly 
reduce the Oz//3/3-decay half-life by as much as a factor of 2-3. The axial- vector coupling 
constant g*j[ or m other words, the treatment of quenching, is also a source of differences 
in the calculated Oi//3/3-decay NMEs. M' 0u is a function of squared g\ , which appears 
by vector and weak- magnetism terms following the definition of Eq. (U4ip . 

In |325] three independent lifetime data (2i/ / 9 / 9-decay, EC, /3-decay) were accurately 
reproduced in the QRPA by means of two free parameters (g pp , g A ), resulting in an 
overconstrained parameter space. The general trend in favor of g 6 ^ < 1 was confirmed. 
This novel possibility to reconcile QRPA results with experimental data, which deserves 
further discussions and tests, warrants a reconsideration of the quenching problem from 
a new perspective. 
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As it was manifested above nuclear NMEs for Oi//3/3-decay are affected by relatively 
large theoretical uncertainties. Within the QRPA approach, it was shown that, within 
a given set of nuclei, the correlations among NME errors are as important as their size 
[326J. This represents a first attempt to quantify the covariance matrix of the NMEs, 
and to understand its effects in the comparison of current and prospective Oz//3/3-decay 
results for two or more nuclei. It would be useful if other theoretical groups in the 
Ovfifi field could present "statistical samples" of NME calculations as well, in order to 
provide independent estimates of (co)variances for their NME estimates. A covariance 
analysis like the one proposed [326] represents a useful tool to estimate correctly current 
or prospective sensitivities to effective Majorana neutrino mass {m v ). 

Currently, the uncertainty in calculated Oz//3/3-decay NMEs can be estimated up to 
factor of 2 or 3 depending on the considered isotope, mostly due to differences between 
the ISM results and the results of other approaches (QRPA, PHFB, EDF, IBM) and 
also due to unknown value of ■ 

10.2. Anatomy of NMEs 

The anatomy of the O^/3/3-decay NME was performed in [2361 |29_3] . was decomposed 
on the angular momenta and parities J"* of the pairs of neutrons that are transformed 
into protons with the same J" 71 ' . It was found that the final value of reflects two 
competing forces: the like particle pairing interaction that leads to the smearing of 
Fermi levels and the residual neutron-proton interaction that, through ground state 
correlations, admixes "broken-pair" (higher-seniority) states. The function C 0u (r) that 
describes the dependence of the on internucleon distances r, 

Jo 

was subject of interest. It was shown that the above competition implies that only 
internucleon distances r < 2 — 3 fm contribute to M®" [293]. The maximum value 
of C 0u (r) occurs around r = 1 fm, which means that almost the complete value of 
M® u comes from contribution of decaying nucleons which are close to each other. This 
distance correspond to a neutrino momentum of q ~ 200 MeV, twice larger value as 
was expected before. This finding, which explains a small spread of results for different 
nuclei, was confirmed also by the ISM [236] and a similar behavior for C 0u (r) was 
obtained also within the PHFB p9]. The QRPA and ISM functions C 0u (r) differ only 
by a scaling factor, which is expected to be related with the ratio of the average number 
of pairs in both calculations. 

The largest component of M® u is the GT part. We have 

M / = M°J t (1 + X f + Xt), (161) 

where xf and \gt are matrix element ratios that are smaller than unity and, 
presumably, less dependent on the details of the applied nuclear model. In [327] it 
was shown that M(¥ T is related to the closure 2i//3/3-decay NME M^. That relation is 
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revealed when these matrix elements are expressed as functions of the relative distance 
between the pair of neutrons that are transformed into a pair of protons. We have 



where H(r, E) is the neutrino exchange potential in nucleus and C d u (r) is defined as 



While the matrix element M d u get contributions only from 1 + intermediate states, the 
function C d u {r) gets contributions from all intermediate multipoles. 

The Eq. f II 6 2 j) represents the basic relation between the Qvf3f3- and 2z//3/3-decay 
modes. An analysis of this relation allowed to explained the contrasting behavior of 
Mqj. and M d v when A and Z is changed, namely that M(Y T changes slowly and smoothly 
unlike M d u , which has pronounced shell effects |327j . 

In [52S] a connection of the Fermi Oi//3/3-decay NME Mf with an energy-weighted 
double Fermi transition matrix element was presented. It is arg ued that M F V can be 
reconstructed, if the isospin-forbidden Fermi transition between the ground state of the 
final nucleus and the isobaric analog state in the intermediate nucleus can be measured, 
e.g. by means of (n,p) charge-exchange reactions. By knowing Mp u one can evaluate 
M 0u by assuming an approximate relation Mp u /Mq^ —2.5, which follows from the 
QRPA calculations [173] . 

11. Distinguishing the (V/3/3-decay mechanisms 

Many extensions of the SM generate Majorana neutrino masses and offer a plethora 
of Oi//3/3-decay mechanisms. Among these we should mention the exchange of heavy 
neutrinos, the exchange of SUSY superpartners with R-parity violation, leptoquarks, 
right-handed W bosons, or Kaluza-Klein excitations, among others, which have been 
discussed in the previous sections or can be found in the literature |68j. 

An unambiguous detection of 0i//9/9-decay will prove that the total lepton number 
is broken in nature and neutrinos are Majorana particles. However, after neutrino 
oscillations have established that the neutrinos are massive, as we have already 
mentioned, the observation of Oz//5/3-decay is expected to play a crucial role in 
determining the neutrino mass scale. This prospect generates the questions: What 
is the mechanism that triggers the decay? What happens if several mechanisms are 
active for the decay? 

11.1. Dominance of a single mechanism 

Usually, the Oz/^/S-decay is discussed by assuming that one mechanism at a time 
dominates. Then the half-life in a given nucleus i = (A, Z) can be written as 



C°J T (r)=H GT (r,E)C 2 J(r), 



(162) 




(163) 




(164) 



80 



Here, rj K , M' ^ , G° K U (A,Z) are the LNV parameter (k denotes a given mechanism of 
the Oz//3/3-decay), associated NME and kinematical factor, respectively. The calculation 
of M' ^ = (g^ /gA) 2 M® u (M® u in some cases depends also on g e /^) allows to deduce 
constraint on r] K from the measured lower bound on the O^/3/9-decay half-life. The 
definition of M' ^ in (11641) allows to display the effects of uncertainties in g e /* and to 
use the same phase factor G° K U when calculating the O^/3/9-decay rate [2911 1292]. 

In connection with the neutrino oscillations much attention is attracted to the light 
neutrino mass mechanism of the O^/3/3-decay (r] v = (m u )/m e ) (see section H]). Small 
neutrino masses and neutrino mixing are commonly considered as a signature of physics 
beyond the SM. Several beyond the SM mechanisms of neutrino-mass generation were 
proposed. The most viable and plausible mechanism is the famous see-saw mechanism 
which is based on the assumption that the total lepton number L is violated at a scale 
much larger than the electroweak scale. 

The Oz//3/3-decay is ruled by the light Majorana neutrino-mass mechanism in the 
case of the standard see-saw mechanism of neutrino-mass generation, which is based on 
the assumption that the lepton number is violated at a large (10 15 GeV) scale. In |170] 
it was shown that if O^/3/3-decay will be observed in future experiments sensitive to the 
effective Majorana mass in the inverted mass hierarchy region, then a comparison of 
the derived ranges with measured half-lives will allow us to probe the standard see-saw 
mechanism (see section [3]), assuming that future cosmological data will establish the 
sum of the neutrino masses to be about 0.2 eV . 

A primary purpose of type I see-saw, see section El which is the simplest extension 
of the SM, is to account for light neutrino masses in a renormalizable gauge model. Only 
heavy sterile neutrino states are added to the spectrum of the SU (3)c x SU {2) L x U(l)y 
theory. These heavy states might lead to measurable effects also for the O^/3/3-decay. 
The possible contribution of sterile neutrino dominated Majorana mass eigenstate 
with mass to the O^/3/3-decay was examined in |329] . From the most stringent lower 
bound on the O^/3/3-decay half-life of 76 Ge upper limits on the neutrino mixing matrix 
element \U e h\ 2 in wide region of values of m h (below and above TeV scale) were derived. 
It was assumed that the value of | (m^) | is significantly smaller than the current limit 
on this quantity (\(m u )\ < 0.2 - 0.3 eV). 

Recently, the 01/,9/3-decay associated with the exchange of virtual sterile neutrinos, 
that mix with ordinary neutrinos and are heavier than 200 MeV 0, was revisited |330] . 
The question of having a dominant heavy sterile neutrino contribution in Oz//3/3-decay 
was explored in detail. Due to the improved result of the NMEs |3311 1191] . the bounds 
on active-sterile mixing coming from 0i//9/9-decay has become one order of magnitude 
stronger. The possibility that the sterile neutrino contribution become dominant over 
the light neutrino contribution was addressed for the two flavor and the three flavor 
scenarios |330j . The dominant sterile neutrino contribution in Ou/3/3 process provide 
a way to overcome the conflict between cosmology and the claim for evidence of the 

+ The case of light sterile neutrinos has already been discussed in section S) 
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Oi//9/3-decay by Klapdor and collaborators [JTJH2]. 

There is a possibility that the total lepton number is violated at TeV scale 
[332, 1333] 1334] 1335] , which is accessible at the Large Hadron Collider. The Large Hadron 
Collider can determine the right-handed neutrino masses and mixings. In [333] it was 
manifested that the discovery of left-right (LR) symmetry at the Large Hadron Collider 
would provide a strong motivation for 0^/3/3 searches. By exploiting the LR model with 
type-II see-saw (see section [3]) it was shown that the exchange of heavy neutrinos may 
dominate the Oz//3/3-decay rate depending on the mass of right-handed charged gauge 
boson and the mixing of right-handed neutrinos [333J (see Eq. (|67|) ). A complementary 
study of lepton- flavor violating processes (e.g., /i — > ej), which can provide constraints 
on masses of right-handed neutrinos and doubly charged scalars is of great importance 
[33311336] . 

The LR symmetric models [T71 [TEl [337J are popular models of particle physics due 
to restoration of parity at high energy scale and because they can naturally account for 
the smallness of neutrino masses. They allow not only the light and heavy neutrino mass 
mechanisms of the O^/3/3-decay but also those associated with effective neutrino mass 
independent parameters (rj) and (A), see section [5721 As it was showed already before 
there is an exchange of light neutrinos between two /3-decaying nucleons in nucleus in 
this case. 

The three terms (m u ), (rj) and (A) in the 0i//9/5-decay rate show different 
characteristics in the angular correlations and energy spectrum [63]. By knowing the 
single electron energy spectrum and the angular correlation of the two electrons with 
sufficient accuracy, one could distinguish between decays due to coupling to the left 
handed and right-handed hadronic currents [63J. This possibility was studied in the 
context of 76 Ge and SuperNEMO (Isotopes under consideration for SuperNEMO are 
82 Se, 150 Nd and 48 Ca pHU [338]) detectors. In [3391 EUl USD IMS] the expected pulse 
shapes to be observed for the Ozz/S/S-decay events in a big 76 Ge detector have been 
calculated starting from their Monte Carlo calculated time history and spatial energy 
distribution. The conclusion was that with the spatial resolution of a large size Ge 
detector for the majority of 0^/3/3 events it is not possible to differentiate between 
the contributions of (raj) and the right-handed weak current parameters (rj) and (A). 
Contrary, the SuperNEMO experiment [180] 1338] has a unique potential to measure 
the decay electron's angular and energy distributions and thus to disentangle these 
possible mechanisms for 0v(3(3 decay [18QJ 1343] . We note that the planned experiment 
SuperNEMO, which allows the measurement of Oz//3/3-decay in several isotopes to both 
the ground and excited states, is able to track the trajectories of the emitted electrons 
and determine their individual energies. We should mention here that a measurement 
to both states can also be useful in reducing the background [344] . Other planned 
experiments that will be able to measure the energy and angular distributions are EXO 
[T77] . MOON [315] and COBRA [33B] . 

There is a motivation to consider the Oz//3/3-decay rate in a general framework, 
parameterizing the new physics contributions in terms of all effective low-energy currents 
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allowed by Lorentz-invariance |2011 1347] . e.g. within the effective field theory [348J. 
This approach allows to separate the nuclear physics part of the (V/3/3-decay from 
the underlying particle physics model, and derive limits on arbitrary lepton number 
violating theories. A general Lorentz invariant effective Lagragian for leptonic and 
hadronic charged weak currents was used to perform a comparative analysis of various 
0^/3/9-decay long-range mechanisms in |3491 [350J. It was shown that by measuring of 
angular correlations of emitted electrons in the 0i//9/3-decay together with the ability 
of observing these decays in several nuclei, would help significantly in identifying the 
dominant mechanism underlying this process. 

There is a class of 0i//9/5-decay mechanisms, which one cannot distinguish from each 
other kinematically. The light (r] u ) and heavy (rj^, rj§) Majorana neutrino mass, the 
trilinear R-parity breaking mechanisms - both the short-range mechanism (rjy) with the 
exchange of heavy superpartners (gluino and squarks and/or neutralinos and selectron) 
[27] l88"t l89l I351] and the long-range mechanism (rjg) involving both the exchange of heavy 
squarks and light neutrino |202] (called squark-neutrino mechanism), constitute such a 
group. A discussed possibility to distinguish between these mechanisms is a comparison 
of results for O^/3/3-decay in two or more isotopes [352[ I353[ [354J. 

Under the assumption of a dominance of a single mechanism of the Oz//3/3-decay the 
LNV parameter r] K drops out in the ratio of experimentally determined half-lives for two 
different isotopes. This ratio depend on the mechanism of the O^/3/3-decay due nuclear 
matrix elements and kinematical factors, but is free of LNV parameter. Thus it can 
be compared with the theoretical prediction for different mechanisms. In addition, it is 
assumed that in ratio of nuclear matrix elements theoretical uncertainties are reduced 
due to cancellations of systematic effects. Relative deviations of half-life ratios for 
various new physics contributions, which were normalized to the half-life of 76 Ge and 
compared to the ratio in the light neutrino mass mechanism, were studied in [354J. It 
was found that the change in ratios of half-lives varies from 60% for supersymmetric 
models up to a factor of 520 for extra-dimensional and LR-symmetric mechanisms. It 
is concluded that complementary measurements in different isotopes would be strongly 
encouraged |354l I355] . 

Another possibility to distinguish between the various Oz//3/3-decay mechanisms is 
a study of the branching ratios of Oz//3/3-decays to excited + |352[ I356j and 2 + [&U I357] 
states and a comparative study of the 0z//3/3-decay and neutrinoless electron capture 
with emission of positron (OuEC f3 + ) [358J. Unfortunately, the search for the 0i>EC(3 + - 
decay is complicated due to small rates and the experimental challenge to observe the 
produced x rays or Auger electrons, and most double beta experiments of the next 
generation are not sensitive to electron tracks or transitions to excited states. 

11.2. Nuclear matrix elements of exotic mechanisms 

Recently, the most interest was paid to the calculation of NMEs associated with the light 
neutrino mass mechanism and ground state to ground state transition. Less progress 
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was achieved in the calculation of NMES of exotic O^/9/3-decay mechanisms. 

Experimental studies of transitions to an excited Of and 2f final states □ allow us to 
reduce the background by gamma-electron coincidences. Drawbacks are lower Q values 
and suppressed nuclear matrix elements. The theoretical studies of the corresponding 
nuclear transitions were performed within the ISM [236] , HartreeFockBogoliubov {35 7] 
and QRPA [3211 [3561 EEHJ [360] approaches. In the ISM the Oz//3/3-decays of 48 Ca, 76 Ge, 
82 Se, 124 Sn, 130 Te, and 136 Xe to Of excited final state were found at least 25 times more 
suppressed with respect to the ground state to ground state transition in the case of light 
neutrino mass mechanism. A similar conclusion was found for the 0^/3/3-decays of 76 Ge, 
82 Se, 100 Mo and 136 Xe to the excited collective + state suggesting a suppression 10-100 
larger than that of the transition to ground state [3241 13561 1359] . In addition to light 
neutrino mass also right-handed current [3521 1360] and R-parity breaking mechanisms 
[356] were considered. 

Quite the opposite is claimed in a different study [360j, namely it was found that 
the transition rate of the 0i//9/3-decay of 76 Zr to first excited + state is favored by the 
enhanced transition matrix elements attributed to the monopole-vibrational structure 
of this state. 

The O^/3/3-decay of 76 Ge and 100 Mo to 2f final state was investigated for light 
neutrino mass and right-handed current mechanisms by taking into account recoil 
corrections to the nuclear currents in |357j . The initial + and final 2f nuclear states 
were described in terms of the Hartree-Fock-Bogoliubov type wave functions, which 
were obtained by a variation after particle-number and angular-momentum projection 
[361] . By the numerical calculation of relevant NMEs, it was found that the relative 
sensitivities of + — > 2 + decays to (m u ) and (77) are comparable to those of + — > + 
decays. At the same time it was noted that the + — > 2f decay is relatively more sensitive 
to (A). We should remind the reader that the observation of + — > 2 + transition does 
not establish the presence of right handed currents. For a more complete analysis one 
should consider not only the right handed currents, but supersymmetric contribution as 
well (see the comment at the end of section |6|) 

The right-handed current mechanisms are associated with many different NMEs. 
Within the nuclear shell model they were evaluated just for Oz//3/3-decay of A8 Ca to 
the final ground state [73J. The VAMPIR approach was exploited to calculate them 
in the case of 76 Ge [361] . Many calculations of NMEs related to right-handed current 
mechanisms were performed within the QRPA and for all nuclei of experimental interest 
[3591, 13601 13621 1363[ 1190] . However, they do not include recent improvements concerning 
the fixing of parameter space of nuclear Hamiltonian [2911 12921 12931 13113] and concerning 
the description of two-nucleon short-range correlations [174] . 

There is a revived interest to heavy neutrino mass (t]^' r ) and R-parity breaking 
supersymmetric (r]y, f]q) mechanisms of the 0z//3/9-decay. The NMEs governing these 

* As we have already mentioned, transitions to non zero angular momentum final states can only occur 
via the leptonic jL,jn interference term associated with the (A) and (77) parameters. In this section we 
will refer to it as right handed current contribution. 
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Table 4: Nuclear matrix elements M'J (light Majoran neutrino mass mechanism), M' ^ 
(heavy Majorana neutrino mass mechanism), M'°y (trilinear R-parity breaking SUSY 
mechanism) and M' ~ (squark mixing mechanism) for the 0^/3/3-decays of 76 Ge, 82 Se, 
100 Mo, 130 Te and 136 Xe within the Selfconsistent Renormalized Quasiparticle Random 
Phase Approximation (SRQRPA). R = 1.1 fm A 1 / 3 is assumed. 
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mechanisms were calculated only within the QRPA [2811101] with exception of the PHFB 
calculation for the heavy neutrino mass mechanism [364J, which, however, neglects the 
role of induced hadron currents. 

Recently, nuclear matrix elements M'°" (light neutrino mass mechanism), M' ^ 
(heavy neutrino mass mechanism) , M'°y (trilinear R-parity breaking SUSY mechanism) 
and M' ~ (squark-neutrino mechanism) were calculated for the Oi//3/3-decay of 76 Ge, 82 Se, 
100 Mo, 130 Te and 136 Xe within the SRQRPA [I9H I53T]. Unlike in previous calculations 
the particle-particle strength was adjusted to the 2i/ / 9 / 9-decay half-life and the two- 
nucleon short-range correlations derived from same potential as residual interactions, 
namely from the CD-Bonn and Argonne potentials [174] . were considered. These 
refinements affect mainly heavy neutrino mass NMEs, which became significantly larger. 
In the case of NMEs related to LNV parameter rjsusY the finite nucleon size effect was 
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Table 5: Upper bounds on the lepton number violating parameters (m u ), i]^f R , rjy 
and rjq deduced from current lower bounds on the half-life (T^J^ exp ) of Oi//3/3-decay 
for 76 Ge [36], 82 Se, 100 Mo [38], 130 Te [39] and 136 Xe 00]. Limits on % SUSY 
coupling A' m and on the products of the trilinear ]jip -couplings A' llfc A' lfcl (k=l,2,3) 
for Asusy = 100 GeV are given by assuming that the gluino and squark masses and 
the trilinear soft SUSY breaking parameters are approximately equal to a common 
SUSY breaking scale Asusy — 100 GeV. Nuclear matrix elements calculated within the 
Selfconsistent Renormalized Quasiparticle Random Phase Approximation (CD-Bonn 
potential, qa = 1.25, see Table Hj) are considered. 
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taken into account, which plays an important role in the case of one-pion exchange. 
For large model space and quenched and unquenched value of weak coupling constant 
NMEs of these four mechanisms are presented in Table HI We note that a large model 
space is important to describe reliably especially tensor matrix element contribution to 
the full matrix element. 

The lepton number violating parameters (m„), f]^ R , r]\i and T]q deduced from current 
lower bounds on the half-life (Ty~ exp ) of Oi//3/3-decay for 76 Ge [36], 82 Se, 100 Mo [38], 130 Te 
[39] and 136 Xe [10] are shown in Table El The SRQRPA NMEs of HI in particular those 
evaluated with CD-Bonn potential and qa = 1.25, were considered. We see that upper 
limits on \(m v )\ and \rjg\ from CUORICINO ( 130 Te) [39] and KamLAND-Zen ( 136 Xe) 
[40] experiments are already comparable with those from the Heidelberg-Moscow ( 76 Ge) 
experiment [36] . The running KamLAND-Zen experiment is even slightly more sensitive 
to the 0^/9/3-decay signal as already finished Heidelberg-Moscow experiment in the case 
of the gluino exchange mechanism. 

The rjx' and T]g parameters are related with the with the Ift p -coupling A' m and 
products of the trilinear .^-couplings A 11A ,A lfcl (k=l,2,3), respectively. The current limits 
on them, presented in Table El have been derived under the conventional simplifying 
assumptions. We assumed all the squark masses and the trilinear soft SUSY breaking 
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parameters to be approximately equal to a common SUSY breaking scale Asusy- 
Thus we approximately have [202J 

w w ^ 1 1 / A S USY \ 3 (1KK\ 

AllfcAlfcl - Ek / T oj-e XPG01 | M /0,| UoOGevJ 1 j 



1/2 

with tfc = (1.8 x 10 3 ; 94.2; 3.9) calculated for the current quark masses = 9 MeV, 
m s = 175 MeV and rrib = 4.2 GeV. In the case of gluino and neutralino 1/ip SUSY 
mechanisms of the 0i//9/9-decay we obtain 

11/ rriq \ 2 ( m- g \V 2 
' LU ~ LS /Ty- exp G 01 \M'y\ VlOOGeV J VlOOGeV ) 

' l!J ^ ' ijQv-expQM \ M f^\ VlOOGeV/ VlOOGeV; 1 ; 



1/2 

with nig ~ ~ ~ m x = A$usy- m q, m g, m e and m x are masses of squark, 
gluino, selectron and neutralino, respectively. This approximation is well motivated by 
the constraints from the flavor changing neutral currents. 

It goes without saying that the calculated NMEs of light neutrino mass and exotic 
mechanisms of the O^/3/3-decay depend on the assumption about the nuclear model. In 
order to improve their reliability and reliability of the upper limits on the Oz//3/3-decay 
lepton number violating parameters, further investigations are necessary. 



11.3. Two or more competing mechanisms 

There is a general consensus that a measurement of the 0i//9/9-decay in one isotope 
does not allow us to determine the underlying physics mechanism. Complementary 
measurements in different isotopes is very important especially for the case there are 
competing mechanisms of the Oz//3/3-decay. 

In the case of coexisting mechanisms with identical phase space factors, the Eq. 
fll64p is generalized as 

(T^y 1 = G°»(z) £ I^M^)! 2 . (167) 

K 

Here, g e / = 9a is assumed. The parameters r\ K may take either sign leading to 
constructive or destructive interference in the decay amplitude, if CP conservation is 
assumed. In general case of CP violation they include complex phases. By exploiting 
the fact that the associated nuclear matrix elements are target dependent, given definite 
experimental results on a sufficient number of targets, in principle one can determine or 
sufficiently constrain all LNV parameters including the light neutrino mass term. 

In |331] up to four coexisting mechanisms for the 0z//3/9-decay, mediated by light 
Majorana neutrino exchange (rj v ), heavy Majorana neutrino exchange [t]^), R-parity 
breaking supersymmetry {f]\')), and squark-neutrino (f]q) were considered. Both, 
constructive or destructive interference in the decay amplitude and the Oz//3/3-decay 
in four different candidate nuclei ( 76 Ge, 82 Se, 100 Mo, 130 Te) with NMEs given in Table 
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Figure 25: (Color online) The effective Majorana mass of neutrinos in the case of 
two active mechanisms of the 0t//3/3-decay, namely light and heavy neutrino exchange 
mechanisms, as function of parameter £ (see Eq. ( 11691) ). Here Ty 2 ( 76 Ge) = 2.23lo; 31 x 
10 25 y |H1 H2] is assumed. Solutions were obtained for equal and opposite signs (++, 
+-) on the left hand side of Eqs. (j!68p . The bold point indicates the value of \(m u )\, if 
the light neutrino exchange is the only active mechanism. The dashed regions showed 
the uncertainty of the obtained predictions for if a 3a experimental error of the 

measured half-lives is considered. We also see that a value of | (m„) | > 2.2 eV is excluded 
due to Mainz tritium /3-decay experiment [166J. 



HJwere assumed. It was found that unfortunately, current NME uncertainties appear to 
prevent a robust determination of the relative contribution of each mechanism to the 
decay amplitude, even assuming accurate measurements of Oi//3/3-decay lifetimes. 

Another important feature in analysis of two or more competing mechanisms was 
pointed out in }365j . For example, in the case of two active mechanisms represented by 
the LNV parameters rj u = (m u ) /m e and rj^, assuming the measurement of the 0^/3/3- 
lifetime of two isotopes (i = 76 Ge, 130 Te) and CP conservation, one obtains four sets 
(phases ++, H — , — h, ) of two linear equations: 



Ty 2 (i) G°"(i) ™ e 



AC(i) + vm^ii). (168) 



It was found that this improved analysis leads to completely different results compared 
to those of one mechanism at a time. By making additional assumption that the 0^/3/3- 
decay of 76 Ge was measured with half-life given in [HI H2] the two different solutions 
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for \(m u )\ are plotted as function of £, where 



|M^( 130 Te) | ^T^poTe) G -( 130 Te) 
^~ |M^( 76 Ge)| v /T 1 °; 2 ( 76 Ge) G°"( 76 Ge) ' (169) 



in Fig. [25j The parameter £ represents the unknown half-life of the Oz//3/3-decay of 130 Te. 
We note that for £ = 1 the solution for active only light neutrino mass mechanism is 
reproduced and that £ = means non-observation of the 0^/9/9-decay for a considered 
isotope. By glancing the Fig. [25] the obtained results allows to conclude: 

i) One of the solutions leads to small values of \(m u )\, when all mechanisms add up 
coherently. This is compatible also with inverted (m, < 50 meV) or normal (m^ ~ few 
meV) hierarchy of neutrino masses. 

ii) The second solution allows quite large values of \(m u )\, even larger than 1 eV. It can 
be excluded by cosmology and tritium /3-decay. However, if the claim for evidence will 
be ruled out by running (GERDA |l3l I366[ [367J, EXO |177j ) and future experiments 
[3T1 1368] the values of two solutions will become smaller and perhaps it will not anymore 
be possible to exclude this solution. 

iii) There is possibility that the non-observation of the Oi//3/3-decay for some isotopes 
could be in agreement with a value of | (m u ) | in the sub-eV region. 

iv) The obtained results are sensitive to the accuracy of measured half-lives and to 
uncertainties in calculated nuclear matrix elements. 

Other possibilities of getting information about the different LNV parameters in 
the case of competing O^/3/3-decay mechanisms were discussed in [191] . First, two 
competitive mechanisms, namely light left handed Majorana neutrino exchange and 
heavy right-handed Majorana neutrino exchange, were considered. As the interference 
term is negligible the 0^/9/9-decay half-life for a given isotope is written as 

(l^COG^Ci))" 1 = kM^)| 2 + |r^M°^)| 2 , (170) 

where the index i denotes the isotope. As we have mentioned in subsection 15.21 the 
interference between the left and right handed currents is small. Given a pair of nuclei, 
solutions for \i] u \ 2 and |?7^| 2 can be found by solving a system of two linear equations. 
From the "positivity" conditions (\rj u \ 2 > and \rj§\ 2 > 0) it follows that the ratio of 
half-lives is within the range [19 lj 

G^)|M°T(z)| 2 iy 2 (j) G^)|M°^)| 2 

0»(j)\M$U)\ 3 - " G^{j)\M^{j)f 1 ' 

Surprisingly, the physical solutions are possible only if the ratio of the half-lives, in 
particular of three considered isotopes 76 Ge, 100 Mo and 130 Te, takes values in very narrow 
intervals [191] . 

The 0^/9/9-decay can be triggered also by two competitive mechanisms whose 
interference contribution to the decay rates is non-negligible. As an example the light 
Majorana neutrino mass and gluino exchange mechanisms were considered in [191] . We 
have 

(r^G^ii))' 1 = \v.M?(i)\* + \ Vx/ M°r(i)\* 
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+ 2cos a |^||77 A ,||Mf(z)||M A °r(z)|. (172) 

Here, a is the relative phase of r\ v and r/y. From (j!72p it is possible to extract the values 
of \t) v \ 2 and \rj\>\ 2 and cos a setting up a system of three equations with these three 
unknowns using as input the data on the half-lives of three different nuclei. Results 
of Ref. |191] show that by using of prospective upper bounds on the absolute scale 
of neutrino masses stringent constraints on some of new physics mechanisms, which 
interfere destructively with light neutrino mass mechanism, can be found or even these 
scenarios can be excluded. 

12. Resonant neutrinoless double electron capture 

As it has already been mentioned in section EJ the resonant Oz/ECEC, was considered by 
Winter [50] already in 1955 as a process that would demonstrate the Majorana nature 
of neutrinos and the violation of the total lepton number. The two asterisks denote the 
possibility of leaving the system in an excited nuclear and/or atomic state. The energy 
excess given by the Q-value of the initial atom is carried away by emission of x rays 
(or Auger electrons) as the daughter atom has two electron holes and by emission of a 
single or few photons due to de-excitation of final nucleus. 

The possibility of a resonant enhancement of the Oz/ECEC in case of a mass 
degeneracy between the initial and final atoms was pointed out by Bernabeu, De Rujula, 
and Jarlskog as well as by Vergados about 30 years ago [HH [52] . The half- life of the 
process was estimated by considering non-relativistic atomic wave functions at nuclear 
origin, simplified evaluation of corresponding NME and assuming that the degeneracy 
parameter A = Ma,z — M^ z _ 2 , being the difference of masses of the initial and final 
excited atoms with masses M^z and M™ z __ 2 , varies from zero to 10 keV (representing 
the accuracy of atomic mass measurement at that time). The range of A induced 
uncertainty of about 5 orders in magnitude in calculated OuECEC half-life. A list 
of promising isotopes based on the degeneracy requirement associated with arbitrary 
nuclear excitation was presented. The U2 Sn — > 112 Cd resonant Oz/ECEC transition was 
identified as a good case. 

In 2004 Sujkowski and Wycech [69] and Lukaszuk et al. |369j analyzed the resonant 
Oz/ECEC process for nuclear + — > + transitions accompanied by emission of a single 
photon. By assuming |(?n„)| = 1 eV and 1 a error in the atomic mass determination 
the resonant Oz/ECEC rates of six selected isotopes were calculated by considering the 
perturbation theory approach. The lowest Oz/ECEC half-life was found for 152 Gd. 

The main limitation in identifying promising isotopes for experimental search of 
Oz/ECEC has been poor experimental accuracy of measurement of Q- values which until 
recently were known with uncertainties of 1 - 10 keV only [370J. The resonance 
enhancement can increase the probability of capture by many orders of magnitude. 
Therefore, accurate mass difference measurements are of great importance in order to 
narrow down the possibilities. Progress in precision measurement of atomic masses with 
Penning traps [56J, I371[ 1372] has revived the interest in the old idea on the resonance 
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Oz/ECEC capture. Recently, the accuracy of Q-values at around 100 eV was achieved 
[551 EZ31 EH EZ3 EEH EZ3 EIS EZ3 ESS EHB EH3 EBB], which has already allowed to 
exclude some of isotopes from the list of the most promising candidates (e.g., 112 Sn and 
164 Er) for searching the Oz/ECEC. 

Recently, a significant progress has been achieved also in theoretical description of 
the resonant Oz/ECEC [54"l 1701 1384] . A new theoretical framework for the calculation 
of resonant Oz/ECEC transitions, namely the oscillation of stable and quasi-stationary 
atoms due to weak interaction with violation of the total lepton number and parity, was 
proposed in [M} 170] . The Oz/ECEC transition rate near the resonance is of Breit-Wigner 
form, 

■pOi/ECEC/ ttt\ |Ki&(J r ' r )| -p fl>7Q\ 
L ab [ J ) - a 2 I 1 "P2 af " y 1 ' 6 ) 

^ "+" 4 1 ab 

where J 71 denotes angular momentum and parity of final nucleus. The degeneracy 
parameter can be expressed as A = Q — B ab — E T Q stands for a difference between 
the initial and final atomic masses in ground states and E 1 is an excitation energy of 
the daughter nucleus. B ab = E a + E b + E c is the energy of two electron holes, whose 
quantum numbers (n,j,l) are denoted by indices a and b and Ec is the interaction 
energy of the two holes. The binding energies of single electron holes E a are known with 
accuracy with few eV [385J . The width of the excited final atom with the electron holes 
is given by 

r a6 = r + r 6 + r*. (m) 

Here, T atb is one-hole atomic width and T* is the de-excitation width of daughter nucleus, 
which can be neglected. Numerical values of T ab are about up to few tens eV [386J . 

For light neutrino mass mechanism and favorable cases of a capture of s\/2 and 
P1/2 electrons the explicit form of lepton number violating amplitude associated with 
nuclear transitions + —¥ J n = ±l , l 1 *" 1 is given in [54]. By factorizing the electron shell 
structure and nuclear matrix element one get 

V ab (r) = ]- <%m#A < E ab > M^ ECEC {r). (175) 

47T K 

Here, < F ab > is a combination of averaged upper and lower bispinor components of the 
atomic electron wave functions [54 J and M 0uECEC (J 7T ) is the nuclear matrix element. We 
note that by neglecting the lower bispinor components M 0uECEC (0 + ) takes the form of 
the Oz//3/3-decay NME for ground state to ground state transition after replacing isospin 
operators r~ by r + . 

There is a straightforward generalization of the LNV potential V ab (0 + ) in (I175P 
for the heavy neutrino exchange, the trilinear R-parity breaking with gluino and 
neutralino exchange and squark-neutrino mechanisms. It is achieved by replacements 
rj v = (m u )/m e with r] K (k = N, A', q) and M^ ECEC (0+) with M° uECEC (0 + ). The 
Oz/ECEC leading to final states different than + , possible only in the presence weak 
right-handed currents due to the leptonic of jh ~ 3r interference, has been discussed in 
[381]. 
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New important theoretical findings with respect of the Oz/ECEC were achieved in 
[54J. They are as follows: i) Effects associated with the relativistic structure of the 
electron shells reduce the Oz/ECEC half-lives by almost one order of magnitude, ii) The 
capture of electrons from the npi/ 2 states is only moderately suppressed in comparison 
with the capture from the nsx/2 states unlike in the non-relativist ic theory, iii) For 
light neutrino mass mechanism selection rules appear to require that nuclear transitions 
with a change in the nuclear spin J > 2 are strongly suppressed. We note that, if 
right-handed currents are considered, selection rules are modified allowing also J ^ + 
[384J. iv) New transitions due to the violation of parity in the Oz/ECEC process were 
proposed. For example, nuclear transitions + — > O 1 * 1 , l 1 * 1 are compatible with a mixed 
capture of s- and p-wave electrons, v) The interaction energy of the two holes Eq has 
to be taken into account by evaluating a mass degeneracy of initial and final atoms, vi) 
Based on the most recent atomic and nuclear data and by assuming M 0uECEC '(J 7T ) = 6 
the Oz/ECEC half-lives were evaluated and the complete list of the most perspective 
isotopes for further experimental study was provided. Some isotopes such as 156 Dy have 
several closely-lying resonance levels. A more accurate measurement of Q- value of 156 Dy 
by Heidelberg group confirmed the existence of multiple-resonance phenomenon for this 
isotope [380]. vii) In the unitary limit some Oz/ECEC half-lives were predicted to be 
significantly below the Oz//3/3-decay half-lives for the same value of (m u ). A probability of 
finding resonant transition with low Oz/ECEC half-life was evaluated, vii) The process of 
the resonant neutrinoless double electron production (Oz/EPEP), i.e. neutrinoless double 
beta decay to two bound electrons, namely 

(A, Z) ^(A,Z + 2)** + e 6 " + e 6 ", (176) 

was proposed and analyzed. This process was found to be unlikely as it requires that 
a Q-value is extremely fine tuned to a nuclear excitation. The two electrons must be 
placed into any of the upper most non-occupied electron shells of the final atom leaving 
only restricted possibility to match to a resonance condition. 

A detailed calculation of the Oz/ECEC of 152 Gd, 164 Er and 180 W associated with 
the ground-state to ground-state nuclear transitions was performed in [551 1387[ [388J. 
Improved measurements of Q-value for these transitions with accuracy of about 100 eV 
[381] 13831 1382[ 1387] were considered. The nuclear matrix elements of 152 Gd — > 152 Sm, 
164 Er — > 164 Dy and 180 W — > 180 Hf transitions were calculated within spherical and 
deformed QRPA [3871 [388J . The obtained results excludes 164 Er and 180 W from the list 
of prospective candidates to search for the Oz/ECEC. The Oz/ECEC half-life of 152 Gd is 
2-3 orders of magnitude longer than the half-life of 0z//3/3 decay of 76 Ge corresponding 
to the same value of (m u ) and is the smallest known half-life among known Oz/ECEC 
transitions at present. 

The transition of 106 Cd to an excited state of 106 Pd with the nuclear excitation 
energy of 2717.59 keV was calculated in Ref. [389] by making assumption that this is 
the + state. However, it was noted in [54] that, as long as this level 7-decays by 100% 
into the 3 + state at 1557.68 keV, this possibility is excluded. 
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Table 6: A comparison of the neutrinoless double beta decay and the resonant 
neutrinoless double electron capture for light neutrino mass mechanism. The lepton 
number violating amplitude V a b{J n ) is given in Eq. (11751) . 



0^/3/3-decay 



OvECEC 



definition 
formalism 
half-life 

nucl. trans, 
fav. at. syst. 



uncert. in T® u 



/2 



exp. sign. 



rpOu—exp 
1 l/2 

exp. act. 



background 



(A,Z) -> (A,Z + 2) + e~ + e 
perturbation field theory 



rp[)u 

1/2 



' GV\M° V (J*] 



0+ -> 0+,2+ 
large Q-value (3-4 MeV) 



48 Ca, 76 Ge, 76 Se, 100 Mo, 

116 Cdj 130 Te; 136 Xe; 150 Nd; 

factor ~ 4 — 9 due 
to calc. of NME 

peak at end of sum 
of two el. energy spectra 

> 10 24 -10 25 y 
const, of (0.1-1 ton) exp. 
with sensitivity to inverted 
hierarchy of neutrino masses 
2i//3/3-decay 
upon resolution of exp. 



(A, Z)+ e - + e - -+(A,Z-2)** 
oscillation of atoms [701 ET 



In 2 IVqbCJ^)! 2 
TfffVEC {M A ,z 



»,/-** \2 T lp2 



A,Z-2 

o+ ->■ o+,o-,i + , 1" 

mass difference 
of few tens of eV 
unknown yet ( 106 Cd, 124 Xe, 
152 Gd, 156 Dy 168 Yb,... [M]) 

many orders in magn. 
up to measured mass diff. 
and due to NMEs 
x rays or Auger el. 
plus nucl. de-excitation 
> 10 19 -10 20 y 
small exper. yet 



2vECEC is strongly 
suppressed 
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There is also an increased experimental activity in the field of the resonant Oz/ECEC 
[7TJ I390J E9J1 13921 E931 EH2 [395] . The resonant Oi/ECEC has some important advantages 
with respect to experimental signatures and background conditions. The de-excitation of 
the final excited nucleus proceeds in most cases through a cascade of easy to detect rays. 
A two- or even higher-fold coincidence setup can cut down any background rate right 
from the beginning, thereby requiring significantly less active or passive shielding |54j. 
A clear detection of these 7 rays would already signal the resonant 0u ECEC without 
any doubt, as there are no background processes feeding those particular nuclear levels. 
It is worth noting that lepton number conserving ECEC with emission of two neutrinos, 

(A, Z) + e- h + e 6 - -> (A, Z - 2)** + v e + u e , (177) 

is strongly suppressed due to almost vanishing phase space [511 [521 [54]. The ground 
state to ground state resonant Oz/ECEC transitions can be detected by monitoring the 
x rays or Auger electrons emitted from excited electron shell of the atom. This can be 
achieved, e.g., by calorimetric measurements. 
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Till now, the most stringent limit on the resonant Oz/ECEC were established for 
74 Se [391], 106 Cd [393] and 112 Sn [391] • The ground state of 74 Se is almost degenerate to 
the second excited state at 1204 keV in the daughter nucleus 74 Ge, which is a 2 + state 
[396]. The 27-ray cascade has been searched for by using the low-radioactivity detector 
setup at the Comenius University in Bratislava and 3 kg of natural selenium. A lower 
limit for the half-life of Tfj^ CEC > 4.3 x 10 19 y was determined |394j . which is slightly 
larger than the value reported in |390j . The resonant Oz/ECEC transition to the O3" 
excited state in 112 Cd (1871.0 keV) has been investigated in an experiment performed 
with natural tin in the Modane Underground Laboratory. A lower bound on half-life of 
0.92 x 10 20 y was established. It is worth noticing that a new mass measurement [376] 
has excluded a complete mass degeneracy for a 112 Sn decay and has therefore disfavored 
significant resonant enhancement of the Oz/ECEC mode for this transition. Within the 
TGV experiment in Modane [393] an interest has also arisen in the Oz/ECEC resonant 
decay mode of 106 Cd (KL-capture) to the excited 2741 keV state of 106 Pd. The spin 
value of this final state was unknown and it was assumed to be J = (1, 2) + . After 
measurements had begun a new value for the spin of the 2741 keV level in 106 Pd of 
J = 4 + was adopted, but , following recent theoretical analysis [54] , this channel is now 
disfavored. Nevertheless the most stringent limit on the Oz/ECEC half-life of of 1.1 x 10 20 
y was reported [393J. 

A comparison of the Oz/ECEC with the 0z//3/3-decay is presented in Table [6j It 
is maintained that these two lepton number violating processes are quite different 
and at different levels of both theoretical and experimental investigation. Precise 
measurements of Q-values between the initial and final atomic states, additional 
spectroscopic information on the excited nuclear states (energy, spin and parity) and 
reliable calculation of corresponding NMEs are highly required to improve predictions 
of half-lives of the resonant Oz/ECEC. It is expected that the accuracy of 10 eV in 
the measurement of atomic masses will be achievable in the near future. The electron 
binding energy depends on the local physical and chemical environment. An interesting 
question is whether it possible and, if so, how to manage the atomic structure in such a 
way as to implement the degeneracy of the atoms and create conditions for the resonant 
enhancement, as discussed in a recent work [54] . 

13. Concluding remarks 

In this review we discussed in some detail the lepton number violating neutrinoless 
double beta decay and other similar transitions, involving various nuclear isotopes for 
which ordinary beta decay and e-capture are forbidden or highly suppressed. Both 
theoretical and experimental aspects were considered. 

We have seen that this is a process with long and interesting history with important 
implications for physics and cosmology, but its observation is still elusive. It is an exotic 
process, which requires physics beyond the SM. At present a complete theory is missing 
and, thus, to motivate and guide the experiments we examined a number of reasonable 
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viable models, beyond the SM, in particular in connection with the neutrino mass matrix 
and mixing (see sections 131 and . Such models predict that lepton number violation, 
and consequently neutrinoless double beta decay, must occur at some level, implying 
that the neutrinos are Majorana particles. These models, however, cannot provide 
a precise determination of the parameters involved, such as the absolute scale of the 
neutrino mass. So they must be extracted from the experiments, if and when reliable 
accurate results become available (see section H] for the neutrino mass). The observed 
values may, then, be used to differentiate between such models and, hopefully, lead to 
the ultimate theory. 

In order to achieve this goal first such processes must be definitely observed. Then 
the obtained results must be analyzed by considering the various mechanisms implied by 
the above models, see sections [5] and [6] for mechanisms involving intermediate neutrinos 
and other particles respectively. This, however, can only be done, if the corresponding 
nuclear matrix elements are evaluated with high precision, accuracy and reliability. We 
have seen that this is a formidable task, since the nuclei that can undergo double beta 
decay have rather complicated structure. 

The evaluation of the nuclear matrix elements involves two steps. In the first step 
the effective transition operators for each mechanism (see section [9]). Special attention 
must be paid to the proper treatment of these operators at short distances (short range 
correlations, nucleon current corrections, inclusion of hadrons other than nucleons etc). 
The second step consists of selecting the proper nuclear model for constructing the 
wave functions involved in the evaluation of the nuclear matrix elements. Practically 
all models available in the nuclear theory artillery have been employed. The most 
prominent are the large basis shell model, the various refinements of the quasi particle 
random phase approximation (QRPA) and the interacting boson model (IBM). The 
essential features of these models and the numerical values of the obtained nuclear 
matrix elements have been summarized in section UM We have seen that great progress 
has been made in this direction in recent years and it is encouraging that the nuclear 
ME obtained with these vastly different nuclear models tend to converge. 

We have discussed in section [8] the ongoing, planned and future experiments. We 
have witnessed great progress in tackling the various background problems, improving 
the energy resolution and preparing large masses of the needed isotopes. It is, thus, 
expected that half lives of the order of 10 26 y can be achieved and, consequently, a 
sensitivity of a few tens of meV for the average neutrino mass can be reached. This 
may be sufficient to differentiate between the normal and inverted hierarchy scenarios 
(see section H]). Furthermore we have seen that various nuclear charge changing nuclear 
reactions can be employed in an effort to experimentally extract useful information or 
provide checks for the nuclear matrix elements. 

It is clear that the observation of neutrinoless double beta decay will be a great 
triumph for physics and experimental physics in particular. It will demonstrate that the 
neutrinos are Majorana particles and there exist lepton number violating interactions 
in the universe. This, however, will not be the end of the story. The data should be 
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analyzed in such a way to determine the mechanism responsible for this process and, in 
particular, to extract the most important parameter, which is scale of the neutrino mass. 
Great progress has been in this direction has recently been made as briefly exposed in 
section [Til In order to unambiguously accomplish this goal, however, the accuracy of 
the nuclear matrix elements must be further improved. 

Finally recent developments, towards the accurate determination of atomic masses 
as well as the evaluation of the inner shell atomic wave functions and energies, have 
stimulated interest in experiments involving the resonant neutrinoless double e-capture, 
see section [121 This new process, if observed, especially in case it leads to negative 
parity final nuclear states, will greatly facilitate the analysis of determining the dominant 
mechanism involved in neutrinoless double beta decay. 

It is clear that theoretical attempts in the determination of the nuclear matrix 
elements and experimental efforts towards achieving the observation of neutrinoless 
double beta decays, involving as many as possible nuclear isotopes and utilizing all 
available techniques, should be encouraged and supported. 
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